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PREFACE. 


IN  this  volume,  an  attempt  has  been  made  to  produce  a 
Treatise  on  the  Planetary  Theory,  which,  being  elementary 
in  character,  should  be  so  far  complete,  as  to  contain  all  that 
is  usually  required  by  students  in  this  University.  But  it 
is  not  without  diffidence  that  I  submit  my  volume  to  their 
notice.  In  the  earlier  part  of  it,  the  methods  which  have 
been  adopted  are  to  some  extent  original*,  and  the  general 
arrangement  of  the  second  chapter  will,  it  is  believed,  be 
found  to  be  new.  Through  the  kindness  of  the  Publishers, 
a  portion  of  Pratt' s  Mechanical  Philosophy  has  been  placed 
at  my  disposal.  Of  this  I  have  availed  myself,  particularly 
in  the  chapter  on  the  Stability  of  the  Planetary  System ;  but, 
on  the  whole,  comparatively  little  has  been  reprinted  ver- 
batim from  that  work.  Among  other  sources  of  information, 
my  obligations  are  mainly  due  to  Ponte*coulant's  TMorie 
Analytique  du  Systeme  du  Monde,  Airy's  Mathematical 
Tracts,  and  Frost's  Planetary  Theory  in  the  Quarterly 
Journal  of  Mathematics:  but  I  have  also  referred  to  the 

*  Some  of  these  have  already  appeared  in  Mathematical  Journals. 
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VI  PREFACE. 

Mecanique  Celeste,  the  Mecanique  Analytique,  Mrs  Somer- 
ville's  Mechanism  of  the  Heavens,  (a  work  forming  a  complete 
Mathematical  treatise  on  Physical  Astronomy,)  a  Memoir 
by  Prof.  Donkin  on  the  Differential  Equations  of  Dynamics, 
Phil.  Trans.  1855,  &c.  A  collection  of  Problems  has  been 
added,  taken  chiefly  from  the  Smith's  Prize  and  Senate- 
House  Examination  Papers  of  the  last  twenty  years.  In 
conclusion,  I  would  express  my  sincere  thanks  to  Messrs 
A.  Freeman,  P.  T.  Main,  and  other  friends,  of  St  John's 
College,  for  the  valuable  assistance  which  they  have  afforded 
me,  and  would  venture  to  hope  that  the  work  will  be  found 
useful. 

C.  H.  H.  CHEYNE. 


ST  JOHN'S  COLLEGE, 

October  9,  1862. 
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THE    PLANETARY    THEORY. 

CHAPTER    I. 

INTRODUCTION. 

1.  To  determine  the  motion  of  a  system  of  bodies  sub- 
ject only  to  their  mutual  attractions,  is  a  problem  the  mathe- 
matical difficulties  of  which   have   not  yet  been  overcome: 
hence  in  the  particular  cases  of  this  problem  which  the  Lunar 
and  Planetary  Theories  present,  recourse  must  be   had   to 
methods  of  approximation.     Happily  the  arrangement  of  the 
Solar  System  renders  approximate  methods  possible,  and  in 
the  skilful  hands  of  the  Mathematicians  of  the  last  century, 
they  have  been  brought  to  a  high  state  of  perfection. 

2.  If  the  Sun  were  the  only  attracting  body,  the  planets 
would  describe  exact  ellipses,  agreeably  to  Kepler's  first  law  ; 
but   in  consequence  of  the  attractions  of  the  planets  them- 
selves, slight  deviations  from  elliptic  motion  are  produced. 
The  method  of  calculating   these  deviations,    to  which  our 
attention  will  chiefly  be  directed,  is  due  to  Euler ;  it  consists 
in  supposing  the  planets  to  move  in  ellipses,  the  elements  (or 
arbitrary  constants)  of  which  are  continually  though  slowly 
changing*. 

*  The  legitimacy  of  this  hypothesis  will  appear  when  we  come  to  treat  of 
the  equations  of  motion.     See  Arts.  21  and  22. 
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2  PLANETARY   THEORY. 

3.  Now  the  elements  of  an  elliptic  orbit  are  (i)  the  mean 
distance,    or  semi-axis  major,    (ii)   the  excentricity ,  (iii)   the 
longitude  of  perihelion,  i.  e.  of  the  point  of  the  orbit  nearest  to  the 
Sun,  (iv)  the  longitude  of  the  epoch*,  or  mean  longitude  at  the 
epoch  from  which  the  time  is  reckoned,  (v)  the  inclination  of 
the  plane  in  which  the  orbit  lies  to  some  fixed  plane  of  refer- 
ence,  (vi)  the  longitude  of  the  ascending  node.    Of  these  (i) 
and  (ii)  determine  the  magnitude  of  the  orbit,  (iii)  determines 
its  position  in  its  own  plane,  (v)  and  (vi)  determine  the  posi- 
tion of  this  plane,  and  (iv)  has  reference  to  the  position  of  the 
body  itself  in  its  orbit. 

If  the  planets  moved  accurately  in  ellipses,  these  would  be 
constants  :  we  must  however  be  prepared  to  consider  them  as 
variable  quantities,  which  it  will  be  the  object  of  the  problem 
to  determine.  They  are  tejmed  the  elements  of  the  orbit. 

4.  But  further,  not  only  is  it  found  that  the  true  orbit  of 
a  planet  is  not  an  ellipse,  but  that  it  is  not  even  a  plane  curve, 
although  the  departure  of  the  planet  from  the  plane  in  which 
it  is  at  any  instant  moving  is  extremely  slow.    We  define  as  the 
plane  of  the  orbit  the  plane  containing  the  radius  vector  and 
direction  of  motion  of  the  planet  at  the  instant  under  con- 
sideration. 

5.  We  shall  suppose  the  Sun  and  planets  so  distant  from 
each  other  that  they  may  be  considered  to  attract  as  if  they 
were  condensed  into  their  respective  centres  of  gravity;   a 
supposition  which  would  be  rigorously  true  if  these  bodies  were 
exactly  spherical,  and  either  of  uniform  density  or  composed  of 
concentric  spherical  shells,  the  density  of  each  shell  being  uni- 
form throughout.     The  errors,  however,  thus  introduced  into 
the  motions  of  translation  are  found  to  be  inappreciable  for  the 
planets,  though  not  in  the  case  of  their  satellites.     The  mo- 
tions of  rotation  will  not  be  considered  in  the  present  treatise. 

*  Also  briefly  termed  the  epoch. 
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6.  Moreover,  since  the  masses  of  the  planets  are  extremely 
small  in  comparison  of  that  of  the  Sun,  it  follows  that  in  cases 
where  it  is  not  necessary  to  carry  the  approximation  beyond 
the  first  order  of  these  masses,  we  are  permitted  to  avail  our- 
selves of  the  Principle  of  the  Superposition  of  Small  Motions, 
and  thus  to  reduce  the  problem  to  a  case  of  that  of  the  Three 
Bodies. 

7.  So  far  the  Theory  of  the  Planets  resembles  that  of  the 
Moon,  and  the  same  method  of  treatment  might  be  employed 
in  both  cases.     But  they  differ  in-  this  respect :  the  ratio  of 
the  distances  of  the  disturbed  and  disturbing  bodies  from  the 
central  one*  is  much  smaller  in  the  Lunar  than  in  the  Planetary 
Theory,  so  that  if  in  the  latter  theory  the  approximation  were 
made  by  means  of  series  proceeding  by  powers  of  this  ratio,  it 
would  be  necessary  to  retain  many  more  terms  than  are  re- 
quired in  the  former.     For  this  reason  a  different  method  of 
development  is  employed.     The  perturbations  of  the  Moon, 
however,  are  far  larger  than  those  of  the  planets,  since  in  the 
former  case  the  Sun,  of  which  the  mass  is  enormous,  and  the 
distance  not  proportionately  great,  is 'one  of  the  disturbing 
bodies. 

8.  To  find  an  expression  for  the  component  in  any  di- 
rection of  the  force  which  disturbs  the  motion  of  a  given  planet 
relatively  to  the  Sun. 

Let  M  denote  the  mass  of  the  Sun,  m,  m',  m",  &c.,  those 
of  the  planets,  and  suppose  the  relative  motion  of  m  re- 
quired. 

Let  x,  y,  z,  x,  y',  z',  x",  y",  z",  &c.,  be  the  co-ordinates 
of  m,  m',  m",  &c.,  referred  to  any  system  of  rectangular  axes 

*  By  the  central  body  is  meant  that  whose  attraction  exercises  the  greatest 
influence  on  the  body  whose  motion  is  required ;  the  Sun  for  instance  in  the 
Theory  of  the  Planets,  and  the  Earth  in  that  of  the  Moon.  All  the  other 
attracting  bodies  are  called  disturbing  bodies. 
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originating  in  the  centre  of  gravity  of  the  Sun ;  r,  r',  r",  &c., 
their  distances  from  the  origin ;  p',  p",  &c.,  the  distances  of 
m,  m",  &c.,  from  m. 


Now  if  to  every  body  of  the  system  we  apply  forces  equal 
and  opposite  to  those  which  act  upon  the  Sun,  we  shall  reduce 
the  latter  to  rest  without  affecting  the  relative  motion.  Hence, 
considering  the  action  of  only  one  disturbing  planet  m',  the 
forces  acting  upon  m  will  be 

M+m.     ,. 

—  2  —  m  direction  mm, 

m  .     ,.       .  , 

-is  m  direction  mm  , 
P 

m  .     ,.        .          ,,, 
-jj  in  direction  m  M, 

of  which  the  last  two  constitute  the  disturbing  force. 

; 
77? 

Let  F=  —  ,  :  then   on  the  hypothesis  of  Art.  5,  Fwill  be 

the  potential  of  m',  and  the  components  parallel  to  the  axes  of 
the  disturbing  force  due  to  the  action  of  m,  will  be 


dV 


mx 
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dV     my 

dy       r«  ' 

dV     mz 

dz       r*  ' 

Let  s  denote  the  length  of  the  arc  of  any  curve  measured 
from  some  fixed  point  up  to  m*  :  then  the  resolved  part  of  the 
disturbing  force  parallel  to  the  tangent  at  m  to  this  curve 
will  be 

dV     /m'x'  dx     m'y1  dy     m'z  dz" 

~ds  ~~  (~^  ~ds~ 

or  restoring  to  V  its  value 
d 


which  may  be  written  -^—  ,  if 


m 


-n,  ,     t         i         f\ 

R  =  j--(xx+yy  +zz). 

If  we  express  in  like  manner  the  disturbing  forces  due 
to  the  actions  of  m",  m",  &c.,  we  shall  have  for  the  whole 
component  in  this  direction 


ds  "  ds       '"' 

or  ^,  if  B  =  R  +  K'  +  .... 
ds 

The  function  E  is  called  the  disturbing  function. 

77-> 

9.     From  the  manner  in  which  -^-  has  been  introduced, 

it  appears  that  E  is  supposed  to  be  expressed  in  terms  of  s 
and  quantities  which  do  not  vary  with  s.     It  must  however 

*  This  arbitrary  curve  we  shall  term  the  carve  of  reference. 
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be  borne  in  mind  that  in  -7-  the  variation  is  purely  hypo- 
thetical, and  has  nothing  whatever  to  do  with  the  actual 
variation  of  R  due  to  the  motion  of  the  planet. 

For  example,  suppose  the  curve  of  reference  a  straight 
line  parallel  to  the  axis  of  x,  and  let  R  be  expressed 
in  terms  of  x,  y  and  z  •  then  in  this  case  x  only  will 
vary,  and  the  disturbing  force  parallel  to  the  axis  of  x  will 

be  denoted  by  -j— ,  y  and  z  being  considered  constant  in  the 
differentiation.  Similarly,  the  disturbing  forces  parallel  to 
the  axes  of  y  and  z  will  be  expressed  by  -7-  and  -7—  respec- 
tively, the  differential  coefficients  being  strictly  partial. 

Again,  suppose  the  curve  of  reference  a  circle  with  its 
plane  parallel  to  that  of  xy,  and  its  centre  in  the  axis  of  z, 
and  let  R  be  expressed  in  terms  of  the  polar  co-ordinates 
(r1?  0J  of  the  projection  of  the  planet  on  the  plane  of  xy,  and 
its  distance  (z)  from  this  plane ;  then  in  this  case  0t  only  will 
vary,  and  the  disturbing  force  perpendicular  to  the  projected 

/77? 

radius  vector  will  be  expressed  by  ,-  ,  r^  and  z  being  con- 
sidered constant  in  the  differentiation.  Similarly,  the  forces 
parallel  to  the  projected  radius  vector  and  to  the  axis  of  z, 

/77? 
will  be  expressed  by  the  partial  differential  coefficients  -=— , 

dR 

~j~  respectively. 

10.  The  disturbing  function,  like  the  potential,  is  inde- 
pendent of  any  particular  system  of  co-ordinates  that  may  be 
employed.  For 

jy      mm.,          ,         ,. 
7  ~       ^XX      yy      ZZ  ^ 
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m      in  r  fx  x      y  y      z  z 
p        r*   \r  r'      r    '      r  r 


m       mr 

— ,.2-  cos  ft>, 

p        r 


if  co  denote  the  inclination  of  r'  to  r. 

11.  To  express  B/  in  terms  of  the  polar  co-ordinates  of 
the  projections  of  m  and  m'  on  a  fixed  plane,  and  of  their 
distances  from  it. 

Take  the  fixed  plane  for  that  of  xy :  let  rv  r/  be  the  pro- 
jections of  r,  r  upon  it,  and  0^  6^  the  inclinations  of  rl9  r^  to 
the  axis  of  x ;  then 

x  =  rx  cos  015      y=r^  sin  6V 

x  =  r/  cos  0/,    y  =  r/  sin  0/  ; 
therefore      cca?'  +  yy  +  ««'  =  r^  cos  (^  —  ^/)  +  zz, 

p'*=(X-X>y+(y-y>Y  +  (Z-z'y 

=  r*  +  r2  -  2'  cos  0  -  0     +      -  « 


Hence  by  substitution, 
R'  = 


,2  +  <2  -  2V/  COS  (^  -  0/)  +  (0  -  ^T} 


12.  In  a  subsequent  chapter  we  shall  consider  the  de- 
velopment of  R  in  terms  of  the  time  and  the  elements  of  the 
orbit,  in  a  series  ascending  by  powers  and  products  of  the 
excentricities  and  inclinations,  which  in  the  Planetary  Theory 
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are  very  small.  At  present  we  shall  content  ourselves  with 
shewing  how  R  may  be  expressed  in  terms  of  these  quantities. 
We  shall  assume  that  the  equations  connecting  the  co-ordi- 
nates, the  time,  and  the  elements  in  an  elliptic  orbit,  hold 
also  in  the  case  of  a  disturbed  planet. 

13.  To  explain  how  R  may  be  expressed  in  terms  of  the 
time  and  the  elements  of  the  orbit. 

Let  r,  0  denote  the  radius  vector  and  longitude  of  the 
disturbed  planet,  the  latter  being  measured  on  a  fixed  plane 
of  reference  as  far  as  the  node,  and  thence  on  the  plane  of  the 
orbit :  let  the  elements  be  a  the  mean  distance,  e  the  excen- 
tricity,  TX  the  longitude  of  perihelion,  e  the  longitude  of  the 
epoch,  (the  last  two  being  measured  in  the  same  way  as  6,) 
fl  the  longitude  of  the  node  measured  on  the  plane  of  refer- 
ence, and  i  the  inclination  of  the  plane  of  the  orbit  to  the 
plane  of  reference.  Our  object  is  to  express  R  in  terms  of  t 
and  these  elements. 

Again,  let  #0,  «TO,  e0,  O0  denote  the  longitudes  of  the 
planet,  of  perihelion,  of  the  epoch,  and  of  the  node,  measured 
entirely  on  the  plane  of  the  orbit. 

Let  a. sphere  be  described  with  its  centre  coinciding  with 
that  of  the  Sun,  and  its  radius  of  any  magnitude:  let  the 
planes  of  reference  and  of  the  orbit  cut  it  in  the  great  circles 
NM,  NP,  then  the  line  of  nodes  will  cut  it  in  N]  let  the 
radius  vector  of  the  planet  cut  it  in  P,  the  projection  of  this 
radius  on  the  plane  of  reference  in  M,  and  the  lines  from 
which  6,  00  are  measured  in  L,  0  respectively.  We  shall 
suppose  L  to  be  the  same  origin  as  that  from  which  0X  is 
measured  in  Art.  11. 

Then  in  the  figure  LM=dv  LN  +  NP  =  0,  OP  =  0Q, 
LN=  ft,  the  angle  PNM=  i,  and  PM=  the  latitude  of  the 
planet  which  we  shall  denote  by  X. 
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Hence  from  the  right-angled  triangle  PNM, 

tan  (0t  -  H)  =  cost  tan  (0-H) (1), 

sin X  =  sin i sin  (0  -  ft) (2): 

also  r1  =  rcosX (3), 

z  =  r  sinX (4). 


Again,  from  the  formulas  of  elliptic  motion*, 
r  =  a{l  +  %ez  -  e  cos  (nt  +  e0  -  OTO)  -  Je2  cos  2  (nt  +  e0  -  -BJO)  -  ...}, 
00=  w<+  e0  +  2e  sin  (w«  +  e0  -  tsr0)  +  Je2  sin  2  (w«  +  e0  -  «TO)  +  .  ..  t  • 
but  0-60  =  LN-  ON=e-e0  =  <ar-'&0', 

therefore  6  =  00  +  e  —  e0,     60  —  OTO  =  e  —  «r, 

and  our  formulae  become 

r  =  a  }  1+  Je2-  e  cos  (nt  +  e  -  w)  -  |e2cos  2  (w^+  e  -  w)  —  .  .  .  }  .  .  .  (5)  , 

-'C7)+  ......  (6). 


In  Art.  11  we  have  expressed  R'  in  terms  of  r1?  61  and  «  ; 
hence  by  equations  (1)  to  (4)  it  may  be  expressed  as  a  func- 
tion of  r,  0,  H,  and  ^  :  we  may  then  substitute  for  r  and  6 
from  equations  (5)  and  (6),  and  R'  will  be  expressed  in  terms 

*  See  Tait  and  Steele's  Dynamics,  Arts.  114  and  122. 

"\"  n  is  termed  the  mean  motion,  and  is  connected  with  the  mean  distance  by 
the  equation  W2a3=/t,  where  p.  = 
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of  t  and  the  elements  of  the  orbit.  Similarly  R",  R",  &c., 
and  therefore  R  may  be  expressed  in  terms  of  t  and  the 
elements. 

14.  We  proceed  to  investigate  certain  relations  which 
subsist  between  the  partial  differential  coefficients  of  R  with 
respect  to  the  co-ordinates  of  the  disturbed  planet,  and  its 
partial  differential  coefficients  with  respect  to  the  elements  of 
the  orbit.     These  will  be  useful  in  obtaining  the  formulae  by 
which  the  values  of  the  elements  are  calculated. 

We  premise  that  when  we  speak  of  the  partial  differ- 
ential coefficient  of  R  with  respect  to  one  of  the  elements,  we 
suppose  R  expressed  in  the  manner  indicated  in  the  last 
article,  and  that  the  time  as  well  as  the  other  elements  are 
considered  constant  in  the  differentiation:  when  we  speak 
of  the  partial  differential  coefficient  of  R  with  respect  to  r  or  #, 
we  suppose  R  expressed  in  terms  of  r,  0,  i  and  fl  ,  which  may 
be  done  by  equations  (1)  to  (4). 

dR     dR  ,  dR 

15.  To  shew  that     -3^-  =  -5-  +  -3—  . 

d0       de      dtzr 

Equations  (5)  and  (6)  of  Art.  13  may  be  written 


r  «**  +  «-.« 


whence  it  follows  that 


dr      dr  _ 
*  ' 


dO     d0^_ 
de      der 


Now  since  e  and  «r  enter  into  R  only  through  r  and 


dRL_dRL  dr     dR  M 
~fa~~dr~  d~e+~d§   de  ' 
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dE  _dR  dr_     dR^de_^ 

d^r      dr  dm      dO  dtz  ' 


therefore,  by  addition, 

de      dix      dO 

OL 


1fi      „,     7        ,       dR     dR     dR     dR 

16.     To  shew  that    -^  =  -v-  +  -\ h  ^7^  • 

do'j       de       dtir      all 

From  equations  (1)  to  (4)  of  Art.  13  we  obtain 


wliere  ^>,  %,  ^  are  symbols  of  functionality. 

It  follows  that 

dr       dr±_ 
+        ~ 


Now  since  by  Art.  11,  R  is  a  function  of  rlt  0l9  and 

dR_dBLdrl     dR  dd±     dR  dz^ 
d6~dr^  ~dd  +W±  ~dd  +  dz   dd  ' 

dRdRdr      dB  dO      dR  dz 


therefore,  by  addition, 

dR     dR     dR 
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whence,  by  the  last  article, 

dR     dR     dR     dR 

W^~fa+lfa+~d5' 

„,  '    4  .    dR   .  ,  dE       ,  dR 

17.     To  obtain  -3—  in  terms  of  -3—  and  -^  . 
de  J    dr  d# 

If  u  denote  the  excentric  anomaly,  we  have* 

r  =  a  (1  —  e  cos  u)  .....................  (1)  , 


nt  +  e—  TX  —u  —  e  sin  u  .....................  (3), 

from  which  r  and  6  may  be  expressed  in  terms  of  t  and  the 
elements  by  eliminating  u.     Assuming  r  arid  6  so  expressed, 

,  ,  .    dr      .  dO 
we  proceed  to  obtain  ^-  and  -r-  . 

-r«         /-,\     dr.          f     .       du  \ 

Irom  (1),    -j-  —  a  le  sin  u  -j  —  cos  uj  , 

and  from  (3),    -7-  (1  —  ecos  w)  —  sin  u  =  0  .............  (4)  ; 


eliminating       -j- ,  we  have 


de 

du 
de 


dr        (    esmu 

-r-  =  a\ — 

de        [1— ecosw 

(  e  —  cos  u  } 
\l—e  cos  u)  ' 


1  —  e  cos  u  —  (1  —  e2) 

6  (1  — ^ 


See  Tait  and  SteeWs  Dynamics,  Arts,  no,  in,  and  122. 
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=  —  a  cos  (6  —  VT), 
by  the  polar  equation  to  the  ellipse. 

Again,  differentiating  the  logarithms  of  equation  (2), 
1          dO      I       1  1    \         1     du 


_  I  f    1  1 

~  2  \1  +  e  +  f 


sin  (0  -  «r)  de  ~  2  \1  +  e      f^e/      sin  u  de  ' 
eliminating  -y-  by  means  of  (4), 


1  1 


sn      —  ^el  —  e       I  —e  cos 


if  h2  —  fjia  (1  —  e2)  ;  therefore 


Now  since  R  is  a  function  of  e  only  because  it  is  a  function 
of  r  and  0, 

de       dr  de      d6  de  7 

.  dR        />      1\    .    ta       .  dB 

r=  —  6t  cos  ( tf  ~~  ffS] 1-  cn  \  —  ~j —  1  sin  ((/  ~~  iff )  — ^  . 

df         \ti       T  /  dfj 

Since  6  —  «r  =  00  —  <GTQ  ,  (see  Art.  13,)  this  equation  may 
be  written 

dR  ^dR 


under  which  form  it  will  be  useful  in  the  next  chapter. 


CHAPTER  II. 

FORMULAE  FOR  CALCULATING  THE   ELEMENTS   OF  THE  ORBIT. 

18.  WE  now  proceed  to  form  equations  of  motion, 
taking  the  Sun's  centre  for  the  origin  of  co-ordinates,  the 
radius  vector  of  the  planet  for  the  axis  of  x,  a  perpendicular 
to  it  in  the  plane  of  the  orbit  for  the  axis  of?/,  and  a  normal 
to  this  plane  for  the  axis  of  z.  With  this  system,  it  will 
be  shewn  that  two  of  the  resulting  equations  can  be  ex- 
pressed in  the  same  forms  as  if  the  planet  moved  in  one 
plane. 

Let  a?,  y,  z  be  the  co-ordinates,  u,  v,  w  the  velocities  of 
the  planet  with  reference  to  three  rectangular  axes  originat- 
ing in  the  Sun's  centre,  and  moving  with  angular  velocities 
0i  ?  02  >  03  about  their  instantaneous  positions :  let  X,,  Y,  Z  be 
the  accelerations  due  to  the  impressed  forces  in  the  directions 
of  the  axes.  Then  (Eouth's  Rigid  Dynamics,  Art.  114), 

dx  ,    -. 


ds 

w=Tt-^ 

and  the  equations  of  motion  are 
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duo 


In  these  equations^,  <£2,  <£3  are  arbitrary;  we  propose 
so  to  determine  them  that  the  axis  of  x  may  coincide  with 
the  radius  vector  of  the  planet,  and  the  plane  of  xy  with  the 
plane  of  the  orbit. 

In  order  that  the  axis  of  x  may  coincide  with  the  radius 
vector  of  the  planet,  we  must  have 

always ;  and  therefore 

dx  _  dr        dy  _         dz 

~di  ~~  dt  '     ~dt  ~    '     di  ~ 

and  in  order  that  the  plane  of  xy  may  coincide  with  the  plane 
of  the  orbit,  we  must  have 

always  ;  and  therefore 

dw  _ 

Hence  equations  (1)  give 
-  — 

and  equations  (2)  become 
,72 

x=w~^' 
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In  order  to  reduce  the  first  two  of  these  equations  to  the 
forms  they  would   take  if  the   motion  were  in  one   plane, 

let  <fe  =  ^|:   thus 

X-^- 
~  dt* 

l^/     ^\ 
~rdt(     dt)9 


19.  If  we  measure  on  the  plane  of  the  orbit  from  the 
planet's  radius  vector  in  a  direction  contrary  to  that  of  motion 
an  angle  equal  to  00,  we  arrive  at  what  may  be  considered 
as  the  origin  from  which  00  is  measured.     Since  this  will  be 
a  point  having  no  angular  velocity  about  the  axis  of  z,  which 
is  normal  to  the  plane  of  the  orbit,  it  is  said  to  be  fixed  in  the 
plane  of  the  orbit*. 

20.  We  shall  for  the  present  confine  our  attention  to  the 
first  two  of  the  above  equations. 

In  order  to  find  the  components  of  the  disturbing  force 
parallel  and  perpendicular  to  the  radius  vector  of  the  planet, 
let  us  take  first  the  radius  vector  as  the  curve  of  reference  ; 
then  s  =  r,  and  R  being  supposed  expressed  as  a  function  of 
r,  0,  fl,  and  i  (see  Art.  13),  we  have 

dE_dEdr_     dR<W      dR^dti      dR^di 
ds       dr  ds      dd  ds      d£l  ds       di  ds 

_dR 

"  dr  > 

*  It  is  necessary  to  define  this  term,  since  the  definition  of  Art.  4  is  not 
sufficient  completely  to  regulate  the  motion  of  the  plane  of  the  orbit,  though 
affording  it  a  distinct  geometrical  position. 
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since  0,  £1  and  i  do  not  vary  with  s.     Hence  the  disturbing 

force  in  direction  of  the  radius  vector  =  -,—  . 

dr 

Again,  let  us  take  as  the  curve  of  reference  a  circle  in  the 
plane  of  the  orbit,  with  its  centre  coinciding  with  that  of  the 
sun  ;  then  $s  =  r£0,  and  we  have 

dR  _ldR 
~3a"~r'd0' 

since  r,  H,  and  i  do  not  vary  with  s.     Hence  the  disturbing 
force  perpendicular  to  the  radius  vector  =  -  -    . 


We  have  then 


Y        jjb      dR       v     1  dR 
~/  +  ^'         =rd8' 


and  the  equations  become 


df  dt     ~       r* 

d  f    d8\      dR 


21.  These  equations  do  not  admit  of  rigorous  integration, 
but  we  may  reduce  them  by  the  method  of  the  Variation  of 
Parameters  to  a  system  of  differentia!  equations  of  the  first 
order.  The  principle  of  this  method  may  be  explained  as 
follows.  Suppose  it  required  to  integrate  the  equations 

dx    dv    d*x    d' 


dx  dy 


where  Pv  P2  are  functions  of  t.     The  solution  of  these  equa- 

tions can  be  made  to  depend  upon  that  of  the  equations 

C.  p.  T.  2 


18  PLANETARY  THEORY. 

<k  =o,  <£2  =  0.     Suppose  the  four  first  integrals  of  these  last 
equations  to  be 

dx 


Cdx     dy\ 
x>  y>  *'  di>  ft) 

t    dx     dy\ 


dx 


where  c1?  c2,  c3,  c4  are  arbitrary  constants  or  parameters. 
The  method  of  the  Variation  of  Parameters  consists  in  so 
determining  cx,  c2  ,  c3  and  c4  as  functions  of  2,  that  these  inte- 

grals (and  therefore  the  two  final  integrals  of  the  equations 
<k  =  o,  <j>2  =  0,  which  can  be  obtained  from  equations  (ii)  by 

ft  rp  //  7/\ 

eliminating  —r-  and  -^J  shall  satisfy  equations   (i).      That 

cr  c2,  c3,  and  c4  can  be  so  determined,  may  be  seen  as  follows  : 
by  the  solution  of  equations  (i),  values  of  x  and  y  and  there- 

fore of  -V  and  -i  can  be  found  as  functions  of  t  and  constant 
dt          dt 

quantities  ;  if  these  be  substituted  in  equations  (ii)  the  requi- 
site values  of  cx,  c2,  c3  and  c4  will  be  obtained.  For  an  ex- 
ample of  the  application  of  this  method,  see  Boole's  Differential 
Equations,  Chap.  IX.  Art.  11. 

22.     If  in  equations  (1)  and  (2)  of  Art.  20  we  put  ^  =  0, 
and  then  integrate  them,  we  obtain 

i  =     {l  +  e  cos  (0,-flr,)}  ..................  (3), 
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where  h,  e,  OTO  are  the  constants  of  integration.  Equation  (3) 
indicates  motion  in  an  ellipse,  of  which  e  is  the  excentricity, 
-sr0  the  longitude  of  perihelion,  and  h  twice  the  area  described 
in  an  unit  of  time.  If  the  mean  distance  in  this  ellipse  be 
denoted  by  a,  we  have  in  addition 

h*  =  tM(l-<?)  ........................  (6). 

We  shall  assume  (in  accordance  with  the  principles  of  the 
method  of  the  Variation  of  Parameters)  the  first  and  second 
integrals  of  equations  (1)  and  (2),  together  with  equation  (6), 
to  retain  the  same  forms  when  R  is  not  zero;  A,  e,  ^Q)  and  a 
being  in  this  case  considered  variable*. 

The  values  of  these  elements  are  to  be  obtained  from  the 
condition  that  the  above  integrals  shall  satisfy  equations  (1) 
and  (2). 

23.  If  their  values  as  calculated  for  any  given  time  be 
substituted  in  equation  (3),  it  will  represent  an  ellipse  having 
a  contact  of  the  first  order  with  the  actual  orbit,  since  the 

values  of  -,-  and  —  7  °  at  the  common  point  will  be  the  same 
at  at 

for  both  curves.  It  is  termed  the  instantaneous  ellipse,  since 
the  planet  may  for  an  infinitely  small  time  be  supposed  to 
move  in  it.  Moreover,  the  velocity  and  direction  of  motion 
of  the  planet  will  be  the  same  as  if  it  moved  in  this  ellipse, 
so  that  if  at  any  time  the  disturbing  force  were  to  cease,  the 
planet  would  continue  to  move  in  the  instantaneous  ellipse 
constructed  for  that  time.  This  is  accordingly  sometimes 
given  as  the  definition  of  the  instantaneous  ellipse. 

*  We  stall  also  for  convenience  suppose  the  equation  %2a3=/i  to  hold  in  the 
disturbed  orbit,  n  being  of  course  considered  variable. 

*-* 
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24.  To  obtain  formulae  for  calculating  the  elements  of  the 
instantaneous  ellipse  at  any  time. 

Suppose  the  value  of  c  required,  where  c  denotes  any  one 
of  the  elements.  From  equations  (3),  (4),  and  (6)  we  may, 
Tby  eliminating  the  other  elements,  obtain  c  as  a  function  of 

r.  00,  TT-  ,  and  h*  :  let  then 


//')'* 

where  r'  is  written  for  -j-  .     Differentiating,  we  have 

dc^dfdr      dfdO,      dfdr^     dfdh 
dt  ;  dr  dt  +  d6Q  dt       dr  dt  +  dh  dt' 

dr      d*r        fd00\*     p      <7R 
Tt  =  df  =  r  (dt)  "  ?  +  Tr 

dh      d       d0\      dR 


dc      dfdr  ,    dfd0Q  ,   df 
therefore  -^  =  -f-  ^7  +  33-  -77  +  x~/          j 

^      <fr  dt      d00  dt      dr   (    \dt  J       r 

dfdR     dfdR 
dr'  dr      dh  dd  " 

But,  since  by  hypothesis,  if  R  were  zero  and  c  constant, 
our  assumed  integrals  would  still  satisfy  the  differential  equa- 
tions, we  have  (making  R  zero  and  c  constant) 


~~drdtdeQ  dt      dr'         dt         r2' 

dc      df  dR     df  dR 
therefore  --j-  =  -f-,  -j-  +  -]  ^  -^  . 

dt      dr  dr      dh  dd 

*  We  retain  h  for  convenience,  in  preference  to  replacing  it  by  r2  —~  , 
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Hence  in  obtaining  the  formulae  for  calculating  the  ele- 
ments of  the  orbit,  we  may  proceed  as  follows.  From  equa- 
tions (3),  (4),  and  (6)  we  may  express  the  element  required 

dr 

as  a  function  of  r,  00,  -y- ,  and  Ji.  We  may  then,  by  differen- 
tiating the  resulting  equation  with  respect  to  t  as  if  r  and  00 

.  .       dR  f      d*r  dR   ,      dli 

were   constants,  writing  -y-  for  -=-5-,  and  -y~-  for  -,- ,  and 

ar  dt  dv  cLt 

eliminating,  if  necessary,   ^-  and  h  by  means  of  equations 

(4)' and  (6),  obtain  the  differential  coefficient  of  the  element 
required  in  terms  of  the  elements,  the  co-ordinates  'of  the 
planet,  and  the  disturbing  force.  The  result,  however,  will 
in  every  case  admit  of  being  expressed  in  terms  of  the  ele- 
ments and  of  the  differential  coefficients  of  R  with  respect  to 
them. 

25.  To  obtain  a  formula  for  calculating  the  mean  dis- 
tance. 

From  equations  (3),  (4)  and  (6),  if  e  and  w0  be  eliminated, 
we  shall  find 


Jt. 

Differentiating  as  if  r  were  constant,  and  writing  -y—  for 

d*r    dR  ,     dh 
,T,  -y-j  for  -y-,  we  have 
dt2 '  dti          dt' 

f*  ^L  -  o  /'^  ^r      ^  ^c 


Now  since  r  —f(nt  +  e  —  *?) , 

C70  —  tc",,  =  6  —  'sr  =  0  (?z£  -f-  e  —  -sr) , 

and  the  forms  of  -y-  and  -y-0  are  the  same  as  if  the  elements 
dt  dt 

were  constant,  we  have 
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and  similarly,  — °  ==  n  -^  =  n  -y- ; 

dt          de         de 

dEdd\ 
dO  de) 


therefore  ^t  =  2    /«*-«'* 


dB 


or        =          dR 
dt        p    de  * 

26.  This  formula  may  also  be  obtained  as  follows.  If  s 
denote  an  arc  of  the  actual  path  of  the  planet  measured  from 
some  fixed  point  to  its  position  at  time  t,  we  have  the  equa- 
tion of  motion 

dfs  dr     dR 


and  by  a  known  formula 

AM*  _  2/A     fju 
\di)  ==T~a* 
Differentiating  the  latter  we  obtain 

r          da 


and,  multiplying  the  former  by  2  -£  , 

cLt 

d<sd^__2/j,dr         dR  ds 
~dtd?~  ~~r*~dt*  2  da  dt* 

,  ,       /,  ^  da     ^  dR  ds 

therefore  ~  -y-  =  2  -y-  -j- 

a  dt         ds  dt 


where      j  '    denotes  the  differential  coefficient  of  R  with 
at 

respect  to  t,  only  so  far  as  it  involves  t  through  involving  the 
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co-ordinates  and  elements  of  the  disturbed  planet.    Now  sup- 
posing E  expressed  as  a  function  of  t  and  the  elements  (Art. 

7  /  -ry\ 

13),  the  form  of  —^  will  be  the  same  as  if  the  elements  were 

invariable  (see  Art.  35)  :  since  then  t  is  always  coupled  with  e 
in  the  expression  nt  +  e,  we  have 

d(R)  _          dR  dR 

dt        U  d(nt  +  e}~  H  de'y 

da      2na?  dR 

therefore  -j-  =  -  •  -j-  . 

at        /ju     de 

27.     To  obtain  a  formula  for  calculating  the  excentricity. 
From  equations  (3)  and  (4),  if  CTO  be  eliminated,  we  obtain 


dt     "  A» 

7TT) 

Differentiating  as  if  r  were   constant,   and   writing  -7- 

" 


72         7.1     '  7  3  ~7Zl  *  2      «7ZJ    ' 

/i     at  n/          a\j      r    ctu 

dh     dR 
since  -T.  =  -TQ  ; 

^ede^^dl^d^     dR  d00     ^(1-e^dR 

(-e2}  dR 


__ 
ndr   de+  dO   de 


.       .              de     ntf  dR     l-e*fdRdR\ 
therefore          -j-  =  —5-  -^ 7 —  (  ^ — 1-  -y- J  ,     (Art.  15.) 


-  e2    dR 


na        -  e 


yu-e          c?e  /ite 

since  ^2  =  ^a  (1  —  e2) ,  and  nV  =  /-t. 
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28.     This  formula  may  also  be  deduced  from  that  of  the 
mean  distance,  by  means  of  the  equation 


™r    i          IT  /-»       2\  aa 

We  have  2A       =/,(!-*)       - 


or 

d 


therefore  __ 

LlltJlCIUiC  i     —  9  ~"V~  "1       7  • 

tw  /*e         ae  pe         \de      far) 

29.     Jb  obtain  a  formula  for  calculating  the  longitude  of 
perihelion. 

From  equations  (3)  and  (4),  if  e  be  eliminated,  we  obtain 
dr        in          .      Ji     u, 

*cot^-^=~i- 

Differentiating  as  if  r  and  6Q  were  constant,  and  writing 

dRf     d*r  ,,  . 

-r-  for  -77  ,  we  obtain 
dr         df 

dr  ,  .  d^  -  dR      fl      u\  dh 


,«  cr  //i  \ 

therefore       -     cosec  (^0  —  -CTO) 


=  —  COS  (VQ  — 

1  dR 


;     (Art.  17.) 


but  from  equation  (4)  -^  cosec  (0Q  —  «TO)  = 

,       ,,  d'sr,.       ri    dR     na  v  (1  — 

therefore          -^  = j-  = * -,— 

at      pea  de  j*e         de 
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Now  if  -CT  denote  the  longitude  of  perihelion  measured  on 
the  plane  of  reference  as  far  as  the  node,  and  thence  on  the 
plane  of  the  orbit,  £1  the  longitude  of  the  node  on  the  plane 
of  reference,  X10  its  longitude  on  that  of  the  orbit,  we  have 


,,       f  n  n 

therefore  _  =  _"  +  ___o. 

Now  -j-  is  the  angular  velocity  of  the  line  of  nodes  on 

/7O 

the  plane  of  reference,  --  r-°  its  angular  velocity  on  the  plane 
of  the  orbit  ; 

«  c?n0    da 

therefore  —  -—-  =  -j-  cos  i  : 

dt        dt 

,1        r  dvr      dor.       .  ..   d£l 

therefore  -y-  =  —^  +  (1  —  cos  i)  -7-  : 

dt       dt       ^  J   dt 

or,  substituting  for  --  from  Art.  31  or  35  y 


dt 


To  obtain  formulce  for  calculating  the  longitude  of  the  nodet 
and  the  inclination. 

30.     We  now  return  to  our  third  equation  of  motion, 


or,  as  it  may  be  written, 

%=^ (!)• 

We  have  seen  (Art.  18),  that  <£2  =  0;  hence  the  motion  of 
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the  plane  of  xy,  which  coincides  with  the  plane  of  the  orbit, 
is  compounded  of  the  angular  velocities  fa  about  the  axis 
of  a?,  and  fa  about  the  axis  of  z.  Now  the  former  is  equiva- 
lent to  an  angular  velocity  fa  cos  (6  —  £1)  about  the  line  of 
nodes,  and  an  angular  velocity  fa  sin  (6  —  H)  about  an  axis 
perpendicular  to  it  in  the  plane  of  the  orbit :  but  the  angular 

velocities  of  the  plane  of  the  orbit  about  these  axes  are  ~  and 

.d£l 
sin  i  —r-  respectively ;  therefore 


dt 

di 

fh    sin  IH  —  {).}  =  sin  i' 

dt 


fa  cos  (0  -  H)  =  -j  ,     fa  sin  (0  -  H)  =  sin  i 


Hence,  by  equation  (1), 

fJi 
hjt  =  Zrcos(0-ty (2), 

h  sin/  -j-**2r  sin  (0- 11) (3). 

dt 

31.  In  order  to  determine  Z  we  must  suppose  the  curve 
of  reference  perpendicular  to  the  plane  of  the  orbit.  If  we  de- 
note by  s  an  arc  of  this  curve  measured  from  some  fixed  point 

up  to  the  planet,  we  have  by  Art.  8,  Z  =  -r- .  Now  the  po- 
sition of  a  point  on  the  curve  of  reference  may  be  determined 
by  its  polar  co-ordinates  r,  6  on  a  plane  passing  through  it 
and  the  sun,  the  inclination  i  of  this  plane  to  the  plane  of 
reference,  and  the  longitude  H  of  its  node.  Since,  however, 
an  infinite  number  of  planes  can  be  drawn  through  two  given 
points,  we  must  introduce  some  further  condition  to  fix  the 
position  of  that  on  which  r  and  6  are  measured. 

First,  then,  let  it  pass  through  SN  the  line  of  the  nodes : 
let  P  be  the  position  of  the  planet,  and  suppose  the  curve  of 
reference  a  circle  AP  with  its  centre  G  in  SN.  Through  SN 
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draw  a  plane  inclined  at  a  small  angle  Si  to  the  plane  of  the 
orbit,  cutting  the  circle  AP  in  p,  and  let  Pp  =  $s ;  then 


or 
therefore 


&s  =  r  sin  (6  —  fl)  Bi ; 

di  _  1 

ds~~  r  sin  (0  —  O)  * 


Now,  (see  Art.  13),  E  may  be  expressed  in  the  form 


in  which,  if  the  series  for  r  and  0  be  substituted,  E  will  be 
expressed  as  a  function  of  t  and  the  elements.     Hence 

dE_dEdr     dE^  d9     dE  d&     dE  di 
ds       dr   ds      d0   ds     d£l    ds       di  ds ' 


,  dr  dQ  d£l  di  1 

and  -T-  =  0,     -7-  =  0,  —r-  =  0,     -7-  =  — : — -TQ — T\\  > 

ds  ds  ds  ds     r  sm  (0  —  O) 

JT>  i  JT> 

therefore 


r  sn      — 


On  substituting  this  value  for  ^  in  equation   (3),  we 
obtain 

,    .     .d£l     dE 

h  sm  %  -j-  =  -JT  , 
dt       di  9 


28 

or  since 
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=  pa  (1  —  e2),  and  ri2a?  —  /JL, 


da 

dt 


na 


dR 


~~       sm 


32.  Again,  suppose  the  plane  on  which  r  and  6  are  mea- 
sured, instead  of  passing  through  8N9  to  pass  through  a  line 
8G  in  the  plane  of  the  orbit  perpendicular  to  SN;  and  take 
for  the  curve  of  reference  a  circle  AP  with  its  centre  C  in  SO. 
Through  8G  draw  a  plane  inclined  at  a  small  angle  to  the 
plane  of  the  orbit,  cutting  the  sphere  in  the  great  circle  nmp9 
and  the  circle  AP  in  p.  Draw  Nm  perpendicular  to  np. 


Then   Nm  which   measures   the    inclination  of  the   two 
planes  =  —  SO  sin  i,  and 


hence  if  Pp  =  Ss, 

§s  —  —  r  cos  (0  —  H)  811  sin  i ; 


therefore 
Again, 


cfe      r  cos  (0  —  O)  sin  i " 
*  =  Ln  +  np-(LN+NP) 

=  nm  —  nN 
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therefore  -7-  =  (1  —  cos  i)  -y- : 

ds  'as 

dr  di 

also  -y-  =  0,    -y-  =  0. 

ds  ds 

dR     dR  de     dR  dfl 

Hence  -j-  =  -TTT  -y-  4-  ^7^  - /- 

as       de/  ds      ail  as 

fd^  .N  d^)  1 


rcos  (0—  H)sin  2" 
On  substituting  this  value  for  ^in  equation  (2),  we  have 
di  1        dR 


-  cos 


) 

jTT  f 

du) 


•       r  j 

A  sin  *  (ail  x  du 

I       (dR         .  2  *  /^     dR\     .  . 

«-T-S  —  ••bo-Msin1-  hr+^")f  (Art 

A  sm  ^    c/ll  2 


° 


33.  The  only  remaining  element  is  the  epoch,  but  before 
proceeding  to  obtain  a  formula  for  its  calculation,  we  shall 
give  another  method  of  obtaining  the  results  of  the  last  two 
articles. 

34.  To  obtain  a  formula  for  calculating  the  inclination. 
(Second  method.) 

If  the  motion  of  the  planet  be  referred  to  the  polar  co-ordi- 
nates of  its  projection  on  the  fixed  plane  of  reference  and  its 
distance  from  this  plane,  we  have  the  equation 

1  d  (  2  dO\      1  d 

r^^J-r 

d_  /  2  d£\  _  dR 
cft'V!  ~dtJ~dOC 
Now  if  8^4  ,  S^j  denote  the  vectorial  areas  swept  out  in 
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the  time  St  on  the  plane  of  the  orbit  and  the  plane  of  reference 
respectively,  we  have 

8^=  $A  cos  i: 

but  M-  :'***    Li- 


2  2      - 

therefore  r^   ^  =  r  z  -?-  tos  ^  =  h  cos  ^. 

Hence  our  equation  of  motion  becomes 
d  ,       dR 


,    .     .di  .dh     dR 

-Asn^  +  ccm^^- 

_   .    .di     dR  .dh 

therefore  —  Asnu  -    —  -      —  cos  ^  - 


__ 
"" 


na         (  1     dR  i  (dR     d±. .. 

-.-777  + tan-    -T-  +  -T-  r- 


jj,  V(l  —  e2)  [sin  tdQ,  2  \ 

35.  To  obtain  a  formula  for  calculating  the  longitude  of 
the  node.  (Second  method.) 

Since  the  velocity  of  the  planet  at  any  time  can  be  ex- 
pressed in  terms  of  the  co-ordinates  and  elements  of  the  instan- 
taneous ellipse  constructed  for  that  time,  in  the  same  form  as 
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if  it  moved  iii  this  ellipse,  its  component  in  any  direction  can 
also  be  so  expressed.  Hence,  considering  R  as  a  function  of 

r1,  Bv  z  (see  Art.  11),  the  values  of  -~  ,    -^  ,   -7-  ,  and  there- 

fore of     ,     (where  —  ^  denotes  that  R  is  to  be  differentiated 
at    ^  at 

with  respect  to  t  only  so  far  as  it  involves  t  through  involving 
the  co-ordinates  of  the  disturbed  planet)  may  be  expressed  in 
the  same  forms  as  if  the  elements  were  invariable. 

Now  we  have  seen  (Art.  13)  that  R  may  be  thus  ex- 
pressed : 

#-/(*,*,  o>t)i 

or  since  evidently,         6  —  H  =  0Q  —  I10, 


We  have  then,  considering  the  elements  variable, 
d  (R)  _  dR  dr     dRd(6Q+Q,-£l0)     dRdti     dR  di 


dt.    ~d       dd0~     ~3T 
and,  considering  them  invariable, 


dRdO 


^ 

~  + 


^ 
dt     ~  dr  dtd0Q  dt  "    * 

Equating  the  two  values  of  ,  we  obtain 


„ 

ow 


_ 

da~di+~dJdt~ 
dR     dR 


*  We  have  made  this  transformation,  because,  although  the  value  of  -^  is 

the  same  in  form  as  if  the  elements  were  invariable,  this  is  not  the  case  with 

dd 

dt' 
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,  ,        .   4    _N    d&0     d£l 
and  (see  Art.  29),   -^-°  =  -=-  cos  i  ; 

(dR     „  -  fdB     dB\\  d£l  ,  dR  di 

therefore  \-j^  +  (1  -  cos  i)   -j-  +-=-  }\  -j-  +-rr  ^=0. 
\d£l  '  V  de      d^i]    dt       di  dt 

Substituting  for  -r-  its  value 

d£l  na  dR 

dt  ~~  fj,  V(l  —  e2)  sin  i  di  " 

36.  Jb  obtain  a  formula  for  calculating  the  longitude  of 
the  epoch. 

If  R  be  expressed  in  terms  of  t  and  the  elements  (see 
Art.  13),  since  nt  +  e  always  occurs  as  one  symbol,  we  may 

write 

R—f(nt  +  e,  a,  e}  TV,  O,  i). 

Differentiating,  the  elements  being  considered  variable^ 
we  have 

d(R)  _dRd(nt+e)      dRda     dRde     dRd^ 
dt         de         dt  da  dt       de  dt      d^  dt 

djRdQ     dRdi 
+  ~d&  dt  +~di  dt' 

and  differentiating  as  if  the  elements  were  invariable,  which 
is  permissible  for  the  reason  explained  in  the  last  Article, 

d(B)  _     dR 
dt     ~n  de' 

Equating  the  two  values  of 


dt 

dRL_dRLf          dn     de\      dRda     dR  de     dR 
~n'    t+ 


dRda     dRdi 
d&  dt  +  di  dt' 
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Substituting  for  -y-  ,   -7-  ,  &c.,  their  values 

(it        Out 


—     —  _L     na      ~ 


_ 
de       dt      dt          p     de  da  pe        de   de 

-e2)  fdR     dR\dR     na  V  '(I  -  e*)  dE  dR 
\de      d^ff)  de  fjue          de  dix 


nqtan2    dRdR  na  dR  dR 


-e<<J    di  d^jul-ez&mi  di 


na  dRdR  2    fdR     dR\dR 

-e2smida  di'  * 


dn     de\      2na  dR  dR 
de  V~di  +  dt)'{'~ir~de~~da 

-  na  V(1 "  *2)  {1  -  V(l  -  e2)}  ~  ~ 
l*&  J  de  de 

• 

m2    ^^ 


-e')  de  di' 

rJTf 

Dividing  every  term  by  -y- ,  and  transposing,  we  obtain 
de  dn     Ina*  dR 


dR 


37.     Of  the  formulas  which  have  been  obtained  for  cal- 
culating the  elements  of  the  orbit,  that  of  the  preceding  article 
is  the  only  one  which  contains  a  term  proportional  to  the 
C.  P.  T,  3 


34  PLANETARY  THEORY. 

time*.  It  may,  however,  be  replaced  by  one  in  which  no 
such  term  exists.  For,  let  f  denote  the  mean  longitude, 
then 


d£  dn     de 

therefore 


or 


Now  let  f  =  fndt  +  e', 

then  by  the  formula  of  the  last  article 


de          2na*  dR     na*J(l-e*}  .          ,,         ~,  dR 
=  *      ~~  "V(l<    e]] 


dR 


Since  in  the  elliptic  formulae  e  never  occurs  except  when 
coupled  with  nt,  in  the  expression  nt  +  e,  it  will  be  altogether 
eliminated  if  for  nt  +  e  we  write  Jndt  +  e'. 

4 

Considered  as  replacing  the  element  e,  e'  is  called  the 
epoch  *.  Since  however  we  shall  never  have  occasion  to  em- 
ploy the  formula  of  the  preceding  article,  the  accent  will  in 
future  be  omitted. 

fndt  is  termed  the  mean  motion  in  the  disturbed  orbit, 
and  is  denoted  by  £ 

38.     To  obtain  a  formula  for  calculating  f,  we  have 
<T£_dn 
d?  ~~  dt  ' 

*  The  reader,  if  acquainted  with  the  Lunar  Theory,  will  have  already  seen 
the  inconvenience  of  such  terms. 

•f1  It  may  be  noticed  that  if.w  were  constant,  e'  would  be  identical  with  e. 
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and  differentiating  the  equation  wV  =  p 

»dn      ,  „  9  da 
2wa3  -j-  +  3nV  -j-  =  0  ; 


therefore  -77  =  —  —  -77  =  ---  -y-  > 

c?^          2a  at  fj,     de  ' 

d*%         Ztfa  dR 

or  -r    =  --     -  --  . 


39.     We  will  here  recapitulate  the  formulae  which  have 
been  obtained  for  calculating  the  elements  of  the  orbit. 

da  _  2na?  dR 
dt**    ft   ~3e' 

de  _na(l-  e2)  dR     na  V(l  -  e2)  fdR     dR\ 
^  '      dt~        fie         de  ^          \de      dvrj' 


~dt~  ie          de+  i^/l-e2    di' 


de         Zna^dR     ^V(l-e2)  ,,         -.dR 


dR  i(dRi 

+  Ln+ 


We  have  also  (vii)  the  equation 

d?%__     ZtfadR 
df  '      ~T  de' 
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but  this  forms  no  new  relation,  since  it  has  been  deduced 
from  (i). 

40.  When  the  elements  have  been  calculated  by  means 
of  the  above  formulae,  the  position  of  the  planet  will  be  given 
by  the  equations 


0  =  f  +  e  +  26  sin  (f  +  e  -  «r)  +    e2  cos  2  (f  +  e  -  «r)  + 


CHAPTER   III. 

DEVELOPMENT   OF  THE   DISTURBING  FUNCTION; 

41.  IN  the  first  chapter  we  have  obtained  equations  by 
means  of  which  R  may  be  expressed  in  terms  of  the  time 
and  the  elements  of  the  orbit  ;  we  now  proceed  to  shew  how 
the  actual  development  may  be  effected  in  a  series  ascending 
by  powers  and  products  of  the  excentricities  and  the  tangents 
of  the  inclinations.    In  the  Planetary  Theory  these  are  ex- 
tremely small,  and  the   series  will  converge  rapidly.     Ac- 
cordingly in  the  present  treatise  small  quantities  of  orders 
higher  than  the  second  will  be  neglected*. 

42.  If  we  recur  to  Art.  11,  it  will  be  seen  that,  consider- 
ing only  one  disturbing  planet, 


== 

1 


2  +  r*  -  2r^  cos  (0,  -  0/)  +  (z  - 


The  first  step  towards  the  required  development  will  be 
the  expansion  of  rv  r^,  6V  0/,  z  and  z  in  terms  of  the  time 

*  We  may  remark  that  to  this  order  of  approximation  the  inclinations,  their 
sines,  and  tangents  will  be  equal. 
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and  the  elements  of  the  orbit.  For  this  purpose  we  may 
employ  the  equations  which  have  already  been  obtained  in 
Art.  13,  viz.  : 

tan  (0X  -  O)  =  cos  i  tan  (9  -  O)  , 
sin  X  ==  sin  i  sin  (0  —  Q), 
rx  =  r  cos  X,         s  =  r  sin  X, 

—  OT)  —  -e2cos2  (n  +  e  —  OT)  —  ...k 

0  =  nt  +  e  +  2e  sin  (nt  +  e  —  *r)  +  -  e2  sin  2  (ntf  +  e  —  «•)  +  ..., 

with  similar  equations  involving  the  co-ordinates  and  elements 
of  the  disturbing  planet. 

(i)     To  expand  rv.    We  have 

rt  =  r  cos  X  =  r  (1  —  sin2  X)^ 


=  r  ji  _  1  sin't  sin*  (0  -  H)  +  ...i 
=  r  |l  -  itan»  t  sin*  (0  -  fl)  +  ...I 


to  the  same  order  of  approximation, 


or  substituting  the  expansions  for  r  and  0, 

rx  =  a  •!  1+  -  e2  —  -tan2^  —  e  cos  (nt+  e  —  «r)  —  -e2  cos  2  (w^+e 


suppose. 
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Similarly,  r^  =  a  (1  +  u). 

(ii)     To  expand  0t.     We  have 
tan  (0,  -  0)  =  tan  {(0t  -  U)  -  (0  -  U)} 

tan  (^-H)- tan  (0-11) 
"  1  +  tan  (6l  -  U)  tan  (0  -  11) 

_  (cos^-1)  tan  (0- 11) 
"  l  +  cositan2(0-H) 

-2  sin2 1  tan  (0-11) 

=    •••M-—-— -.-.—»•••=— »«—™->——.-«™—— — ——-———. 

1  +  tan2  (0  -  to)  -  2  sin2|  tan2  (0  -  H) 


=  -  sin2|  sin  2(0-11)  -  ...  ; 
therefore        0X  -  0  =  -  sin2 1  sin  2  (0  -  U)  -  ... 

=  -  tan2|  sin  2  (0  -  U)  -  ..., 

to  the  same  order  of  approximation;   or,    substituting  the 
expansion  for  0, 

5 

0i  =  nt  +  e  +  2e  sin  (w£  +  e  —  •&)  +-r&  sin  2  (nt  +  e  —  ty) 

4 

-  tan2|  sin  2  (nt  +  e  -  U)  +  ••• 
=  nt  +  €  +  v,  suppose. 
Similarly,  0/  =  n't  +  e  +  v'. 
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(iii)     To  expand  z.    We  have 

z  —  r  sin  X  =  r  sin  i  sin  (6  —  O) 
=  r  tan  i  sin  (^  —  O)  —  ... 

to  the  second  order;  or,  substituting  the  expansions  for  r 

and0, 

z  =  a  (tan  /sin  (TZ£  +  e  —  O)  -f  ...]. 

A  similar  expression  may  be  found  for  z. 

43.  Having  obtained  the  expansions  of  rt,  r',  6t,  0', 
z,  z'  we  must  now  substitute  them  in  the  expression  for  E. 
This  may  be  effected  as  follows. 

Let  R'  be  the  value  of  R  when  u,  u,  v,  v  are  severally 
zero  :  then,  writing  <j>  for  nt  4-  e  —  (nt  +  e')  we  have 

R'  =  m'  [{a2  +  a'2  -  2aa'  cos  <£  +  (z  -  z'Y}^ 


=  m'  [(a2  +  a'2-  2aa  cos  <£)"*  -  -^-2  cos 


m'       (a2  +  a'2  -  2aa'  cos  ^>)-   («  -  z')z 

,  -1      »     3  a    „ 

" 


AT  T>  Tlf  /      r  / 

IN  ow  R  =  R  +  -j-  au  +  -r-7  au  +  -=7-  (v  —  t?  ) 

v 


aa 


,     ,  ,  , 

+  -7—7-7  aa  uu  +.    -,,au  (v—v)+  ,  ,  ,  ,  aw  (v  -  v  ) 

s         '      da  d 
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44.     It  will  be  shewn  in  a  subsequent  article  that 

(ft2  +  a/2  -2ftft'  cos  <£)-', 
can  be  expanded  in  a  series  of  the  form 


Assume  then 

"*  = 


(a2  +  a'2  -2a«'  cos  <£)"*  =  -  <70+  (^cos^-f-  <72cos2<£  +  ... 


Thus 


m  au 


ldCQ      fdC.       l\  dC2  } 

-  -j-*-  +  [-J-  1  -  -J  cos  <f>  +  -j-2-  cos  26  +  ....  X 
2  rfa       \da      a  J  da  J 
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-  m'au  (v  -  »')  j(^-*  - -^  sin  <j>  +  2  -j±  sin 

-  »W  (,  - O  jgj  +  g  sin  *  +  2  g*  sin  2*  +  ...} 


,  (zz       3  a 


45.     By  Art.  42, 

1.1  1 

w  =  -  e  —  j  tan  i  —  e  cos  (nt  +  e  —  w)  —  e2  cos  2  (w£  +  e  —  -sr) 

^5  4  2 

1 
+*rtan«  cos  2  (?i£-{-e  —  fl)  +  ... 

4 

v  =  2e  sin  (w£  +  e  —  «r)  -f  -  e2  sin  2  (n«  +  e  -  w) 
4 

-  tan2 1  sin  2  (n*  +  e  -  fl)  +  . . . , 
2  =  «  (tan  i  sin  (n£  +  e  —  H)+  . . .}, 
with  similar  expressions  for  u ,  v',  2'. 
Hence  w2  =  e2  cos2  (w^  +  e  —  «r)  +  . . 


(v  -  v')2  =  4e2  sin2  (n<  +  e  -  «r)  +  4e'2  sin2  (w'«  +  e'  -  w1) 

—  See'  sin  (nt  +  e  —  -or)  sin  (n'<  +  e'  —  «r')  +  . . . 
=  2  (e2  +  e'2)  -  2e2  cos  2  (n*  +  e  -  w)  -  2e'2  cos  2  (w'*+er  -  «r') 
-  4ee'  cos  ((/>  -  «r  +  «•')  +  4ee'  cos  {(/i  +  n')  «-fe  +  e'—or— sr'}  +. .., 
ww'  =  eef  cos  (n^  +  e  —  -sr)  COS(TI^  +  e  —  «•')  +  ... 
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2 
u  (v  _  #')  =  —  e2  sin  2  (nt  +  e  —  vr) 

+  2ee  cos  (nt  -f  e  —  OT)  sin  (rc'£  +  e'  —  -BJ')  +  . .. 
=  —  e2  sin  2  (w$  +  e  —  «r)  —  ee'  sin  (^  —  w  +  «r') 

+  ee  sin  {(w  +  w')  t  +  e  +  e'  —  tar  —  *r'}  +  . . ., 

>.  9     •  FOj  9*^  2 

'  tan  ^     a   tan  ^      a 


A1*  tan**1 

cos2(^+6f-Qy)-^tan^tan^cos(^~Q+^/) 


&c. 

46.  If  these  values  be  substituted  in  Art.  44,  it  will  be 
seen  that  cosines  will  be  multiplied  only  by  cosines,  and  sines 
by  sines.    Hence  the  series  will  consist  of  two  parts,  one  inde- 
pendent of  t  explicitly,  and  the  other  consisting  of  periodical 
terms  of  the  form 

Pcos{(pn±qn')  t+  Q}, 

where  p  and  q  are  any  positive  integers  or  zero,  P  is  a  function 
of  the  mean  distances,  excentricities,  and  inclinations  ;  and  Q 
a  function  of  the  longitudes  of  perihelia,  nodes,  and  epochs. 
The  former  part  is  denoted  by  the  symbol  F:  we  proceed  to 
determine  its  value  as  far  as  the  second  order  of  small  quan- 
tities. 

47.  To  determine  that  part  of  R  which  is  independent  of 
the  time  explicitly. 

If  those  terms  only  be  written  down  which  either  are,  or 
after  reduction  will  become,  independent  of  t,  we  have 

,(C0     adC0  fe*     tan2 1\      a'  dC0  (e2     tan2 

"~       ~~~"~~ 
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azd*C0e*  ,  ^^o^_A/altanV     o^tanV\ 
4    da2    2  +  4   da'2   24^2  2       / 


+  -  (  Cj  --  -2   cos  (j)  4ee'  cos  (^  —  -sr  +  «•') 


f  -7^  --  ^  J  sin  <^>  ee  sin  (0  —  w  +  «•') 
f 


-5    sn   >  ee  sn       — 
a 


+  —  -  cos  ^  ad  tan  z  tan  ^'  cos  (<j>  —  11  +  fl')  +  ...!-. 


Now  cos  <^>  cos  ((j)  —  &  +  iff')  and  sin  <f>  sin  (<£  —  vr  —  -BT')  con- 
tain the  term  J  cos  (w  —  «•'),  cos  <£  cos  (^  —  O  +  fl')  contains 
the  term  J  cos  (fl  —  II')  ;  hence 


i 

2       4 


4-  j  aa'  Dx  tan  i  tan  i'  cos  (H  —  ft')  +  .  .  .  r  . 

We  shall  hereafter  be  able  to  simplify  this  expression. 

48.     We  have  seen  that  the  remaining  terms  of  E  are 
of  the  form  Pcos  {(pn±qn)  t+  Q}:  if  then  values  of  p  and 
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q  could  be  found  such  that  pn  ±  qn  —  0,  this  term,  being  in- 
dependent of  t  explicitly,  would  form  an  additional  term  in  F. 
No  instance  of  this,  however,  occurs  among  the  planets. 

49.  In  consequence  of  the  extreme  smallness  of  the  ex- 
centricities  and   inclinations  of  the  planetary  orbits,   terms 
in  R  of  orders  higher  than  the  second  may  in  general  be 
neglected  :  but  it  sometimes  happens,  as  in  the  Lunar  Theory, 
that  higher  terms  become  sensible  through  the  process  of 
integration.     This  we  shall  consider  in  a  subsequent  chapter, 
but  the  following  proposition  has  an  important  bearing  on 
the  subject. 

50.  The  principal  part  of  the  coefficient  of  a  term  in  II 
of  the  form  P  cos  {(pn  —  qn)  t  +  Q]  is  of  the  order  p  ~  q. 

DEF.  By  the  principal  part  of  the  coefficient  is  meant 
that  part  of  P  which  is  of  lowest  dimensions  in  e,  e',  tan  i, 
tan  i'. 

If  we  return  to  the  expression  for  R  in  Art.  44,  it  will 
be  seen  that  in  order  to  obtain  the  general  term  it  will  be 
necessary  to  multiply  the  product  of  the  general  terms  of  the 
expansions  for  ua,  u'P,  v"*,  v's,  2e,  z'£  by  cos  k  (f>  or  sin  k  <f>. 

Now  (1)  in  the  expansions  of  u,  u,  v,  v ',  z,  z'  the  follow- 
ing law  is  observed  to  hold : — The  number  which  multiplies 
nt  +  e  or  n't  +  e'  in  the  argument  of  any  term  represents  the 
order  of  the  principal  part  of  the  coefficient  of  that  term. 

(2)  The  same  holds  good  in  any  power  of  w,  u',  v,  v', 
z,  or  z'.  For,  consider  a  term  P  cos  (put  +  q)  in  uz.  It  can 
only  have  arisen  in  the  following  ways;  partly  from  the 
multiplication  of  two  terms  in  u  of  which  the  arguments  are 
Int  +  X  and  mnt  +  y^,  where  I  +  m  —p ;  and  partly  from  such 
as  have  the  arguments  Int  +  V  and  mnt+ftj  where  I  ~m=p. 
In  the  former  case  the  order  of  the  coefficient  will  be  l  +  m, 
which  equals  p,  in  the  latter  it  will  be  I'  +  m,  and  this  is 
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greater  than  p.     Hence  the  principal  part  of  the  coefficient 
of  a  term  P  cos  (pnt  +  g)  in  w2,  will  be  of  the  order  p. 

Since  then  the  law  holds  in  w2,  it  may  be  shewn  in  like 
manner  to  hold  in  the  product  of  uz  and  u,  i.e.  in  us.  Thus 
it  may  be  proved  for  any  power  of  u.  In  like  manner  it 
may  be  shewn  to  hold  for  any  powers  of  u,  v,  v',  z,  or  z. 

(3)  The  same  law  is  true  for  the  product  of  any  powers 
of  u,  v,  z ;  and  likewise  for  the  product  of  any  powers  of 
u,  v,  z .     This  may  be  proved  by  a  method  similar  to  that 
of  (2). 

(4)  In  the  product  of  any  powers  of  u,  u,  v,  v',  z  and 
z',  the  order  of  the  principal  part  of  the  coefficient  is  the 
arithmetical  sum  of  the  multipliers  of  nt  and  n't. 

For  let  us  consider  a  term  If  cos  {(In  ±  I'ri)  t  +  N}.  Now 
this  must  evidently  have  arisen  from  the  multiplication  of 
L  cos  (Int  +  X)  with  L'  cos  (I1  n't  +  X'),  or  of  L  sin  (Int  +  X) 
with  L'  sin  (l'rit  +  \r),  where  by  (3)  L  is  of  the  order  I  and 
L'  of  the  order  I'.  Hence  H  will  be  of  the  order  I  +  I'. 

Now  atiy  term  in  the  development  of  R  of  the  form 
Pcos  {(pn  —  qn]  t+  Q]  must  have  arisen  partly  from  the 

cos 
multiplication  of  Pt    .    &(/>,  or  as  it  may  be  written 

S1U 


with  P2  «  '  {[(p  -Jc}n-(^-Tc}n'}t  +  <y, 

Dili 

and  partly  from  its  multiplication  with 


sin 


where  Tc,  is  any  positive  integer  or  zero,  PA  is  a  function  of  a 
and  a  only,  and  P2,  P8  are  functions  of  the  excentricities  and 
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inclinations,  such  that  the  orders  of  their  principal  parts  are 
given  by  law  (4).  Hence  the  order  of  the  principal  part  of 
P  will  be  equal  to  the  lesser  of  those  of  P2  and  P3. 

Now  the  order  of  the  principal  part  of  P2  will  be  the 
least  value  of  which  the  arithmetical  sum  ofp~Jc  and  q  ~  Jc 
is  susceptible,  for  different  values  of  k. 

(i)     Suppose  k  intermediate  to  p  and  q  ;  then  this  sum 
=p~k  +  q<-k=p~q', 

(ii)  Suppose  Jc  not  greater  than  the  smaller  of  p  and  q  ; 
then  this  sum  —p  +  q  —  2k,  the  least  value  of  which  (by 
putting  Jc  equal  to  the  smaller  of  p  and  q)  =p  ~  q  ; 

(iii)  Suppose  Jc  not  less  than  the  greater  of  p  and  q; 
then  this  sum  =  2Jc—  p  —  q,  the  least  value  of  which  (by 
putting  k  equal  to  the  greater  of  p  and  q)  =p~q. 

Thus  p  ~  q  is  the  order  of  the  principal  part  of  P2.  That 
of  P3  will  be  the  least  value  of  which  p±k  +  q  +  k  is  sus- 
ceptible, i.e.  p  +  q. 

Hence  it  appears  that  the  order  of  the  principal  part  of 
P  is  p  ~  q. 

51.  The  principal  part  of  the  coefficient  of  a  term  in  E, 
of  the  form  P  cos  {(pn  +  qn')  t  +  QJ  is  of  the  order  p  +  q. 

This  term  arises  from  the  multiplication  of  such  terms  as 


P2  .n  [Up-Qn  +  te  +  K)  n']  t+Q,}, 


and  P 


and  as  in  the  last  Article,  the  order  of  the  principal  part 
of  P  will  be  equal  to  the  lesser  of  those  of  P2  and  P3. 
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Now  the  order  of  the  principal  part  of  P2  will  be  the 
least  value  which  the  arithmetical  sum  of  p  ~  Jc  and  q  +  k 
can  assume,  for  different  values  of  k. 

(i)     Suppose  Jc  less  than  p  ;  then  this  sum 

=p  —  Jc  +  y  +  Jc  =p  +  q. 
(ii)     Suppose  Jc  not  less  thanp;  then  this  sum 

=  k  —p  4-  g  +  Jc, 
the  least  value  of  which  (by  putting  Jc  equal  to  p)  =p  +q. 

Similarly  it  may  be  shewn  that  p  +  %  will  be  the  order 
of  the  principal  part  of  P3. 

Hence  it  follows  that  p  +  g:  will  be  the  order  of  the 
principal  part  of  P. 

In  Art.  44  we  have  assumed  that  (a2  4-  a'2  -  2aa  cos  <£)"* 
can  be  expanded  in  a  series  of  cosines  of  </>  and  its  multiples, 
we  shall  now  give  a  proof  of  this  and  shew  how  the  coeffi- 
cients may  be  calculated. 

52.  To  shew  that  (a2  +  a'2  —  2aa'  cos  $)~*  can  be  expanded 
in  a  series  of  cosines  of  multiples  of<$>. 

Suppose  a  greater  than  a',  and  for  —  write  a  ;  then 

Ob 

(a«  +  a'2  _  2aa  cos  ft'8  =  <T28  (1+  a2  -  2a  cos  <£)-* 


1  +  a2  -  cc  e       + 


'fr+yfr+2)ay*-+... 

U? 
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,   30Vri          } 

"'I 


where   the   coefficients  of  e**       and  e-'1  will   always  be 
equal.     Hence  we  may  write 


=  -  A0  +  A1co$<l)+A2Qos2<l>+  ...  +^4fc  cos  &<£  +  ...  , 

where  AQ,  A±,  &c.,  are  functions  of  a  and  a.  The  series 
which  they  represent  will  be  always  convergent  provided  a  is 
less  than  unity,  or  a  greater  than  a.  If  a  be  less  than  a',  we 
have  only  to  interchange  a  and  a  in  the  above,  so  that  a  will 
then  denote  the  ratio  of  a  to  a. 

53.     To  calculate  C0  and  Ct. 

In  the  preceding  article,  let  s  =  -  ;  then 


C.P.T. 
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Unless  a  be  small,  these  series  will  converge  too  slowly 
to  be  practically  useful.  More  convergent  series  might  be 
obtained,  but  according  to  Pontecoulant  (Systeme  du  Monde, 
Tome  III.  p.  81),  it  is  more  convenient  to  employ  elliptic 
integrals  for  the  purpose,  in  the  manner  we  proceed  to  ex- 
plain. We  have 


,     ft  »«  *  §  \  ±  ,-*.     ^-^n  v-wt*    f     *      £-k          1    i 

(a  +  a   —  2aa  cos  0)  '2      2  2 

+  -  (72  (cos  30  +  cos  0)  +  . 

Integrating  both  sides  of  these  equations  with  respect  to 
between  the  limits  0  and  2-Tr,  we  obtain 

r  2-rr  7  -t  "if*  %1T  7  -L 

„,       I                      d(b                       1  I                   am 
7rC0=l     — z— ; 7T4  =  aJ    7 a —     ^1 ' 


PO^ 

TrO, 


_  r2"-        cos  ^  d$ 

1~J0  (a*  +  d*-2aa 


o    ^w.   1-  tt    —  zio/t*   cue  «^>y  ~        —  "  o    (1  +  &.   —  20C  COS  0J  'J 

These  integrals  may  be  reduced  to  the  standard  forms  of 
elliptic  functions  by  assuming 

sin  (0-0)  =  a  sin  6 (1), 

whence  tan#  =  —      ^     (2). 

cos  (p  —  a 


From  (1)    cos  (0  -  (jf 

,       ,,  c?0     cos  (9  —  0)  —  oc  cos 

therefore  -£  = — — 7:5 — -, 

du  cos  (c/  —  <p) 


DEVELOPMENT  OF  THE   DISTURBING  FUNCTION.  51 

Now 

cos  (0  _  <£)  _  a  Cos  6  =  cos  6  (cos  <j>  —  a)  +  sin  0  sin  <£ 

/cos20        •    /A    •     •    t.     /<w 
**     7  (  )? 


in2<#,by  (2), 
=  V(l  +  a2-  2a  cos  <£)...  .........  (3). 

Also  cos(0-<£)=V(l-a2sin260  ...................  (4). 

d6 


Again,  from  equations  (3)  and  (4) 

V(l  +  a2  -  2a  cos  ^>)  =  V(l  -  a2  sin2  6}  -  a  cos  Q  ; 
therefore   1  +  a2  —  2a  cos  (£  =  1  -  a2  sin2  0 

+  a2  cos2  0  -  2a  cos  19  V(l  -  a2  sin2  0), 
2a  cos  <£  =  2a2sin2  0  +  2a  cos  6  V(l  -  a2  sin2  0), 
or  cos  <£  =  a  sin2  0  +  cos  0  V(!  —  «2  sin2  ^)» 

Now  as  <£  increases  from  0  up  to  2  TT,  0  also  increases  from 
0  to  2-7T  ;  hence 


—  2acos 


=-  f 

»•*'„ 


4-2 
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1    [2«     asin2<9dfl  1 

=  --  -77^  -  a     •    */i\  H  --- 

a^Jo   V(l-«  sm'0)      a7r 

i    ff27r        do  c 

=  —  *7  . 

aa?r  V0   V(l  -a"  sm20)     J 


COS 


Hence  with  the  usual  notation  for  elliptic  integrals  (see 
Todhunter's  Integral  Calculus,  Art.  222), 


The  numerical  values  of  F[a9^}  and  E[a,  ^)  may  "be 

\     «/  \     *J 

found  from  Legendre's  tables  of  elliptic  functions. 

54.     Given  Ck  and  C^  to  obtain  Ck+1. 
We  have  .  • 

(a2  +  a'2  -  2aa  cos  0)-i  =  1  (70  +  C7X  cos  ^  +  .  .  .  +  (7,  cos  ^  +  .  .  .  ; 

differentiating  with  respect  to  </>, 

aa  sin  </>  (a2  +  a'2  -  2aaf  cos  0)"1  =  Q  sin  ^>  +  2  Ca  sin  2(^>  +  .  .  . 


therefore  aa1  sin  <f>  (-  CQ  +  <7X  cos  </>  +  .  .  .  J 

=  (a2  +  a'2  -  2aa'cos  ^)  (^  sin  ^+  2C2  sin 
equating  coefficients  of  sin  7c<f), 
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2k 
whence         C    =  -r-  Ck  - 


55.     6rivew  Ck  and  Ck+1  to  obtain  Dk. 
As  in  the  last  article,  we  have 

ad  sin  </>  (a2  4-  a"2  -  2aa'  cos  <£)~f  =•  (7,  sin  $  +  2  <72  sin  2<£  +  .  .  . 

+  fc  Ok  sin  1c$>  +  .  . 

therefore  aa'  sin  ^  (-DQ  +  D^  cos  ^  +  Dz  cos  2^>  +  .  .  .1 

=  0,8111^  +  2^  sin  2£  +  ...  ; 
equating  coefficients  of  sin  /£(/>, 

tkC>  =  <ui(D^-D^  ..................  (1), 

writing  Jc  +  l  for  -A:, 

2(A  +  l)Citl  =  aa'(A-l>*«)  ...............  (2). 

Again, 

(a2  +  a'2  -  2aa  cos  *)"*  =  5  A  +  A  cos  <f>  +  D2  cos 


and 

(a9  +  a'2  -  2aa  cos  ^)~^  =  |  (70  +  (7X  cos  ^>  +  (72  cos 

therefore  -  <7  +  Ci  cos 

4S 


-  2aa'cos  ^)         0+  A  c 

equating  coefficients  of  cos  &0, 

0,=  (a?  +  a'*)  Dt-aa'  (D^  +  flJ  ..........  (3), 

writing  Tc  +  I  for  &, 

GM  =  (a"  +  a'2)  D^  -  aa'  (Dt  +  D^}  .........  (4). 

Eliminating  Dk_:  between  (1)  and  (3), 

1)  Ck  =  (a2  +  a'2)  DH  -  2aa'Dtfl  .........  (5). 
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Eliminating  JD^  between  (2)  and  (4), 


(6). 

Finally,  eliminating  Dk+l  between  (5)  and  (6), 
(a2+  a'*)  Ck- 


or  Dt  =  ,     {(a2  +  a")  ft  - 


56.     To  calculate  the  successive  differential  coefficients  of 
Ck  and  Dk  with  respect  to  a  aw^  a'. 

We  have 
(a2  +  a'2  -  2aa'  cos  0)~*  =  -C0+  Cl  cos  ^  +  (72  cos  2^>  +  ... 

2 

4-  Ci  cos  ^  +  .  .  .  : 
differentiating  with  respect  to  a, 

-(a  -a'  cos  <£)  (a2  +  a'2  -  2aa'  cos  $"*  =         *  +        cos  *  +— 


substituting  for  (a2  +  a**  —  2aa  cos  ^>}~^  its  expression  in  series, 
-  (a-  a'cos<£)  (-Dt  +  Dt  cos  ^+  ...  -f  Dk  cos  ^0+  ...) 


2  c?a        ^a 
equating  coefficients  of  cos  k<f>,' 


By  giving  to  A;  in  succession  the  values  1,  2,  3,  &c.,  those 
of  -~,   ~~j^>  &c*  ma7  ^e  f°lind?  the  right-hand  member 
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being  calculated  by  the  formula  of  the  last  article.    By  equat- 
ing the  parts  independent  of  </>,  we  obtain 


The  value  of  -j-?  may  be  found  by  differentiating  the 
expression  for  Dk  in  Art.  55,  and  substituting  for  —^  and 

da 

^.k+l  their  values  as  given  by  the  present  article. 

The  successive  differential  coefficients  of  Ck  and  Dk  with 
respect  to  a  may  be  obtained  from  the  expressions  for  -~- 

and  —  ~  by  simple  differentiation  and  substitution. 

57.  We  might  determine  in  the  same  way  the  successive 
differential  coefficients  of  Ck  and  Dk  with  respect  to  a  ;  but 
when  those  with  respect  to  a  have  been  found,  the  former 
may  be  derived  from  them,  as  we  proceed  to  shew.  On  ex- 
amining the  expansion  of  (a2  +  a'*—  2aa  cos  </>)"*  in  Art.  52, 
it  will  be  seen  that  Ak  is  a  homogeneous  function  of  a  and  a 
of  —  2s  dimensions.  Hence  Ck  and  Dk  are  homogeneous 
functions  of  a  and  a',  the  former  of  —1,  the  latter  of  —3 

dimensions.     It  follows  that  —~  ,  -^—  will  be  homogeneous, 

functions  of  —  2  and  —  4  dimensions  respectively  ;  and  so  on. 
Now  by  a  known  property  of  such  functions 


dCk       ,  dCk 
a  -T-  +  a  ~ 
da  da 


which  determines  -^ : 
da 
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which  determines 


, 

dor          da  da'  da 

,      *  : 
da  da 


-, 

da2          da  da  da  ' 

d2C 

which  determines     ,  ,2*  :  and  thus  all  the  differential  coeffi- 
da 

cients  of  Ck  may  be  determined. 

In  like  manner  all  the  successive  differential  coefficients  of 
Dk  may  be  calculated. 

We  are  now  in  a  position  to  simplify  the  expression  for 
F.     We  have  (Art.  47.) 


?K 


da'* 


+  j  aa'D^  tan  /  tan  i'  cos  (O  -  flf)  +  .  .  .!  . 

The  following  proposition  will  be  found  useful. 


58.     To  shew  that  ->  ,  =  -  D.,  an^  *Aa*  y—  ^  =  -  D0. 
da  da  da  da 

We  have 

-  (70+  (7,008^)+  <72cos2<£+  ...  =  (a2  +  a'2  -  2aa'  cos  0)~*  ; 

therefore          -  -^  +  -^  cos  <f>  +  -^  cos  26  +  .  .  . 
2  da       c?a  da 

=  —  (a  —  a'  cos  <£)  (a2  +  a'2  —  2aa'  cos  </>)~^J 
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therefore 


+  3  (a  —  a  cos  <£)  (a'  —  a  cos  0)  (a2  +  a'2  —  2aa'  cos  <£)" 

=  cos  (f)  (a2  -f  a'2  —  2aa'  cos  <£)"* 
3  [ad  (1  +  cos2  <£)  -  (a2  +  a'2)  cos  0}  (a2  +  a'2  -  2aa'  cos 

=  cos  <j>  (a2  +  dz  —  2aa  cos  </>)^ 
3  {aa  sin2  </>  —  cos  <j>  (a2  +  a'2  —  2aa'  cos  ^>)} 

(a2  +  a'2  -  2aa  cos  </>)~ 
=  —  2  cos  <j>  (a2  +  a'2  —  2aa  cos  </>)"^ 

sin2  ^>  (a2  +  a'2  -  2aa'  cos  c^)-1- 


Now*     (a2  +  a'2  -  2aa'  cos  ^)"l  =  1  J)0  +  D^  cos  ^  +  ...  ; 

differentiating  with  respect  to  <£ 

3  ad  sin  <j>  (a2  +  a'2  —  2aa'  cos  <£)~^ 

=  D  sin  <>  +  2D  sin  2>  +  ... 


therefore        -j — ^  +  -= — f-,  cos  6  +  ... 
o  .7,.  j,.'     fa  da1 

=  —  2  cos  ^  [-  D0  +  DI  cos 


-f  sin  </>  (Z)x  sin  <£  +  2Z>2  sin  2<^>  +  ..»), 

whence,  equating  the  parts  independent  of  ^>,  and  also  the 
coefficients  of  cos  <j> 


2dada'~         1       2  "        2  ' 


o 

or  ,    /, 

da  da 


and  ,      *, 

da  da 
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59.     Since  -y-5-  is  a  homogeneous  function  of  a  and  a  of 
da 

—  2  dimensions, 


da2  da  da  da  ' 


^      f  dC0  ,  a*d*C0        aa' 

therefore  a  —r-^  +  —  -  -7  2°  = 

da       2    oU2  2 


a 


,       0  0 

Sim  July,        a5^  +  -^=5«a'l>1. 

Hence  the  coefficients  of  e2  and  ez  in  the  expression  for  F 

are  each  equal  to  -  aa'D  . 

8 

Again,  since  Ct  is  a  homogeneous  function  of  a  and  a  of 
—  1  dimensions, 


^     ,i         ^ 

a  -T-1  +  a  -T-  ;  *  =  -  (7.  ; 
aa  aa 

hence  the  coefficient  of  eef  cos  («•  —  -cr') 


4 

but  (Art.  55)          2k  Ck  =  aa'  (D^  -  Dk+l) ; 

therefore,  making  k=l,    . 

hence  the  coefficient  of  eef  cos  («r  —  -nr') 

-  -  aa     2. 
Again,  (Art.  56) 
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therefore  a?DQ  +  a  -j-°  =  aa  D^  . 


Similarly,  a'2  D0  +*'jf**  aa'  A  • 

Hence  the  coefficients  of  tan2*  and  tanV  are  each  equal 

to  -\aa'D^ 
8 

Finally,  the  expression  for  F  becomes 

F=  m'  {§»  +  1  aa  D,  (e2  +  e'2)  -  \  aa  D^  ee  cos  (tr  -  w') 
[28  4 

-     aa'Dj  (tan2  i  +  tan2  1")  +     aa  ^  tan  t  tan  i'  cos  (£1  -  flt') 


CHAPTER  IV. 

SECULAR  VARIATIONS   OF  THE   ELEMENTS   OF  THE   ORBIT. 
STABILITY   OF  THE   PLANETARY  SYSTEM. 


60.  WE  have  seen  in  the  preceding  Chapter,  that  the 
disturbing  function,  when  developed,  consists  of  two  parts ; 
the  one   independent  of  the   time   explicitly,  the  other  in- 
volving it  under  a  periodical  form:  we  shall  consider  sepa- 
rately the  effects  of  these  two  parts.     In  the  present  Chapter 
our  attention  will  be  directed  to  the  first  or  non-periodical 
part  of  jR,  which  we  have  denoted  by  F.     The  inequalities 
thus  produced  in  the  elements  of  the  orbit  are  termed  secular, 
in  consequence  of  their  very  slow  variation. 

61.  By  differentiating  the  expression  for  F  in  Art.  59, 
with  respect  to  the  elements,  we  obtain 

^-0 

d€~    °' 

dF     m       ,_.      ,    .    ,  ,N 

—  —--aaD^ee  sm  (sr  —  OT  ), 

dF     m       ,-„        m       ,-T.    ,        ,  « 

i-  =       aaD.e aa  2JZ  e  cos  (or  —  or  ), 

de       4  4 

<®L  =  _  ^  aa'j)  tan  i tan/  sin  (fl  -  O'), 
x 
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^=  -  ^  aa'Dl  tan  i  +  ~  aa'D^  tan  i'  cos  (£1  -  fl'), 

dF 

-r-  =  an  expression  similar  to  If. 

Substituting  these  in  the  formulas  of  Art.  39,  and  neg- 
lecting small  quantities  of  orders  higher  than  the  second, 
we  have 

da 

3T°' 

de          m'ncfa         ,   .     .  ,. 


,  r.          ^    ,        ,  ,.  , 

{    *e~  D*e™*  (*-«>}}> 

di     m'ncfd 


m'ncf 


.  -^  fl       .  ..        ,~      ^/N-, 

tan  z  -T-  =  --       —  Z\  (tan  ^  —  tan  %  cos  (H  -  H  )} 
wt  ^/^ 

^  =  A  +  A^  (e*  -  tan2  »  )  +  ^8  (e'2  -  tan2  *  ') 

e'  cos  (CT  —  txr')  +  A4  tan  ^tan  i"  cos  (H  —  H'), 


where  in  the  last  expression,  A,  A19  &c.,  have  been  written 
to  denote  certain  functions  of  a  and  a. 

62.  To  calculate  approximately  the  secular  variations  of 
the  elements  of  a  planet's  orbit  ,  in  a  given  time. 

Let  a09  e0,  tsr0,  &c.,  be  the  values  of  the  elements  at  some 
given  epoch  ;  a0  +  Sa*,  e0  +  Se,  OTO  +  SOT,  &c.,  their  values  after 
an  interval  t  :  then  Sa,  £e,  SCT,  &c.,  are  the  required  variations. 
By  Maclaurin's  Theorem, 

*  It  will  be  shewn  in  Art.  64  that  Sa  is  always  zero. 
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u>  +  ••* 


which  may  be  carried  to  any  required  degree  of  accuracy,  but 
in  practice  the  first  two  terms  will  generally  be  sufficient. 

We  have  supposed  the  variations  of  the  elements  required 
at  a  time  t  after  the  epoch  ;  if  they  be  required  at  a  time  t  'be- 
fore the  epoch,  we  have  only  to  change  the  sign  of  t  in  the 
above. 

We  may  remark  that  (  y  )  ,    (~y~)  >  &Ct  are  °f  *ne  order  of 

the  disturbing  force,  since  they  involve  the  first  power  of  m  ; 

fd2e\      fdz^\      0  -I   i         ?    i 

I  -Ta  )  ,     -7^  )  ,  &c.  will  involve  in  2  and  be  of  the  second 

\at  /0     \  at  /0 

order  ;  and  so  on. 

In  the  short  period  of  one  year  all  terms  after  the  first  may 
be  neglected,  so  that  putting  t=I,  we  have 


Hence  the  coefficient  of  t  in  the  above  formulae  is  called 
the  annual  variation. 

63.  Since  the  elements  of  the  planetary  orbits  are  con- 
tinually changing,  it  will  be  interesting  to  shew  that  the 
dimensions  of  these  orbits,  and  their  inclinations  to  the 
ecliptic,  nevertheless  fluctuate  between  very  narrow  limits. 
This  constitutes  what  is  termed  the  Stability  of  the  Planetary 
System  :  in  order  to  establish  it,  it  will  be  necessary  to  prove 
the  stability  (i)  of  the  mean  distances,  (ii)  of  the  excentri^ 
cities,  (iii)  of  the  inclinations. 
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64.  To  prove  the  stability  of  the  mean  distance  of  the  j- 
planets  from  the  Sun,  and  of  their  mean  motions. 

By  Art.  61  —  =  0,  so  that  a  is  constant.  .Now  it  will  be 

shewn  in  a  subsequent  chapter,  (see  Art.  91),  that  to  the  first 
order  of  the  disturbing  force,  the  periodical  terms  of  R  can 
produce  only  periodical  variations*;  consequently  to  this 
order,  the  mean  distance  is  susceptible  of  no  permanent 
change.  The  same  is  true  of  the  mean  motion  n,  since 

Jk 

it  =  —  i  ,  and  p  does  not  alter.    We  are  hereby  assured  of  the 

a'2 

impossibility  of  any  of  the  bodies  of  our  system  ever  leaving  it, 
in  consequence  of  the  disturbances  it  may  experience  from  the 
other  bodies  ;  and  this  secures  the  general  permanence  of  the 
whole,  by  keeping  the  mean  distances  and  periodic  times  per- 
petually fluctuating  between  certain  limits  (very  restricted 
ones)  which  they  can  never  exceed  or  fall  short  of. 

This  result  may  easily  be  extended  to  all  orders  of  the  ex- 
centricities  and  inclinations  :  for  since  nt  +  e  always  occurs  in 
R  as  one  symbol,  e  cannot  occur  in  F  because  t  does  not,  so 

,       dF  f      da. 

that  —  ,  and  therefore  -r-  is  zero. 
de  dt 

65.  To  prove  the  stability  of  the  excentricities  of  the  pla- 
netary orbits. 

We  will  first  consider  the  case  of  two  planets  only.  By 
Art.  Gl, 

de         m'naza  .  . 


*  This  result  is  also  true  when  the  square  of  the  disturbing  force  is  included: 
for  the  demonstration  the  reader  is  referred  to  Ponte'coulant's  Systems  du  Monde, 
or  to  Laplace's  Mecanigue  Celeste. 
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0.    .,    ,      dv          mria'*a  ^  ,     .         ,  < 

Similarly,  -j-  —  --  -  -  Z>2  e  sin  («r  —  w). 

ut  ^/^ 

Now  since  Z>2  is'  the  coefficient  of  cos  2<£  in  the  develop- 
ment of 

(a2  +  a'2  —  2aa'  cos  <£)"*, 

an  expression  in  which  a  and  a'  are  similarly  involved,  it 
follows  that 

A'  =  A- 

Hence,  multiplying  the  above  equations  by  —  e,  —?-,  e, 

na          n  a 

respectively,  and  adding,  we  have 

m     de      m     ,  de 

_____  _  P  _    I    _  __  p    _  _  _  r\  • 

na     dt     rid      dt  ~ 
therefore,  since  a  experiences  no  secular  variation, 


na  na 


A  similar  equation  holds  for  any  number  of  planets.     Ee- 
placing  for  convenience  —  —  -  by  (a,  a'),  we  have 

dp 

—  -  =  —  m'na  (a,  a')  e  sin  («r  —  -cr') 

Cut 

—  m'na  (a,  a")  e"  sin  (-sr  —  -BT")  —  ... 
ju» 

t/c/  f     /    /     f          \  *         /       '  \ 

_  _  =  —  ^^  a  (a  j  a)  e  sin  f«r  —  w) 
a£ 

n    r    t  /    t        n\      tt     '       ft  ti\ 

—  mna  (a,  a  )  e   sin  ^«r  —  -cr  )  —  .  .  . 

——  =  —  mri'a"  (a",  a)  e  sin  («•"  —  in-) 
etc 

t    n    n  r    it        i\     t     *       r     tt  t\ 

—  mna   (a  ,  a  )  e  sin  (&  —  ty  )  —  .  .  . 


Since  D2'  =  J92,  it  follows  that  (a,  a')  =  (af,  a) :  hence  mul- 
ying    1 
we  obtain 


tiplying    these    equations  by  —  e,  -r-,  e,  &c.,  and  adding, 

na      n  a 
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m     de       m     ,  de       m"  ,  „  de" 


whence  by  integration 


\na 

Now  observation  shews  that  all  the  planets  revolve  round 
the  sun  in  the  same  direction,  so  that  the  mean  motions  n,  n ', 
n",  &c.  are  of  uniform  sign.  Hence  all  the  terms  of  the  left- 
hand  member  of  the  above  equation  are  positive. 

We  learn  also  from  observation  that  the  excentricities  of 
the  planetary  orbits  are  at  present  very  small  indeed,  with  the 
exception  of  the  Asteroids,  the  masses  of  which  are  very 
small.  Hence  the  constant  must  be  small.  Since,  then,  all 
the  terms  of  the  first  side  of  the  equation  are  positive,  and 
their  sum  always  equals  a  small  constant;  it  follows  that 
every  term  is  small,  and  therefore  that  the  excentricities  are 
always  small*. 

66.  To  prove  the  stability  of  the  inclinations  of  the  planes 
of  the  planetary  orbits. 


By  Art.  61  .  n  tan  f  sin  (H  - fl'). 

at         4/* 

a.    .,    ,       di'     mricPa  n,         .  .     /rv      _. 
Similarly,  -^  =  — D;  tan  i  sin  (fl  —  O). 

As  in  the  last  article,  it  may  be  shewn  tnat  D't  =  Dt . 
Hence,  multiplying  the  above  equations  by 

m  m'          ., 

-  tan  i ,  — —  tan  i , 
na  na 

*  It  should  be  noticed  that  the  above  is  satisfactory  only  for  those  planets 
whose  masses  are  considerable,  which  is  the  case  with  Jupiter  and  Saturn ;  but 
the  stability  of  the  excentricities  is  not  confined  to  these  planets.  For  a  com- 
plete discussion  of  the  subject  the  reader  is  referred  to  Pontecoulant's  Syst&me 
du  Monde,  or  to  Laplace's  Mecanique  Celeste. 

C.P.T  5 
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respectively,  and  adding,  we  have 

m         .di      m'          .,  di' 

—  tan  i  -j  +  -r-,  tan  i  -j-  =  0, 

na          dt     na  dt 

or  to  the  same  order  of  approximation, 

m         .  d  (tan  {]       m    ,       ..  d  (tan  i'} 
-  tan  i  -±-3  —  -  +  —  tan  i  -^—j-  —  -  =  0  ; 
na  dt          na  at 

i 

therefore  —  tan2  i  +  -j—,  tanV  =  C. 

na  na 

A  similar  equation  would  (as  in  the  case  of  the  excentri- 
cities)  be  true  for  any  number  of  planets.  Now  the  inclina- 
tions of  the  planetary  orbits  to  the  ecliptic  are  at  present  very 
small  :  hence  if  we  take  for  our  fixed  plane  of  reference  a 
plane  coinciding  with  the  present  position  of  the  ecliptic,  it 
follows,  as  in  Art.  65,  that  their  inclinations  to  this  plane  must 
always  remain  very  small*. 

67.     The  stability  of  the  excentricities  and  inclinations 
may  also  be  established  as  follows. 
By  conservation  of  areas 

2  (mh  cos  i}  =  const., 
or  since 

A2  =  fM  (1  -  e2)  =  (M  +  m}  a(l-  e2), 

if  M  denote  the  mass  of  the  sun,  we  have 


m2 
Since  a  is  constant,  if  we  neglect  -7-77,  and  the  fourth 

>y  Jj-L 

powers  of  the  excentricities   and  inclinations,  this  may  be 
written 


*  We  may  remark  that  the  above  demonstration,  like  that  of  the  preceding 
article,  is  applicable  only  to  the  case  of  planets  of  considerable  mass. 
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/ 

or  to  the  same  order  of  approximation 

2  (m  V«e2)  +  2  (m  »Ja  tan2  i)  =  C. 

Since  we  know  from  observation  that  all  the  planets  re- 
volve round  the  sun  in  the  same  direction,  all  the  radicals 
in  this  equation  must  be  taken  with  the  same  sign.  Also, 
since  the  excentricities  and  inclinations  are  at  present  very 
small,  the  constant  must  be  small.  Hence  it  follows,  as  in 
Arts.  65  and  66,  that  the  excentricities  and  inclinations  must 
always  remain  very  small. 

68.  It  may  be  observed  that  the  result  of  the  preceding 
article  proves  the  stability  of  the  excentricities  and  inclina- 
tions as  far  as  the  third  order  of  small  quantities,  while  in 
Arts.  65  and  66  it  was  only  established  to  the  second  order. 
We  will  now  shew  that  if  small  quantities  of  orders  higher 
than  the  second  be  neglected,  the  equation  of  the  preceding 
article  includes  those  of  Arts.  65  and  66. 

On  referring  to  Art.  61,  it  will  be  seen  that  to  the  second 
order,  the  excentricities  and  inclinations  are  given  by  equa- 
tions independent  the  one  of  the  other.  Each  must  therefore 
be  the  same  as  if  the  other  did  not  exist.  Hence  in  the 
equation  of  the  preceding  article,  if  we  make  successively 
i  =  0  and  e  =  0,  we  have 

2  (m  yW2)  =  0,     2(m  V«  tan2  *)  =  <?; 

/      VA& 

or,  since  y  a  =  —  -  , 
na 


na  na 

which  aree  with  Arts.  65  and  66. 


69.  From  the  results  of  the  preceding  articles,  we  draw 
the  following  remarkable  conclusion  :  The  fact  that  the  planets 
revolve  about  the  sun  in  the  same  direction,  ensures  the  stability 

5—2 
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i 

of  the  planetary  system.  The  converse  of  this  would  not 
necessarily  be  true,  as  we  shall  see  in  Art.  75  :  the  numerical 
relations  of  the  dimensions  and  positions  of  the  orbits  of  the 
planets,  might  be  such  as  to  ensure  stability,  although  they 
revolved  in  opposite  directions.  But  the  above  is  independent 
of  particular  numerical  relations. 

70.  The  conclusions  at  which  we  have  arrived  with 
regard  to  the  stability  of  the  planetary  system  are  of  especial 
interest.  In  consequence  of  the  changes  in  the  elements  it 
might  have  been  supposed  that  the  orbits  would  ultimately 
undergo  such  alterations  in  their  dimensions  as  to  bring  the 
planets  into  collision  or  hurry  them  into  boundless  space. 
Or  even  if  no  such  violent  catastrophe  occurred,  a  derange- 
ment of  the  seasons  might  seriously  have  interfered  with  the 
physical  comfort  of  man*.  But  analysis  shews  (and  the 
results  are  confirmed  when  the  approximation  is  carried  fur- 
ther,) that  unless  some  unknown  cause  should  operate  to 
produce  the  contrary,  the  dimensions  and  position  of  the  orbits 
will  for  ages  remain  nearly  the  same  as  they  are  at  present, 
i.  e.  nearly  circular  in  form,  and  but  little  inclined  to  each 
other,  thus  affording  a  beautiful  illustration  of  Gen.  viii.  22 : 
"  While  the  earth  remaineth,  seed-time  and  harvest^  and  cold 
and  heat,  and  summer  and  winter,  and  day  and  night  shall 

*  See  Herschel's  Outlines  of  Astronomy. 


CHAPTER    V. 

SECULAR  VARIATIONS  OF  THE   ELEMENTS   CONTINUED. 
INTEGRATION   OF  THE   DIFFERENTIAL  EQUATIONS. 

71.  IN  Art.  62  we  have  given  a  method  of  calculating 
the  secular  variations  sufficiently  accurate  for  the  practical 
purposes  of  astronomy,  but  the  equations  of  Art.  61  admit  of 
actual  integration.     To  this  we  now  proceed. 

iLa^ 

72.  To  integrate  the  equations  for  the  excentricity  and 
longitude  of  perihelion. 

We  have  (Art.  61)  for  the  planet  m 

de         m'na?a'        ,   .  ,. 


with  similar  equations  for  the  planet  m. 

We  shall  be  able  to  reduce  these  to  a  system  of  linear 
differential  equations  if  we  assume 

u  =  e  sin  -cr,          v  —  e  cos  CT, 
u  —  e  sin  -sr',         v  —  e  cos  -BT'  ; 

.       ,.  du  dor       .       de 

therefore  = 
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Substituting  the  values  of  -^  and  -^  ,  and  writing  a 

,,      m'ncfa  , 

tor  —  -  -     we  nave 
4/* 

--  =  a  (Z^e  cos  v?  —  D2e  cos  w') 


Similarly,          ~  =  a  (Z>X  -  .Zty). 

In  like  manner  for  the  planet  m,  writing  a'  for 

, 
we  nave 


The  forms  of  these  equations  suggest  the  following  par- 
ticular integrals  : 

u=Msm(gt  +  y),         v  =  M  cos  (gt  +  7), 
u  =  Jf  'sin  (#£  +  7)  ,        v  —M'  cos  (#£  +  7)  . 

Substituting  these  in  the  differential  equations,  we  obtain 
from  either  of  the  first  two 


and  from  either  of  the  last  two 

gM1  =  a!  (D^M1  -  DJM)  -, 

eliminating  the  ratio  M  :  M' 


or        - 


and  the  roots  of  this  equation  will  be  real  and  unequal,  real 
and  equal,  or  impossible,  according  as 

(a  +  a')'  A'  -*««'(*»'  -A1) 
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is  positive,  zero,  or  negative.     Now 

(a  +  a')2  A2  ~  *««'  (D?  ~  A2)  =  («  ~  a?  A* 
a  positive  quantity,  since  n,  ri  and  therefore  a,  a'  are  of  like 
sign.  Hence  the  values  of  g  will  be  real  arid  unequal  :  denote 
them  by  ft,  #2  ,  and  let  %,  72  5  ^i>  ^J  ^',  M*  5  he  tlie  cor- 
responding values  of  7,  M,  M'  respectively.  Then  the  com- 
plete solution  of  the  differential  equations  will  be 

u  =  Ml  sin  (ft*  +  %)  +  M2  sin  (gf  +  72), 
v  =  J/x  cos  (ft*  +  7X)  +  Mz  cos  (^/  +  72)5 
u  =  7^'  sin  (^«  +  yj  +  Jf,'  sin  (^  +  7,), 
«'  =  If/cos  (gj  +  7j  +  Jf2'  cos  (gzt  +  7,). 

Of  the  constants  in  these  equations,  four  are  arbitrary  and 
must  be  determined  from  observation.     We  have 


cos  {(ffl  -gz)t  +  7l  - 


tan  «•  =  u  = 


,  cos  (^*  +  7,)  +  M%  cos  (^  +  72)  ' 
with  similar  equations  for  e  and  tzr'. 

73.     Had  we  considered  a  system  of  several  planets,  we 
should  have  obtained  by  a  similar  process 


cos      r  - 
=  M  sin<  Jsin 


with  similar  equations  for  each  of  the  other  planets. 

74  From  the  form  of  the  expression  for  e  in  Art.  72,  the 
stability  of  the  excentricities,  in  the  case  of  two  planets,  may 
]be  inferred.  We  have 
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Consequently  the  excentricity  fluctuates  between  the  limits 
M^  +  Mz  and  M^-M^,  and  since  we  know  from  observation 
that  MI  and M2  are  very  small*,  it  follows  that  the  excentri- 
city must  always  remain  very  small. 

The  period  of  the  changes  in  the  excentricities  =    2?r    , 

and  is  the  same  for  each  planet.  In  the  case  of  Jupiter  and 
Saturn  this  amounts  to  70414  years  !  The  greatest  and  least 
excentricities  which  Jupiter's  orbit  can  attain  are  '06036  and 
•02606,  those  of  Saturn  '08409  and  '01345;  the  maximum  of 
each  excentricity  taking  place  at  the  time  of  the  minimum  of 
the  other.  This  follows  from  the  equation 


na         na 
which  has  been  obtained  in  Art.  65. 

75.  It  appears  from  the  preceding  article  that  the  stabi- 
lity of  the  excentricities  is  a  consequence  of  the  periodical 
form  of  the  solution  of  the  differential  equations,  a  result 
which  depends  upon  the  fact  that  ffl  and  #2  are  real  and  un- 
equal. Now  we  have  seen  that  in  order  that  this  may  be  the 
case,  it  is  only  necessary  that 


shall  be  positive,  a  condition  which  might  be  satisfied  if  the 
signs  of  n,  ri,  and  therefore  of  a,  a'  were  different.  In  this  case, 
then,  the  stability  would  still  subsist.  Let  us  however  con- 
sider what  would  be  the  effect  of  equal  or  impossible  roots  to 
the  quadratic  from  which  g  is  found.  In  the  former  case  a 
term  would  be  introduced  into  u,  u,  v,  and  v  proportional  to 

*  In  the  case  of  Jupiter  and  Saturn,  Sir  John  Herschel  finds  that 
Mv  =  -'01715,  Ma  ='04321,  for  Jupiter; 

MI=     "04877,  J)if2'  =  -03532,  for  Saturn.- 

the  year  1700  being  taken  as  the  epoch.  See  Article,  Physical  Astronomy  in 
the  Encyclopaedia  Metropolitans. 
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the  time,  and  in  the  latter  the  periodical  terms  would  be  re- 
placed by  exponentials.  Consequently  the  excentricities 
would  increase  indefinitely  with  the  time,  and  the  stability 
would  no  longer  subsist. 

76.  We  now  proceed  to  examine  the  expression  which 
has  been  obtained  in  Art.  72,  for  the  longitude  of  perihe- 
lion, viz. 

tan  w  =  MI  sin  (gf  +  7,)  +  M9  sin 

M,  cos  (gf  +  %)  +  MI  cos  (gf  + 


f  2^1/2  cos  {(gl  —  ff2)  t  -f  7X  —  72J 
The  maxima  and  minima  values  of  w,  if  such  exist,  will  be 

found  by  equating  -=-  to  zero.     Thus 
cut 


If  this  (disregarding  sign)  be  not  greater  than  unity,  the 
perihelion  will  oscillate,  the  period  of  a  complete  oscillation 

being  the  same  as  that  of  the  excentricities,  viz.  -    —  ;  but 

n  ~  n 
ffi    fft 

if,  as  is  the  case  with  Jupiter  and  Saturn,  this  be  greater 
than  unity,  the  longitude  of  perihelion  has  no  maximum  or 
minimum,  and  the  perihelion  moves  constantly  in  one  direc- 
tion. 


Again, 


, 

2 


2e2  2 

Hence  when  e  is  a  maximum  or  minimum,  —j-  will  be 

at 

either  a  maximum  or  minimum,  and  the  apsidal  line  will  be 
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moving  most  rapidly  or  most  slowly,  different  cases  occur- 
ring according  to  the  signs  and  magnitudes  of  the  quantities 
involved. 

.  77.      When  the  apsidal  line  oscillates,  to  find  the  extent 

and  periods  of  its  oscillations. 

We  have  (Art.  72) 
tantsr 


sin  (9f  +  %)  +  ^  am  (gf  +  7.)  . 


i  cos  (&*  +  7i)  +  K  cos  (9f  +  7.)  ' 

,       r  ,  s        tan-zzr  —  tan  (g.t  +  y.) 

therefore    tan  [m  —  gj  —  7  ,)  =  -  ^^x  - 

1  +  tan  -57  tan  (gf  +  7^ 


—       sn 


x  +       2  COS 

if       =        - 


Also  by  the  last  article,  if  T  be  the  least  positive  angle 

,  .      .          <2  2 

whose  cosine  is  —  -f 


~dt       'l  l2  ~  °OS  T 

Different  cases  will  occur  according  to  the  signs  of  J/J, 
M^  &c.  Suppose  Mv  Mz  to  have  different  signs,  g^  and  ^ 
positive,  and  g^  greater  than  g2.  Then  ty  increases  as  t  in- 

creases, and  -J-  will  be  negative,  or  the  apsidal  line  will 
regrede,  while  cos  ty  —  cos  r  is  positive,  i.  e.  so  long  as  ^  is 
between  2mr  —  T  and  2n7r  +  T  :  —  T-  will  be  positive,  or  the 

Cut 

apsidal  line  will  progrede,  while  ^  is  between  2n7r  +  r  and 

2  (fl  +  1)  7T-T. 

To  find  the  angle  through  which  the  apsidal  line  regredes 
and  the  period  of  the  regression.  Let  t',  t"  be  the  values 
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of  t,    «r',  -cr"  the  values  of  «r  corresponding  to  the  values 
27Z7T  -  T  and  2w-7r  +  r  of    r  :  then 


~  72  = 

sn  T 


.  9 

tan  (w  —  ?J  —  %)  =  -=7  —  ^7  -  » 
Ifj  +  M,2  cos  T 


/   fi        »       \ 

tan  (  37  —  q.t  —  yj  =  ^  -  ~  —    —  . 
MI  +  Jf2  cos  r 

From  these  equations  the  values  of  t',  t",  vr',  vr"  may  be 
found,  and  thus  -57'  —  w"  the  amount  of  regression  will  be 
known.  The  period  of  regression 


In  like  manner  the  amount  and  period  of  the  progression 
may  be  obtained.     The  latter  will  be  found  to  be  -  —  —  • 
The  period  of  a  complete  oscillation  will  be  the  sum  of 

27T 

the  periods  of  the  regression  and  progression,  that  is  —  —  > 

9\     9z 
which  agrees  with  the  preceding  article. 

78.  The  motion  of  the  centre  of  the  instantaneous  ellipse 
in  consequence  of  the  secular  variations  of  e  and  -cr  may  be 
exhibited  geometrically  as  follows. 

We  have,  by  Art.  72, 

e  cos  «r  =  Jfj  cos  (gf  +  yt)  +  Mz  cos  (gj  +  ya), 
e  sin  tar  =  M^  sin  (gf  +  yj  +  ^T2  sin  (^a<  +  y2)  . 

Let  a  circle  be  described  in  the  plane  of  the  orbit  with 
its  centre  8  coinciding  with  that  of  the  sun,  and  its  radius 
equal  to  J^a,  where  a  is  the  mean  distance.  Let  a  point  P 
describe  this  circle  uniformly  with  a  velocity  g^  starting  from 


76 
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0.  Again,  with  centre  P  and  radius  equal  to  Mza  let  another 
circle  be  described,  and  let  a  point  Q  describe  this  circle 
uniformly  with  a  velocity  g^  starting  from  C.  Let  SL  be 


the  line  from  which  longitudes  are  reckoned,  and  draw  PK 
parallel  to  it:  then  if  the  angle  OSN  be  equal  to  71?  and 
CPK  to  72,  the  angle  PSN  will  be  equal  to  ffj  +  v^  and 
QPK  to  g2t  +  72.  Produce  QP  to  meet  the  circle  again  in  E, 
and  draw  QN  perpendicular  to  SN.  Then,  supposing  Mt  and 
Ma  to  be  both  positive,  we  have 

=  SPcos  PSN+  PQ  cos  QPK 

(gj  +  7X)  +  M^a  cos  (gj  +  72) 


=  ae  cos  iy. 

Similarly,  it  may  be  shewn  that 
QN  =  ae  sin  -or. 

Hence,  the  apse  being  supposed  to  move  from  L  in  the 
direction  contrary  to  that  of  the  hands  of  a  watch,  Q  will  be 
the  centre  of  the  instantaneous  ellipse. 
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If  Ml  be  positive  and  Mz  negative,  it  may  be  shewn  in 
like  manner 'that  the  centre  of  the  ellipse  will  be  E.  If 
Mlt  M2  be  both  negative,  join  QS  and  produce  it  to  Q'  so 
that  SQ'  =  SQ:  then  the  centre  of  the  ellipse  will  be  Q'. 

A  similar  construction  will  of  course  apply  for  the  motion 
of  the  further  focus. 

79.  To  integrate  the  equations  for  the  inclination  and  lon- 
gitude of  the  node. 

We  have  (Art.  61)  for  the  planet  m 

7«       vn/nn^n 

j-  =  — D:  tan  i'  sin  (ft  -  ft'), 

» 

tan  i-jjj-  = —  D^  (tan  i  -  tan  i'  cos  (ft  —  ft')}  ; 

with  similar  equations  for  the  planet  m'. 
To  integrate  these,  assume 

p  =  tan  i  sin  ft,     q  =  tan  i  cos  ft, 
p  =  tan  i'  sin  ft',    q  =  tan  i'  cos  ft' ; 

therefore  -f-  =  tan  i  cos  ft  -T-  +  sin  ft  (1  +  tan2  i}  -y  . 
at  at  at 

Substituting  the  expressions  for  -T-  and  -j ,  and  writing 

cit  a/t 

„      m  na  a       .  „  ,di ,    .         „,      i.-i       -i 

a  lor  — ,  since  tan  i  -y  being  of  the  third  order  may  be 

omitted,  we  have 

-y-  =  ccZ^  (tan  i'  cos  ft'  —  tan  i  cos  ft) 


Similarly, 
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.  ,      r      ,       ,  ,        ...         ,  f 

Also  for  the  planet  w,  writing  a  for 


-  5, 


The  forms  of  these  equations  suggest  the  following  par- 
ticular integrals  : 


(Jit  +  8), 
/  =  JV'sin  (^  +  8),'         q=N'cos  (ht  +  8). 

Substituting  these  in  the  differential  equations,  we  obtain 
from  either  of  the  first  two, 


and  from  either  of  the  last  two, 
hN'  =  v!D, 

eliminating  the  ratio  N  :  N'9 


or 

therefore  h  =  -  (a  +  a')  Z\,  or  A  =  0. 

Denote  the  former  by  ^,  and  let  81?  82,  ^,  J^2,  ^',  JV2', 
be  the  values  of  8,  JV,  ^V'  corresponding  to  h  =  \  and  A  =  0. 
Then  NJ  =  N2,  and  the  complete  solution  of  the  differential 
equations  will  be 


p  =     I  sn      t  +       +     2  sn 
q=Nl  cos  (^<  +  8J  +  ^s  cos 
'  =  JV  sin      t  +  S    +  ^  sin 


2  , 
=  N  cos 
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Of  the  constants  in  these  equations,  four  are  arbitrary, 
and  must  be  determined  from  the  known  values  of  i  and  £t  at 
some  given  epoch. 

We  have  then 


tan2  i  =/  +  f  =  A?  +  ^22  +  2AyV2  cos  (JiJ 

_      p     N,  sin  (hj  +  Bt)  +  Na  sin  82 
tan  11  =  —  —  -^7  -  it.  —  ^~~\      AT         s> 
q     Nt  cos  (^  +  oj  +  jy,  cos  02 


with  similar  equations  for  i'  and  H'. 

Had  we  considered  a  system  of  several  planets,  we  should 
have  obtained  a  result  precisely  similar  to  that  of  Art.  73. 

a* 

80.     From  the  form  of  the  expression  for  tan  i,  the  sta-    . 
bility  of  vthe  inclinations,  in  the  case  of  two  planets,  may  be 
inferred.     We  have 

tan2  *  =  N?  +  N*  +  ZN^  cos  ( 


Consequently  tan  i  fluctuates  between  the  limits  N^  +  ^V2  and 
NI  ~  Nz  ;  and  since  we  know  from  observation  that  N^  and  j\T2 
are  very  small,  it  follows  that  the  inclination  must  always 
remain  very  small. 

Further,  the  periods  of  the  changes  in  the  inclinations  of 

the  orbits  of  the  two  planets  are  the  same,  being  -      ;  and 


as  appears  from  the  equation  of  Art.  66,  the  maximum   of 
each  inclination  will  take  place  at  the  time  of  the  minimum 

of  the  other. 

In  the  case  of  Jupiter  and  Saturn,  the  period  is  50673 
years  ;  the  maximum  and  minimum  inclinations  of  Jupiter's 
orbit  to  the  ecliptic  are  2°  2'  30"  and  1°  17'  10",  those  of  Saturn's 
orbit  2°  32'  40"  and  0°  47'. 

81.     We  now  proceed  to  examine  the  expression  which 
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has  been  obtained  in  Art.  79  for  the  longitude  of  the  node. 
We  have 

n  _  ffi  sin  (\t  +  8t)  +  JV2  sin  82 
' 


The  maxima  and  minima  values  of  O,  if  such  exist,  will 
j 

Out 


be  found  by  equating  -j-  to  zero.     Thus 


. 


2 


If  this  (disregarding  sign)  be  not  greater  than  unity,  the 
node  will  oscillate,  the  period  of  a  complete  oscillation  being 

the  same  as  that  of  the  excentricities,  viz.  —  r  .     But  if  it  be 

±h 

greater  than  unity,  there  cannot  be  any  stationary  positions, 
and  the  node  will  move  continually  in  one  direction. 

It  may  be  shewn,  as  in  Art.  76,  that  the  motion  of  the 
node  will  be  fastest  or  slowest  whenever  the  inclination  is 
either  a  maximum  or  minimum. 

82.  When  the  line  of  nodes  oscillates,  to  find  the  extent 
and  periods  of  its  oscillations. 

It  may  be  shewn  as  in  Art.  77,  that  if  -\Jr  be  written  for 
j^-j.^  _  S2)  and  T  denote  the  least  positive  angle  whose 

.      .        N, 
cosine  is  —  Tn?  , 


—  cos  T  sn 

~~  1  —  COS  T  COS  A|T  ' 

and  tan2  i  —  y-  =  h^N^  (cos  ^r  —  cos  T). 

Different  cases  will  occur  according  to  the  signs  of  N^  , 
Nt  and  ht.     Suppose  N19  Nz  of  like  sign,  and  \  negative: 
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then  ^  decreases  as  t  increases,  and  the  line  of  nodes  re- 
gredes  so  long  as  ty  is  between  2/iTr  +  T  and  2w7r  — T,  and 
progredes  so  long  as  ty  is  between  2nir  —  r  and  2  (n  —  1)  TT  +  T. 

Let  ft',  ft"  be  the  values  of  ft  corresponding  to  the  values 
2n7T  +  T  and  %nir  —  T  of  ty ;  then 

tan  (ft'  -  S2)  =  -  cot  T, 
tan(ft"-S2)=cotT; 

7T 

therefore  ft'  —  S2  =  WTT  +  T  —  —  , 


therefore  ft'  -  ft"  =  2r  -  TT, 

which  is  the  angle  through  which  the  line  of  nodes  regredes. 
Also  the  period  of  this  regression  may  be  shewn  as  in  Art.  77 

2r 
to  be  — j- .     Similarly,  the  angle  through  wThich  the  line  of 

nodes  progredes  may  be  shewn  to  be  2r  —  TT,  and  the  period  of 

the  progression  — • — 7 — -  . 
~~  "i 

The  period  of  a  complete  oscillation  will  be  the  sum  of 

rt__ 
the  periods  of  the  regression  and  progression,  that  is  — =-  , 

~~"i 
which  agrees  with  the  preceding  Article. 

The  remaining  cases  corresponding  to  different  arrange- 
ments of  the  signs  of  Nl9  Nz  and  \  may  be  treated  in  like 
manner. 

The  mean  value  of  fl  is  imr  +  B^,  (this  will  be  found  to 
be  the  case  whatever  be  the  signs  of  Nt ,  N^  and  \ ;)  and 
the  mean  value  coincides  with  the  true  whenever  sin  ty  =  0. 
Since  then  ty  is  the  same  both  for  the  disturbed  and  disturb- 
ing planet,  the  nodes  of  both  orbits  will  arrive  simultaneously 
at  their  mean  positions. 

C.P.T.  6 
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In  the  case  of  Jupiter  and  Saturn  N2  is  for  each  planet 
numerically  less  than  N19  so  that  the  node  oscillates;  the 
extent  of  oscillation  being  13°  9'  40"  in  Jupiter's  orbit,  and 
31°  56'  20"  in  that  of  Saturn  on  either  side  of  their  mean 
position,  the  ecliptic  being  taken  for  the  plane  of  reference, 
and  supposed  immoveable. 

83.  To  shew  that  the  inclination  of  the  orbits  of  two 
mutually  disturbing  planets  to  each  other  is  approximately 
constant. 

If  7  denote  this  inclination,  we  have  by  Spherical  Trigo- 
nometry, 

cos  7  =  cos  i  cos  i'  +  sin  i  sin  i'  cos  (XI  —  XI') 

=  cos  i  cos  i'  {1  +  tan  i  tan  i'  cos  (XI  —  XI')} 

=  (1  +  tan2 1')~*  (1  +  tan2 /')"*  (1  +  tarn  tarn"  cos  (XI  -  XI')} 

=  1  -  i  {tan2  i  +  tan2  i'  -  2  tan  /tan  {  cos  (XI  -  XI')}, 

if  we   neglect    small   quantities   of  orders   higher   than  the 
second. 

Now       tan2  i  +  tan2  i'  —  2  tan  i  tan  i'  cos  (XI  —  XI') 
=  (P-P'}*+(1-CL}" 


therefore  1  -  cos  7  =  -  (^  -  A7/)2, 

A 


or  sn     =  ~; 

whence  it  follows  that  7  is  constant. 

84.     The  equations  which  give   the   secular  variations 
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of  the  node  and  inclination  may  be  interpreted  geometrically 
as  follows*. 

The  equations  to  be  interpreted  are 

p  =  Nl  sin  (hj  +  8J  4-  Nz  sin  S2, 
q  =  JV;  cos  (hjt  +  8J  +  JV;  cos  S2, 
where  p  —  tan  i  sin  fl,       t 

q  =  tan  i  cos  H, 

Since  in  the  differential  equations  from  which  these  have 
been  obtained,  small  quantities  of  the  third  order  have  been 
neglected,  we  have  to  the  same  order  of  approximation 

Nt  sin  (hj  +  8J 


sn    sn 


sn 


sin  i  cos  &  =  Nt  cos  (hf  4-  3J  +  ^2  cos 


.  . 


Let  a  sphere  be  described  with  its  centre  coinciding  with 
that  of  the  sun,  and  its  radius  of  any  magnitude:  let  the 
fixed  plane  of  reference  cut  it  in  the  great  circle  AL,  L  being 
the  origin  from  which  longitudes  are  measured,  and  LA  a 

*  For  this  elegant  geometrical  interpretation,  the  Author  is  indebted  to 
Mr  Freeman,  of  St  John's  College.     See  Appendix,  Art.  7. 

6—2 
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quadrant.  Let  another  fixed  plane  inclined  at  a  small  angle 
I  to  the  former  cut  the  sphere  in  NM,  and  let  LN=  a).  Let 
K,  P,  Q  be  the  poles  of  AL,  NM,  and  the  plane  of  the  orbit 
respectively:  join  KA,  KL,  PA,  PL,  QA,  QL,  PQ,  KP,  and 
produce  the  latter  to  meet  AL  in  0.  Then  N  will  be  the 
pole  of  OK,  and  ON  a  quadrant  ;  therefore 


Let  PQ=  p,  and  the  angle   QPO=6.     Then  from  the 
right-angled  triangles  A  OP,  LOP,  we  have 

cos  AP—  sin  /cos  w,  cos  LP=  sin  Jsin  co, 

sin  AP&in  APO  —  sin  o>,     sin  ^IPcos  APO  =  cos  I  cos  w. 

Similarly 

cos  A  Q  =  sin  i  cos  O,         cos  L  Q  =  sin  /sin  H. 

Now  cos  J.  §  =  cos  AP  cos  P$  +  sin  AP  sin  P$  cos  APQ 
=  cos  JUPcos  P§  +  sin  ^LPsin  P§ 

(cos  APO  cos  0P<2  -  sin  APtfsin  OPQ], 

or         sin  i  cos  O  =  sin  /cos  o>  cos  p 

+  sin  p  (cos  /  cos  &>  cos  0  —  sin  &>  sin  6) 
=  sin  /cos  co  cos  p  +  sin  p  cos  (0  +  o>), 

neglecting  sin2  -  /sin  p. 

A 

Similarly, 

sin  i  sin  fl  =  sin  /sin  w  cos  p  +  sin  p  sin  (0  +  co). 

Now  these  equations  will  be  identical  with  equations  (1), 
if  we  suppose 

sin  p=*Nl9  sin  /  cos  p  =  Nz, 
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We  have  then  the  following  interpretation  :  —  the  normal 
to  the,  plane  of  the  orbit  moves  uniformly  so  as  to  generate  in 
space  a  fixed  right  circular  cone. 

v^x 

85.  To  integrate  the  equation  for  the  longitude  of  the 
epoch. 

We  have  (Art.  61) 
^  =  A+A,  (e*  -  tan2  i)  +  A^  (e'z  -  tan2  t  ') 

+  As  ee  cos  (OT  —  -or')  +  A4  tan  i  tan  i'  cos  (fl  —  £t')  . 

Now  from  the  formulae  of  Art.  72,  we  obtain 

e"  =  M?  +  M'  +  2^,  cos  {(9l  -g,}  t+y.-y,}, 


ee  cos  w  -       = 

+  (MM  +  MJ®  cos  ((g.  -  g)  t  +  7,  - 

In  like  manner,  from  the  formulae  of  Art.  79, 

tan2  i  =  ^2  +  JV22  +  2JV;  N2  cos  fat  +  St  -  8 
tan2  1  '=  AT'2  +  ^2/2  +  2NJNt  cos  {^  +  Sx  - 
tan  i  tan  t  '  cos  (O  -  fl')  =  J^^'  +  ^"22 


If  these  values  be  substituted  in  the  expression  for  ~ 

dt 

it  takes  the  form 

de 
dt 


cos  ((&-&)  *  +  7,  -  yj  +  B2  cos  (^  +  S, -  8 
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where  Bn,  Bt ,  and  B2  denote  certain  constants.     Integrating, 
we  have 

A 

e^+    +_^.      "«  > '^*  ^  ^^ 

We  may  omit  the  term  -Zfotf,  if  we  consider  it  as  furnishing 
a  correction  on  the  mean  motion  n,  which  thus  becomes 
(1  +  B)  n.  With  this  understanding 

B,  ,          jD,, 


If  this  expression  be  developed,  we  may  again  omit  the 
term  involving  the  first  power  of  t,  and  consider  it  as  affording 
a  further  correction  to  the  mean  motion*.  Thus  we  shall  ob- 
tain a  series  of  the  form 


86.  In  the  Theory  of  the  Planets  this  inequality  is  insen- 
sible, but  in  that  of  the  Moon  it  amounts  to  upwards  of  10 
seconds  in  a  century,  forming  what  is  termed  the  secular  ac- 
celeration of  the  Moons  mean  motion.     Thus  it  appears  that 
this  inequality  does  not,  as  its  name  would  seem  to  imply, 
contradict  the  general  theorem  of  the  invariability  of  the  mean 
motions,  since  it  is  due  to  a  variation  not  of  the  mean  motion, 
(as  we  have  employed  the  term  in  the  preceding  pages,)  but 
of  the  epoch.     If  however,  as  in  the  Lunar  Theory,  the  epoch 
be  omitted,  any  variation  in  the  mean  longitude  will  of  neces- 
sity be  thrown  upon  the  mean  motion  ;  only  in  this  case,  n 
will  not  be  given  by  the  equation  n*a*  =  fi. 

87.  We  have  hitherto  supposed  the  planetary  motions  to 
be  referred  to  a  fixed  plane,  but  have  left  the  particular  plane 

*  The  advantage  of  thus  disposing  of  these  terms  arises  from  the  fact  that 
the  mean  motion  as  determined  by  observation  is  the  complete  coefficient  of  t 
in  the  expression  of  the  mean  longitude. 
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undetermined.  In  practice  it  is  usual  to  take  the  position  of 
the  ecliptic  at  some  given  epoch,  as  for  instance  the  year 
1800;  but  since  it  is  to  the  true  ecliptic  that  astronomers  refer 
the  celestial  motions,  we  will  now  obtain  formulas  for  deter- 
mining relatively  to  the  plane  of  the  earth's  orbit,  the  position 
of  that  of  any  other  planet. 

Let  then  m,  m  denote  the  masses  of  the  earth  and  the 
planet  considered,  and  suppose  the  orbits  of  r&  and  m  but 
little  inclined  to  each  other  and  to  the  fixed  plane  of  reference. 
Let  X,  X'  denote  the  latitudes  of  points  in  these  orbits  corre- 
sponding to  the  same  longitude  0l  •  then  (see  fig.  to  Art.  13) 

tan  X  =  tan  i  sin  (0X  -  H),  tan  X'  =  tan  i'  sin  (0X  -  O'). 

Now  since  i  and  i'  are  very  small,  we  may  replace  tan  X, 
tan  X'  by  X,  X'  respectively :  thus 

X'  -  X  =  tan  i'  sin  (^  - H')  -  tan  t  sin  (0t  - 11) 
=  (tan  i'  cos  H'  —  tan  i  cos  H)  sin  0t 

—  (tan  i'  sin  O'  —  tan  t  sin  11)  cos  0t , 
or,  with  the  notation'  of  Art.  79, 

V-X=(2'-2)sin01-(|/-dp)cos01 (1). 

Now  let  7  denote  the  inclination,  v  the  longitude  of  the 
node  of  the  orbit  of  m  relatively  to  that  of  m ;  then  approxi- 
mately 

X'  —  X  =  tan  7  sin  (6^  —  v) 

=  tan  7  cos  v  sin  9^  —  tan  7  sin  v  cos  0± (2). 

Hence  equating  coefficients  of  sin  0A  and  cos  0X  in  equations 
(1)  and  (2), 

tan  7  cos  v  —  Ci—Qi    tan  7  sin  v  =  p'  —  p ; 
whence  tan2  7  =  (p'  —p}*  +  (q  —  #)2, 
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and  io^v  —  —, — -. 

2-2 

These  expressions  determine  the  position  of  the  orbit  of 
m  relatively  to  that  of  ???,,  when  the  values  of  p,  p,  q,  and  q 
are  known.     Differentiating  them,  and  neglecting  small  quan- 
tities of  orders  higher  than  the  second,  we  obtain 

dv      (dp      dp\    .  fdq      dq\ 

-31=(j*  ~~fl}  Sm  v  +  \~jl  —  J*     COS  v> 

,   dt      \dt      dtj  \dt      dtj 

dv  _  (dp      dp\  cos  v      fdq      dq\  sin  v 
di ~  \di ~~ ~di)  tan 7 "  \dt  ~  dt)  tan 7 * 

If  the  values  of  -~ ,  -^- ,  &c.  be  substituted,  these  equations 
give  the  variations  of  7  and  v. 

88.  For  the  theory  of  the  invariable  plane  of  the  solar 
system,  the  reader  is  referred  to  Pontecoulant's  Systeme  du 
Monde. 


CHAPTER     VI. 

PERIODIC  VARIATIONS  OF  THE  ELEMENTS  OF  THE   ORBIT. 

89.  WE  come  now  to  consider  the  variations  produced 
by  the  periodical  terms  of  R.     These  are  called  Periodical 
Variations,  as  opposed  to  the  Secular  Variations  produced  by 
the  non-periodic  terms.     We  have  seen  indeed  that  the  latter 
are  for  the  most  part  periodical  in  form,  but  in  the  Planetary 
Theory,  the  term  Periodical  Variations  is  restricted  to  those 
we  are  about  to  consider  in  the  present  Chapter. 

90.  We  have  seen  (Art.  46)  that  the  general  type  of  a 
periodical  term  is  P  cos  {(pn  ±  qri)  t  +  Q],  where  P  is  a  func- 
tion of  a,  e,  i;  and  Q  is  a  function  of  «r,  e,  fl.      Now  such  a 

.     dR     dR         ,  dR    , 
term  will  produce  a  similar  term  in  -3— ,   -T-  ,  and  -^r ;  out 

a   term   of  the  form  P  sin  {(pn  ±  qri)  t+  Q}   in  ^ ,   -^  , 

/77? 

an(j         .      If  then   these  be  substituted  in  the   equations 
all 

of  Art.  39,  they  will  take  the  forms 


90  PLANETAEY  THEORY. 


-&=?* cos  ((Pn  ±  ?»')  *  +  <2)' 
+  G}, 


where  P1?  P2,  &c.  are  functions  of  the  elements  of  the  dis- 
turbed and  disturbing  planets,  and  involve  the  first  power  of 
the  disturbing  mass. 

91.  In  integrating  these  equations,  we  may  in  general 
consider  the  elements  which  enter  in  the  right-hand  members 
as  constant  and  equal  to  their  values  at  the  epoch  from 
which  the  time  is  reckoned*. 

Let  then  a,  e,  -or,  &c.  denote  the  values  of  the  elements 
at  epoch,  Sa,  Be,  SOT,  &c.  their  periodical  variations  after  an 
interval  t :  then  integrating  the  above  equations 

p 

Ba  — ;  --1 — 7  cos  {(pn  ±  qri)  t  +  Q}, 

pn  Hr  CTfl 
p 

Be  — -* — r  cos  {(pn  ±  qri]  t  +  Q}, 

p 

SOT  = -? — ,  sin  {(pn  ±  qri}  t  +  Q}, 

pn  ±  qn 

*  This  is  equivalent  to  neglecting  the  square  of  the  disturbing  force :  see 
Art.  95. 
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p 


''  =  "  7^+W  C°S  {(PU  ±  ***   *  +  ®' 


Hence   it   appears  that  the  variations  produced  by  the 
periodical  terms  of  It  are  all  periodical  in  form. 

92.  It  will  be  seen  that  all  the  expressions  of  the  last 
Article  involve  the   divisor  pn  ±  qri ,  while  Sf  involves  the 
divisor    (pn  ±  qnf.     If  then  it  should  happen   that  either 
pn  +  qn    or  pn  -  qn'  is  very  small,  a  term  in  R  containing 
(pn  ±  qn)  t  in  its  argument,  though  of  a  high  order,  may 
have  a  sensible  effect  on  the  elements  of  the  orbit.    Now 
since  p  and  q  are  either  positive  integers  or  zero,  pn  +  qn 
cannot  be  small  unless  n  and  n  are  small,  a  case  which  does 
not  occur  with  any  of  the  planets:  but  we  have  instances 
in  which  pn  ~  qn   is  small*. 

Since  the  period  of  such  inequalities  is  very  great  (being 

—, ) ,  they  are  called  inequalities  of  long  period,  or  long 

pn  "**  qn  / 

inequalities. 

93.  To  select  such  terms  in  R  as  will  produce  the  prin- 
cipal inequalities  of  long  period. 

*  If  in  any  case  the  mean  motions  of  two  planets  were  exactly  commen- 
surable and  in  the  ratio  of  p  to  q,  the  corresponding  term  of  R,  as  we  have 
already  remarked  (Art.  48),  would  cease  to  be  periodical  and  would  form  a  part 
of  F,  but  no  instance  of  this  occurs  among  the  planets. 
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We  have  seen  that  the  dimension  of  the  principal  part 
of  the  coefficient  of  a  term  containing  (pn  ~  qri)  t  in  its  argu- 
ment is  p  ~  q  (Art.  50) ;  hence  if  we  can  find  two  integers 
p  and  q  nearly  in  the  ratio  of  n  to  n,  and  having  a  small 
difference,  the  corresponding  term  of  E  will  produce  an  im- 
portant long  inequality  in  the  elements  of  each  planet. 

In  the  case  of  Jupiter  and  Saturn  n  :  n  ::  5  :  2  nearly, 
and  5  —  2  =  3;  hence  there  is  a  long  inequality  arising  from 
a  term  in  R  of  the  form  P  cos{(2/i  —  6ri)  t  +  §},  the  prin- 
cipal part  of  P  being  of  the  third  order.  This  inequality  is 
interesting  in  an  historical  point  of  view,  having  long  baffled 
the  labours  of  mathematicians  and  appeared  inexplicable  on 
the  hypothesis  of  gravitation.  It  was  at  last  successfully 
explained  by  Laplace. 

For  the  Earth  and  Venus,  n  :  n  ::  8  :  13  nearly,  so  that 
there  is  a  long  inequality  arising  from  a  term  in  R  of  the 
fifth  order.  The  discovery  of  this  inequality  is  due  to  the 
Astronomer  Royal. 

Finally,  in  the  case  of  Neptune  and  Uranus,  n  :  n  ::  I  :  2 
nearly,  hence  there  is  a  long  inequality  arising  from  a  term 
in  R  which  is  of  the  first  order. 

94.  Between  corresponding  terms  of  the  long  inequali- 
ties in  the  mean  motions  of  two  planets,  arising  from  the 
near  commensurability  of  n  and  n',  there  is  a  simple  approx- 
imate relation. 

Let  m,  m  be  the  masses  of  the  two  planets,  -5,  R'  their 
disturbing  functions:  then  by  Art.  8,  considering  only  the 
mutual  action  of  m  and  m,  we  have 
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We  shall  distinguish  the  first  and  second  terms  of  R  and 
R  as  the  symmetrical  and  unsymmetrical  parts  respectively, 
since  the  co-ordinates  of  m  and  m  are  involved  symmetrically 
in  the  former  but  not  in  the  latter. 

Since  then  the  symmetrical  parts  of  R  and  R'  differ  only 
in  having  m  and  m'  interchanged,  if 

m'M  cos  {(pn  —  qn}  t  +  Q] 

be  any  term  in  the  symmetrical  part  of  jft,  that  of  R  will  con- 
tain the  term 

mMcos  {(pn  —  qn}  t+  Q}. 

Confining  our  attention  to  these  terms,  we  have  (Art.  39) 

d*$_      3n*adR_      Znad(R) 
dt*  fj.de  jj,      dt 

—^  m'M  sin  {(pn  —  qn}  t  +  Q}  ; 

,T       r  r^         %ri*ap       m'M 

therefore      8f  =      —-        _          sm  {(pn  -  qn')  t+Q}. 


a-    -i    i      **>      3n'2dq       mM         .     ,,  ,.         ^ 

Similarly,  S£  =  —  ^  >    sm  {(pn-q.n}t+  Q}. 


TJ  > 

Hence        ~i  =  -- 

of          w/t/i  a  q 

approximately,  since  qn  is  nearly  equal  to  pn  ;  therefore 


or  since  ^'  differs  from  ^  by  a  quantity  of  the  order  of  the 
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disturbing  force,  the  square  of  which  we  are   neglecting, 

we  have 

5?          m'ija' 

S£~    ~  m  V«  ' 
the  required  relation. 

The  same  relation  is  also  approximately  true  in  the  case  of 
terms  arising  from  the  unsjmmetrical  parts  of  E  and  R'*. 
For  denoting  these  by  J5X  and  E{  respectively,  we  have 

x' 


Now  the  equations  of  motion  of  the  planet  m  referred  to 
rectangular  axes  are 


=      _ 

*  ' 


dt?       r*       dx 

Hence,  the  differential  coefficients  being  taken  as  if  the 
elements  were  constant  f 


_ 

df  ?    r'3~        1  dt* 


m 
therefore          B 


*  For  the  demonstration  of  this  we  are  mainly  indebted  to  The  Theoi'y  of 
the  Long  Inequality  of  Uranus  and  Neptune:  an  essay  which  obtained  the 
Adams  Prize  for  the  year  1850.  By  E.  Pierson,  M.A. 

•j*  This  is  simply  an  analytical  artifice :  we  merely  assert  that  if  the  differen- 

x' 
tial  coefficients  be  so  taken,  then  -^,  &c.,  and  therefore  RI  and  JRi  will  take  the 

above  forms. 
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,     mf,d*xt     ,d*y       ,d* 


,     mf,xt     ,y 
Similarly,       B,  =  -  (x  -^  +  y  w 

Now  any  term  in  E^  containing  (pn  —  qn'}  t  in  its  argu- 
ment can  arise  only  from  the  combination  of  terms  in  #,  yy 


and  z,  containmg^??^  in  their  arguments  with  terms  in  -™-  , 

~  ,  and  -^  containing  qrit.     Suppose  then  x  and  x  when 
dt  cit 

developed  in  terms  of  t  and  the  elements  to  contain  respec- 
tively the  terms 

L  cos  (pnt  +  Z),   L'  cos  (qnt  +  I'). 

dzx 
Hence  the  product  x  -^  will  contain  the  term 

-  i  LL'fn*  cos  {(pn  -  qn}  t  +  1-1'}, 

dzx 
and  the  product  x  --,  2  the  term 


-     LL'  »V  cos  {(jw  -  ^n')  *  +  1  -  Z'}  : 

A 

the  coefficients  are  in  the  ratio  (fri*  to  ^?V.     Similarly,  the 

coefficients  of  the   same  cosine  in  y  —j-»   and  z  -7-5-  are  to 

at  at 

those  in  y'  -—  and  z  -^  in  the  same  ratio. 


Hence  if       —  MqV  cos  {(pn  —  qn}  t  +  Q} 
be  any  term  in  R^  then  E^  will  contain  the  term 
-  J^V  cos  {Q?/i  -  2/0  <  +  Cl- 


96  PLANETARY   THEOEY. 

Confining  our  attention  to  these  terms,  we  have 
s.,         Srfap    m'MqV     . 

•~s 


~^ 

: 


ma  m'na 


TT  , 

Hence  ^  =  --  -?  =  --  —  nearly, 

eg          map          mn  a 


m 


the  square  of  the  disturbing  force  being  neglected. 

By  means  of  this  relation,  when  one  of  the  long  inequali- 
ties is  known,  the  other  may  be  calculated  :  it  may  be  used  as 
a  formula  of  verification. 

95.  We  have  remarked  that  in  integrating  the  equations 
of  Art.  90,  we  may  in  general  consider  the  elements  which 
enter  in  the  right-hand  members  as  constant  and  equal  to 
their  values  at  the  epoch  from  which  the  time  is  measured. 
In  the  case,  however,  of  inequalities  whose  periods  are  very 
long,  the  secular  variations  of  the  elements  in  the  interval 
produce  a  sensible  effect.  In  order  to  take  account  of  these, 
we  may  integrate  our  equations  by  parts,  considering  the 
elements  variable  ;  and  then  substitute  their  values  as  calcu- 
lated by  the  method  of  Art.  62.  For  example,  consider  the 
equation 


=  Psin  X,  suppose. 

Integrating  by  parts,  and  remembering  that  n  is  constant 
with  regard  to  secular  variations,  we  have 


I  • 
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P  1          dP   .    . 

—7  COS  X  +  7 7V2  -JT  sm  * 

pn—qn  (pn  —  qn)    dt 

I          d*P 


(pn-qnf  dt2 

i       *  P  t  •-'•;•«» 

therefore  o?  =  —  7 ^2  sin  X  —  ^ r^  -,- 

(^>n  -  qn  )2  (^>^  —  ^ )    dt 

+ 


-T-  cos  X  + 


(pn-qrif  dt^         T(pn-qn')*  df 


1         <Z2P  .   _ 

smX  + ... 


(pn  - 
dP 


d*P         I    . 

4  df  "  '  '  '  J  S 


df 

smX 


- 

(pn-qn)*dt 

In  this  equation  P,  -7-,   -77-,  &c.  are  functions  of  the 
dt       dt 

elements  ;  their  values  may  be  calculated  by  the  formulae  of 

77} 

Art.  62.     It  may  be  noticed  that  P  is  of  the  first  order,    -  of 


the  second,  and  —^  of  the  third  of  the  disturbing  force: 

dt 
-j-p 

for  -j-  ,  being  found  from  P  by  differentiation,  will  involve 
ctt 

the  differential  coefficients  of  the  elements,  which  are  them- 

dzP 

selves  of  the  first  order  ;  and  similarly  for  -^  . 

96.    Having  now  completed  our  account  of  the  methods 
c.  P.  T.  7 


•  I 
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of  treating  the  secular  and  periodic  variations  of  the  elements 
of  the  orbit,  we  will  say  a  few  words  011  the  distinction  be- 
tween them.  In  the  first  place  we  may  observe  that  the 
periodic  variations  involve  the  mean  longitude  of  the  ^dis- 
turbed  and  disturbing  planets,  and  therefore  depend  chiefly 
upon  the  configuration  of  the  planetary  system.  On  the  con- 
trary the  secular  variations  depend  solely  upon  the  values  of 
the  elements.  The  latter  class  of  variations  take  place  with 
extreme  slowness,  so  that  if  these  only  existed,  a  considerable 
time  must  elapse  before  the  deviation  of  the  planet  from 
elliptic  motion  became  appreciable.  On  the  other  hand,  the 
periodic  variations  (such  at  least  as  are  rapidly  periodic) 
"  are  in  their  nature  transient  and  temporary :  they  disappear 
in  short  periods,  and  leave  no  trace.  The  planet  is  tempo- 
rarily drawn  from  its  orbit  (its  slowly  varying  orbit),  but 
forthwith  returns  to  it,  to  deviate  presently  as  much  the  other 
way,  while  the  varied  orbit  accommodates  and  adjusts  itself 
to  the  average  of  these  excursions  on  either  side  of  it ;  and 
thus  continues  to  present,  for  a  succession  of  indefinite  ages, 
a  kind  of  medium  picture  of  all  that  the  planet  has  been 
doing  in  their  lapse,  in  which  the  expression  and  character 
is  preserved;  but  the  individual  features  are  merged  and 
lost*."  On  this  account  it  is  convenient  to  suppose  the 
planet  to  move  in  an  ellipse,  the  elements  of  which  are  cor- 
rected for  secular  variations  only,  and  to  take  account  of  the 
periodic  variations  by  applying  small  corrections  to  the  radius 
vector  and  longitude  as  calculated  from  the  elliptic  formulae. 

} 

97.     We  will  accordingly  shew  how  by  means  of  the 

periodic  variations  of  the  elements,  the  corresponding  in- 
equalities in  the  radius  vector  and  longitude  may  be  calcu- 
lated. If  we  take  for  our  plane  of  reference  the  position  of 

*  Herschel's  Outlines  of  Astronomy,  5th  edit.  Art.  656. 
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the  plane  of  the  orbit  of  the  disturbed  planet  at  the  epoch 
from  which  the  time  is  reckoned,  the  inclination  will  be  of 
the  order  of  the  disturbing  force,  and  therefore  if  we  neglect 
the  square  of  the  latter,  we  may  also  neglect  the  square  of 
the  former. 

98.  To  calculate  the  periodic  variations  in  radius  vector. 

Let  Sa,  Se,  SOT,  &c.  denote  the  periodic  variations  in 
a}  e,  OT,  &c.,  and  let  oY  be  the  corresponding  variation  in  r  ; 
then 

j,       dr  ~        dr  .        dr  .         dr  ^     dr  . 
or  =  -j-  ba  +  -j-  be  +  -j—  OOT  +  -77,  of  +  -y  oe, 
da          Je          ^OT  c?f          cfe 

in  which  the  square  of  the  disturbing  force  is  neglected, 
since  this  would  be  introduced  by  the  squares  and  products 
of  Sa,  &e,  &c.  The  values  of  8a,  Se,  &c.  have  been  found  in 

Art.  91,  those  of  -y-,    -=-,  &c.  may  be  obtained  from  the 
da     de 

equation  (Art.  40) 

r  =  a  jl+-e2-ecos  (f+e  —  -or)  —  -e2cos  2 

\ 

99.  To  calculate  the  periodic  variations  in  longitude. 

These  might  be  found  in  the  same  manner  as  the  varia- 
tions in  radius  vector,  but  they  may  also  be  deduced  from 
them  :  we  proceed  to  obtain  a  formula  for  this  purpose. 

We  have  =     >     (Art.  22), 


and 

- 


therefore  —  =     +  (1  -  cos  i)  ~     (see  Art.  29) 


7—2 
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since  (1  —  cos  i)  —T-  being  of  the  order  of  the  square  of  the 

disturbing  force  may  be  neglected. 

Let  Sr,  SO,  and  Sh  be  corresponding  variations  in  r,  0 

and  h  ;  then 

d( 

dt 
JO     dW      h  /       Bh 

or     :77  +  -^r==- 

dt       dt       r 


neglecting  the  square  of  the  disturbing  force  ;  therefore 


dt   ~  r*        r*    > 

which  gives  the  variations  in  longitude.     The  value  of  Bh 
may  be  found  from  the  formula 

dh  =  dE     dR 

dt      de      dix  ' 

For  the  periodic  variations  in  latitude,  we  refer  to  Ponte- 
coulant's  Systime  du  Monde,  Tome  I.  p.  492. 

100.     As  an  example  of  the  processes  of  this  chapter,  we 
will  calculate  the  variations  in  radius  vector  and  longitude 
due  to  the  term  m'Me  cos  {(n  —  2ri)  t  +  e  —  2e'  +  -cr}  in  JR. 
Considering  this  term  only,  we  have 

E  =  m'Me  cos  {(n  -  2ri)  t  +  e  -  2e'  +  ^} 
=  m'Me  cos  X,  suppose. 

dR        ,  dM  dR 

Hence      —r=m  ——  e  cos  X,        -j~  =m  Jf  cos  X, 

da  da  de 

dR  ,   _    .   .          dR 

-7-=  —  wiJfesinX,        -j—  =  —  m  Me  sm  X, 

dR  dR 
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Substituting  these  in  the  formulae  of  Art.  39,  and  neg- 
lecting small  quantities  of  orders  higher  than  the  first,  we 
have 

da         2na*    tlt,    . 
~-Ti~ m  Me  sin  X, 

de     na     t^r  . 
3J-  — wJfBinX, 

dm     na    ,,, 
e  -T-  =  —  m  .M  cos  X, 

de         2nd*     ,dM  I  na      ,,, 

-77  = m  -j-ecos\  +  - — mJfcosX, 

at  p         da  2  fj, 

—rr  = m  Me  sin  X. 

of       p 

By  integration  we  have 

.,       2m  M    na*e 

oa  = — ,  cos  X, 


m'M     na 

----  r-7  cos  X, 
p     n  —  2n 

m'M      na 


n  m'M     2m'a  dM\     nae      . 
=  I  ----  -;  —    -  r-/  sin  X, 
\2    fju  fjudajn-2n 


*,,        3m  M  .    . 

6f  =  ---  7  -  —  TTj  smX. 
p      (n  —  2n)z 


Also 


r  =  a    l  +    e2  -  e  cos  (f+  e  -  -or)  -    e2  cos  2  (f+  e  -  «r)  - 
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therefore,  small  quantities  of  orders  higher  than  the  first  being 
neglected, 

aS-l-«coB(?+e-,r),     -    • 

dr 

—  =  a  [e  —  cos  (f  +  e  —  «r)  —  e  cos  2  (f  +  e  -  «r)}, 

-r-  =  —  a(esin(f  +€  — 

WOT 


-T-  =oesm 

XT        «,       Jr  o,       dr  .,        <7r  -,         dr  ^^     dr  ~ 

Now   Sr  =  -j-  3a  +  -j-  Se  +  -j-  SCT  +  -^  Sf  +  -7-  Se 


2m'  M    na?e 


7  cos 


ii        -w        9<w' 
LL         71  —  Alb 

m'M     na*  .  ,„  N  , 
;— 7  cos  X  {e  —  cos  (f  +  e  —  -57)  —  e  cos  2  (f  +  e  — 


-7  sinX  {sin  (f  +  e  -  «r)  +  e  sin  2  (f  +  e  -  «r)J. 


Since  we  are  neglecting  the  square  of  the  disturbing  force, 
the  elements  in  this  equation  may  be  considered  as  constant, 
and  therefore  nt  written  for  f  :  we  have  then,  restoring  to  X 
its  value 

'  (f         n  ,  /) 

cos  2  \(n  —  n  )  f  +  e  —  e  } 


n- 


m'M    ncfe  ,,  , 

— ,  cos  {(n  -  2»)  « +  e  -  2e 

yu.     n  —  2n 
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r/0        ^  /\  .  rt  t        i 

cos  l3?i  —  2n    t  -f  3e  —  2e  -r  «r, 


m'M    ncf 


which  is  the  variation  in  radius  vector. 

101.     To  calculate   the  variation  in  longitude,  we  shall 
employ  the  equation 


Now          =       +        =  - 


.,       -          ^T 

therefore     SA  =    _^  ,  cos  \  ; 


therefore     -,  =   8 ,       a  '\  eos  * 

-^-,  cos  {(n  -  2n')  « -f  e  —  2c' 

Also  — j-  =  — j-  {1  • 


5  m  M    n* 


Hence  by  substitution 
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Sm'M    n*ae 


_  , 

—  2e 


0  , 

—  2e  — 
n  —  2n 

By  integration 


/A     (w  —  2n  )  (n  —  w 
m'M      n*ae       .    (, 

sm   »- 


5m  M 


yu,      (n  —  2»)(3ii—  .2») 
which  is  the  variation  in  longitude. 

In  the  case  of  Uranus  and  Neptune,  since  n  :  n  nearly 
as  2  :  1,  the  term  we  have  been  considering  is  important  in 
the  theory  of  the  long  inequality. 


CHAPTER  VII. 

DIEECT  METHOD  OF  CALCULATING  THE  INEQUALITIES  IN 
RADIUS  VECTOR,   LONGITUDE,  AND  LATITUDE. 

102.  IN  the  calculation  of  the  planetary  inequalities,  we 
have  hitherto  employed  exclusively  the  method  of  the  Varia- 
tion of  Elements,  but  there  is  another  method  of  solving  the 
problem,  which  demands  our  attention.    It  consists  in  obtain- 
ing equations  for  calculating  the  inequalities  in  radius  vector, 
longitude,  and  latitude  directly  from  the  equations  of  motion. 
The  two  methods  are  sometimes  distinguished,  the  former  as 
that  of  Lagrange,  the  latter  as  that  of  Laplace.     In  practice 
both  are  employed,  that  of  Lagrange  chiefly  for  secular,  that 
of  Laplace  for  periodic  variations:  but  the  method  of  La- 
grange  may  also  be  advantageously  employed  in  the  calcula- 
tion of  long  inequalities.     We  proceed  then  to  explain  the 
method  of  Laplace. 

103.  If  rl9  01?  and  z  denote  the  projected  radius  vector, 
longitude,  and  distance  from  the  plane  of  reference,  of  the 
planet,  we  have  (see  Art.  9)  the  equations  of  motion 


d_f  2  dO\  _  dR 
dt\l   dt)~~d§S 


d*z         iz     dR 
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If  we  take  for  the  fixed  plane  of  reference  the  position  of 
the  plane  of  the  orbit  at  the  epoch  from  which  the  time  is 
measured,  the  inclination  (as  we  have  remarked  in  Art.  97) 
will  be  the  order  of  the  disturbing  force,  the  square  of  which 
will  be  neglected.  Now  it  will  be  seen  on  referring  to  Art. 
42,  that  rt  and  0t  differ  from  r  and  0  by  quantities  depending 
upon  the  square  of  the  inclination  :  hence  in  the  above  equa- 
tions, we  may  replace  rl  and  0t  by  r  and  6  respectively.  Also 
if  X  denote  the  latitude  of  the  planet,  we  have 

z  =  r  sin  X. 
Hence  our  equations  of  motion  become 


104.  As  a  first  approximation,  let  values  of  r,  6  and  \ 
be  obtained  from  these  equations  by  neglecting  the  disturbing 
force,  and  let  r  +  Sr,  0  +  S0,  X  +  SX  denote  the  true  values  of 
these  co-ordinates;  then  Sr,  &0  and  SX  will  be  very  small 
quantities,  of  the  order  of  the  disturbing  force:    they  are 
termed  the  perturbations  in   radius  vector,  longitude,   and 
latitude.     We  proceed  to  investigate  equations  by  means  of 
which  these  quantities  may  be  determined. 

105.  To  obtain  the  equation  for  the  perturbation  m  radius 
vector. 

From  equations  (1)  and  (2)  of  Art.  103,  we  obtain 
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Multiply  (1)  by  r  and  add  it  to  (4)  :  thus 


If  the  disturbing  force  be  neglected,  this  equation  be- 
comes 


Let  a  value  of  r  be  obtained  from  this  equation,  and  let 
r  +  Sr  denote  the  true  radius  vector  :  then  if  we  agree  to 
neglect  the  square  of  the  disturbing  force  in  our  next  ap- 
proximation, it  will  be  sufficient  to  write  r  -f  Br  for  r  in  those 
terms  of  (5)  which  do  not  involve  the  disturbing  force  :  also 
since  Br  is  itself  of  the  order  of  the  disturbing  force,  its  square 
may  be  neglected.  Hence  from  (5) 


-      +»+4 

dt  r  +  cr  dr        J     at 

,,       ,         <Z2(r2)      nd*(rSr}      ty     2  ft  ,  dR 

therefore   --A^  +  2     ^     ;  =_^-J^3r  +  gr 

ar  £?<2  r        r  or 


hence  by  (6), 
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which  is  the  equation  for  the  perturbation  in  radius  vector. 
We  may  express  the  right-hand  member  in  a  more  convenient 
form,  for  since 

r  =  a  (1  +  u), 


da      dr  da  dr  ' 

.,       P  dR        dR 

therefore  r  —j-  —  a-^-i 

dr         da 

,  d(R]        dR 

also  — \ — -  =  n  — . 

Hence  our  equation  becomes 

d*(r8r)      *,  dR  ,     _FdB 


106.     To  obtain  the  equation  for  the  perturbation  in  longi- 
tude. 


We  have  from  equation  (1)  of  Art.  103, 

Y  =  i—     £-1^-? 
,dt)      r  df      r3     r  dr  ' 

As  before,  let  a  value  of  0  be  obtained  from  this  equation, 
the  disturbing  force  being  neglected,  and  let  6  +  80  denote 
the  true  value  of  the  longitude:  then  writing  r  +  Sr  for  r 
and  6  +  80  for  6}  and  neglecting  the  square  of  the  disturbing 
force,  we  have 


_.  ___ 

dt~dt~r   d?  "7ZW~^ 

i  Q  dO     ,         dR        dR 

but  rz-j-  =  h,     r-j-  =  a-j-; 

dt  dr         da 
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f          n,  d^r      3   d\      3/4  s  dR 

therefore     2h  -7-  =  r  -7-0  —  or  -=-5-  --  £  or  —  a  -7- 
dt          dr  dt       r  da 

d  /    dSr      ^   dr\      Su,  dR 

=  :*»   r  ~^7  —  &•  -y-   --  T  •  ro?*  —  a-j-  . 
dt\    dt  dt)      r  da 

This  equation  will  become  integrable  if  we  eliminate  the 
term  —  ~  .  rSr  by  means  of  the  equation  for  the  perturbation 
in  radius  vector.  We  have  from  that  equation 

,    3/*         5  ,      d 

-L  +  -Jj.rSr-3a-; 
r3  da 


A          O  ,  5  ,  fJR  J* 

0  =  3  —      -L  +  -J.rSr-3a-;  --  Qn    -j-dt; 

J  de 


therefore  by  addition, 

_ ,  dBO  _     d2  (roV)      d  (^  dSr      ^  dr\      A  _  dR 
~dt~        d?          c 


^       f        «7o/i      n  d&r      ^  dr  CdR 

therefore  2^8^  =  3     v,    }  +  r-^  -Sr-j-  -±a  l-^ 
dt  dt  dt          J  da 


dt  dt          J  da  JJ  de 

0 
or        ^o^  =  2 


, 
dt 


j,   dr  [dR  ,  FfdRj, 

—t)r-r  —  2a\^j-dt  —  3nH-7-dr, 
dt          J  da  JJ  de 


which  determines  the  perturbation  in  longitude,  when  that  in 
radius  vector  is  known. 

107.     To    obtain    the    equation  for    the  perturbation    in 
latitude. 

From  equation  (3)  of  Art.  103,  we  have 
d*  (r  sin  X)      p  (r  sin  \)  __  dR 
~~~     ~"' 
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Since  the  plane  of  reference  is  supposed  to  coincide  with 
the  position  of  the  plane  of  the  orbit  at  the  epoch  from  which 
the  time  is  reckoned,  we  have  at  the  epoch  X  =  0  :  hence, 
denoting  by  SX  the  latitude  at  time  t,  our  equation  becomes 


dz  ' 

which  is  similar  in  form  to  the  equation  for  the  perturbation 
in  radius  vector. 

108.     To  integrate  the  equation  for  the  perturbation  in 
radius  vector. 

The  equation  is  (Art.  105), 

dz(rSr)       a  ,  dB  ,     ..  [  dB 


,  fi  ,  j,  .         dR  [dB  , 

+  5  (ror)  =  a-j-  +  2n  I  -;-  dt. 

r*  x      '         da          J  de 


df 

Let  us  consider  a  term  in  B  of  the  form 
Pcos{(jpn-2»')  t+  Q], 

where  P  is  a  function  of  the  mean  distances,  excentricities, 
and  inclinations,  and  Q  of  the  longitudes  of  the  perihelia, 
nodes,  and  epochs:  then  uniting  this  term  with  the  non- 
periodic  part  of  B,  which  we  have  denoted  by  F,  we  have 


=  F+  Pcos  {(pn  -  gn]  t+Q}; 

dP 
-& 


c        dB     dF    dP       ,.  „  ,     ^ 

therefore    -      =       +  -     cos  {(pn  -  gin  )  t  +  Q}, 


since  F  does  not  contain  e  ;  therefore 
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where  g  is  an  arbitrary  constant,  the  value  of  which  may  be 
any  whatever.  We  might  omit  it,  and  consider  it  as  in- 
cluded in  the  constant  C  of  Art.  105,  but  the  advantage  of 
retaining  it  will  be  seen  hereafter. 

dB  [dR  ,,          ,          dF 

Hence          t**- 


n 

dP  (  *nr  de  \  cos  {(pn  -  qri]  t  +  Q} 

[    da     pn—gii] 


f 
'g  +  a-+Plcos  {(pn-qri)t  +  Q}9 


suppose,  where        P. —  a  -j-  +  , . 

da     pn  —  qn 

Again  (see  Art.  13), 
r  =  a  \  I  +  -  e*  —  e  cos  (nt  +  e  —  «r) 

e2  cos  2  (nt  +  e  —  -or)  —  . . .  !• , 

2  ) 

therefore         -s  =  —*  f  1  • 


snce 


Hence,  by  substitution,  the  equation  for  the  perturbation 
in  radius  vector  becomes 

-{.  Q} 


-  7i2  r$r  {3e  cos  (nt+  e  -  w)  +  ...}. 

109.     This    equation   inust  be  solved  by  successive  ap- 
proximation,  as    in  the  Lunar  Theory.     By   omitting  all 
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small  quantities,  we  obtain  a  first  approximation  to  tlie  value 
of  rSr;  this  being  substituted  in  the  second  member,  and 
small  quantities  of  orders  higher  than  the  first  neglected,  we 
obtain  a  second  approximation,  which  will  be  correct  to  the 
first  order.  In  like  manner,  a  third,  and  higher  approxima- 
tions may  be  obtained. 

On  referring  to  Art.  59,  it  will  be  seen  that  small  quanti- 
ties of  the  second  order  being  neglected, 


hence,  neglecting  all  small  quantities,  the  equation  of  the 
preceding  Article  becomes 


The  integral  of  this  equation  is 
m'     dC 


where  A  and  B  are  arbitrary  constants.  Since,  if  all  small 
quantities  be  neglected,  r  =  a,  we  have  as  a  first  approxima- 
tion, 

1    /      ,       m'     dC0\ 
Sr  =  ^(2^+T^) 

T) 

ir  i^ncos{(^n-Sfn/)*+C} 


a  [nz  —  (pn  —  qn ) 

^ 

+  —  cos  (nt  —  B}> 
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110.  We  may,  however,  omit  the  last  term:  for,  con- 
sidering this  only,  the  radius  vector  of  the  planet  becomes 

a  \  1  —  e  cos  (nt  +  e  —  tzr)  +  ^  cos  (nt  —  B)  +  ...[ 

A 

=  a  [1  -  {e  cos  (e  —  w) §  cos  B}  cos  w£ 

A 

+  {e  sin  (e  —  -or)  +  —z  sin  B}  sin  n£  +  . . .] 

—  a  [1  —  e1cos>  (nt-\-  e  —  wj  +...}, 
if  ex  cos  (e  —  -cj-J  =  e  cos  (e  —  -cr) ^  cos  B, 

€^  sin  (e  —  -crj  =  e  sin  (e  —  ur)  +  -2  sin  i?, 

from  which  e1  and  ^  may  be  determined. 

Now  since  the  ellipse  upon  which  our  approximations  are 
based,  has  been  obtained  by  neglecting  the  disturbing  force, 
we  may  in  the  elliptic  formulae  replace  e  and  CT  by  e^  and  ^ 
respectively,  since  they  differ  by  quantities  of  the  order  of  the 
disturbing  force.  If  this  be  done,  our  first  approximation 
becomes 

Sr  =  4-(2M>  +  £a^) 

n'a  \  2daJ 

p 

H — r~2 — r-1 ^77  cos  {(Pn  —  qn)  I  +  Q}- 

a  \n  —  (pn  —  qn)  \ 

111.  In   order  to  obtain  a  second  approximation,   this 
value  must  be  substituted  for  &r  in  the  right-hand  member  of 
the  equation  of  Art.  108.     Also  since  the  square  of  the  dis- 
turbing force  is  neglected,  we  may  write  e^  and  «rx  for  e  and  ix 
in  this  equation.    We  will  write  for  brevity 

Sr  =  L  +  P2 cos  {(pn  -  qri]  t+Q}. 
C.  P.  T.  8 
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Substituting  this  in  the  equation  of  Art.  108,  and  omitting 
those  terms  which  have  produced  the  first  approximation*, 
we  have 


=  -  3tt\  cos  (nt  +  6-vr1)[L  +  P2  cos  {(pn  -  qri)  t  +  Q}] 
=  — •  3n^e1  L  cos  (nt  +  e  —  -GrJ 

W^PgCOS 

3 

112.  On  the  form  of  this  equation,  we  have  an  important 
remark  to  make.  In  consequence  of  the  term 

—  3n^e:  L  cos  (nt  +  e  —  -crj, 
its  integral  will  contain  the  term 

-  e^ntL  cm^nt  +  e  —  wj. 

Here,  then,  we  are  met  by  a  difficulty:  our  equations 
Irave  been  formed  on  the  hypothesis  that  the  square  of  Sr  is 
small  enough  to  be  omitted,  whereas  here,  we  have  a  term 
capable  of  indefinite  increase.  This  term  then,  if  retained, 
would  ultimately  vitiate  the  whole  approximation.  The  diffi- 
culty might,  as  in  the  Lunar  Theory,  be  obviated  by  writing 
en  for  n  in  the  elliptic  formula,  which  amounts  to  supposing 
the  perihelion  to  be  in  motion.  Its  motion  is  however  better 
found  by  the  method  of  the  variation  of  elements.  Indeed  it 

*  These  terms  are  omitted  for  the  sake  of  brevity:  in  order,  therefore,  to 
obtain  the  complete  second  approximation,  we  must  add  to  the  integral  of  the 
above  equation  the  result  of  the  first  approximation. 
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may  be  shewn  that  such  terms  as  those  we  are  considering, 
lead  to  the  formulae  which  have  already  been  obtained,  for  the 
secular  variation  of  the  elements*.  We  shall  accordingly 
altogether  neglect  such  terms,  and  suppose  the  elements  of  the 
ellipse  on  which  our  approximations  are  based,  to  have  been 
previously  corrected  for  their  secular  variations. 

113.     With  this  understanding,  the  complete  integral  of 
the  equation  of  Art.  Ill,  will  be 

3  *£»  T> 

rBr  =  —  - 


2  »*-{(/?  +  !)  n-qn'}* 


2  n*  -  {(p  -l)n-  qri} 

cos  [{(p-l)n-  qri}  t  +  Q  -  e  +  «rj 
•f  A  cos  (nt  —  B). 

If  this  be  added  to  the  result  of  the  first  approximation, 
we  obtain  for  a  second  approximation 

1   /      ,       m     JOA  P. 

rSr=-i(2mg  +  -a 


cos  {(pn  —  qn)  t  +  Q] 


cos  [{(p  +  l)n-  qn 


cos  [{(p  -  1)  n  -  qri}  t+Q-e  +  vr^ 
+  A  cos  (nt  —  B). 

*  It  is  thus  that  the  Secular  Variations  are  first  obtained  in  the  Mecanique 
Celeste.     See  Ponte'coulant's  Systdme  du  Monde,  Supplement  au  Livre  II. 

8—2 
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The  arbitrary  constants  might  Ibe  disposed  of  as  in  the 
first  approximation,  but  it  is  more  convenient  in  practice  to 
determine  them  otherwise.  In  order  to  obtain  a  second  ap- 
proximation to  the  value  of  Sr,  it  is  only  necessary  to  multiply 
the  right-hand  member  of  the  above  equation  by 

-  {1  +  et  cos  (nt  -f  e  —  vrj}, 
neglecting  e*. 

114.     To  calculate  the  perturbations  in  longitude. 

We  have  (Art.  106) 

z*a      nd.r$r      jj    dr  fdR  ,  f  [dR  _. 

JiW  =  2  —  =  --  Br  -y  —  2a    -7-  dt  —  %n        y-  dr. 
dt  dt          J  da  JJ  de 

Taking  for  simplicity,  the  first  approximation  to  the  value 
of  rSr,  which  has  been  obtained  by  neglecting  the  first  power 
of  the  excentricity,  we  have 

1   /  m      dC\  P. 

rSr  =  -2    277i'^  +  —  a-7-0   +^  —  -.  —  s  --  7^ 
n  \  2       da  /      n  2  —  (p/i  —  gw  ) 

cos  {(pn  —  gn]  t  +  Q]  ; 
therefore 

d.rSr          2Pn—n       . 


also,  neglecting  the  first  power  of  the  excentricity, 


=    . 
dt 

By  Art.  108,  writing  -  m'C0  for  F,  we  have 

j£ 

dR     1    ,  cZC0    rfP 
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dQ 


nP 

PMMQOi 

pn  —  gn 
therefore  —  2a  I  ~-  ^  — 


and        n\--r-dt=  -  ^r 


3nP 


pn  —  qn  y      da     pn  —  gn 

where/  is  an  arbitrary  constant. 
Hence  by  substitution,  we  have 

=/-  (m'a  ^  +  3m' 


(jw  —  <p'      ( jw  —  gn')'z 

4-    a  t   , ^^r  sm  {(Pn  ~  2W')  ^+  C}- 

»  —  (^?w  —  gn} ) 

115.  This  expression  is  open  to  the  objection  of  contain- 
ing a  term  proportional  to  the  time,  which  being  capable  of 
indefinite  increase,  would  ultimately  vitiate  the  whole  approxi- 
mation. Here,  then,  we  see  the  advantage  of  having  a 
quantity  g  which  may  be  determined  at  pleasure :  we  will  so 
determine  it  that  the  objectionable  term  shall  vanish.  This 
condition  gives 


We  may  also  omit  the  constant  /,  and  consider  it  as  con- 
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tained  in  the  epoch.     We  have  then,  writing  for  h  its  value 
-  neglecting  e2, 

dp 


s. 


116.  Before  proceeding  to  obtain  the  perturbations  in 
latitude,  we  will  make  a  few  remarks  on  the  forms  of  the  ex- 
pressions for  Sr  and  80.  If  we  confine  ourselves  to  the  results 
of  the  first  approximation,  it  will  be  seen  on  substituting  the 
value  of  Pj,  that  pn  —  qn  and  n*  —  (pn  —  qn')2  occur  as  divi- 
sors, and  that  the  expression  for  80  contains  besides,  the 
divisor  (pn  —  qn')z.  The  second  of  these  may  be  written 


If  then  either 
(i)  pn  -  gn',     (ii)   (1  -p)  n  +  qn  ',     or  (iii)   (1  -jp)  n  -  qri, 

be  very  small,  the  corresponding  terms  in  &r  and  80,  though  of 
a  high  order,  may  yet  be  sensible.  This  is  especially  the  case 
with  the  first,  since  as  we  have  remarked,  its  square  occurs 
in  the  expression  for  80.  These  are  instances  of  what  in 
the  preceding  chapter  have  been  characterised  as  long  in- 
equalities. 

The  period  of  the  term  Pcos  {(pn  —  qn)  t+  Q],  which  has 
given  rise  to  these  inequalities  is 

27T 

pn  —  qn  ' 

in  the  case  of  (i),  this  is  very  large,  and  in  that  of  (ii)  or  (iii) 

f)—, 
it  is  very  nearly  equal  to  —  ,  since  pn  —  qn'  is  nearly  equal 
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to  ±  n.  Hence  it  appears  that  terms  in  E  whose  period  is 
either  very  large,  or  nearly  equal  to  that  of  the  planet,  may 
give  rise  to  important  inequalities  in  the  radius  vector  and 
longitude.  Their  actual  importance  will  of  course  depend  in 
part  upon  the  order  of  the  principal  part  of  Pt  with  respect  to 
the  excentricities  and  inclinations,  i.e.  (see  Art.  50)  upon 


117.     To  integrate  the  equation  for    the  perturbation   in 
latitude. 

The  equation  is  (Art.  107) 

d*(rS\]      »  dR 


the  position  of  the  plane  of  the  orbit  of  the  disturbed  planet  at 
the  epoch,  being  taken  for  the  fixed  plane  of  reference. 

Differentiating  the  expression  for  R  in  Art.  44,  with  respect 
to  z,  we  obtain 


-m  (z-z')  (JZ>0  +  D,  cos  </>+  ...  +  A  cos 

or,  putting  z  equal  to  0,  and  substituting 

a'  tan  H  sin  (n't  +  e  —  Of) 
for  z  (see  Art.  42), 

-  ma  tan  i'  sin  (n't  +  e'  -  Q')  j-  D0  +  -,3  +  A  cos  </>  +  ...i  . 

This  expression,  after  reduction,  consists  of  terms  of  the 

form 

Q], 


where  p  and  q  are  positive  integers,  and  either  may  be  zero. 
Considering  one  such  term,  our  equation  becomes 

—  -  +    rSx  =  Psin 
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Now  as  in  Art.  108, 

-3  =  ft2  (1  +  3e  cos  (nt  +  e  —  «r)  +  ...}; 
hence,  neglecting  the  product  eS\, 


The  integral  of  this  equation  is 
p 


If  instead  of  taking  for  the  fixed  plane  of  reference,  the 
plane  of  the  orbit  of  the  disturbed  planet  at  epoch,  we  take  a 
plane  slightly  inclined  to  this,  we  may  omit  the  arbitrary 
term.  For,  denoting  the  planet's  latitude  with  respect  to  this 
plane  by  X,  we  have  approximately 

X  =  tan  i  sin  (nt  +  e  —  11), 

and  it  may  be  shewn  as  in  Art.  110,  that  omitting  the  term 
in  question  is  only  equivalent  to  changing  slightly  the  values 
of  i  and  H. 


CHAPTER,  VIII. 

ON    THE   EFFECTS  WHICH   A   RESISTING    MEDIUM   WOULD 
PRODUCE   IN  THE   MOTIONS  OF   THE   PLANETS. 

118.  IN  the  preceding  chapters,  we  have  supposed  the 
planetary  motions  to  take  place  in  free  space,  and  the  results 
of  calculations  based  upon  this  hypothesis  manifest  a  very 
close  agreement  with  observation.  There  is,  however,  a 
remarkable  circumstance  connected  with  Encke's  comet  which 
seems  to  indicate  the  possibility  of  the  existence  of  a  very  rare 
medium,  too  rare  indeed  to  cause  any  sensible  resistance  to 
the  motions  of  the  planets,  but  which,  as  we  shall  presently 
see,  may  yet  influence  the  motions  of  comets,  in  consequence 
of  the  extreme  smallness  of  the  masses  of  these  bodies.  It 
has  been  observed  that  the  comet  above  referred  to  (which 
describes  an  elliptic  orbit  in  a  period  of  about  3j  years,)  has 
since  its  appearance  in  1786,  been  moving  round  the  Sun 
with  an  increasing  mean  motion.  Encke  attributes  this  to 
the  resistance  of  a  medium  pervading  space.  We  shall 
therefore  proceed  to  examine  the  effects  which  such  a  medium 
would  produce  upon  the  elements  of  a  planet's  orbit,  assuming 
the  resistance  to  vary  as  the  product  of  the  density  of  the 
medium  and  the  square  of  the  velocity  of  the  planet.  We 
shall  neglect,  in  the  present  investigation,  all  forces  except 
the  Sun's  attraction  and  the  resistance  of  the  medium ;  conse- 
quently the  planet  may  be  supposed  to  move  wholly  in  one 
plane. 
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119.  Let  r,  6  be  the  radius  vector  and  longitude  of  the 
planet,  s  the  length  of  an  arc  of  its  actual  orbit  measured  from 
some  fixed  point  to  its  position  at  time  £,  and  p  the  density  of 
the  medium.  Then  if  Jc  be  a  constant,  we  may  represent  the 

resistance  on  the  planet  by  kp  (-J-  J  ,  and  the  equations  of  mo- 
tion will  be 

d'r  dO\z         ju  dszdr 


7/1 

If  r2  -j-  =  hj  these  may  be  written 

.d'r          de\z          L      ,    dsdr 


d  (  2  d&\         7    7  ds 

—  (r*-j-)  =-7cph-r  .......................  (2). 

dt\    dt)  r    dt 


These   equations   are   the  .  same   in  form  with   those  of 
Art.  20,  and  may  be  treated  in  a  similar  manner,  —kp-^-  -j 

taking  the  place  of  -7-  ,  and  —  Jcph  -j-  that  of  -^  .    We  have 

from  equation  (2) 

dh  ..  ds 


120.  To  obtain  a  formula  for  calculating  the  mean  dis- 
tance. 

We  might  proceed  as  in  Art.  25,  but  we  shall  here  employ 
the  method  of  Art.  26.  We  have 

d*s          i  dr  ds 
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and  by  a  known  formula  of  elliptic  motion 


dt  r       a' 

Differentiating  the  latter,  we  obtain 

ds  d*s  _     2yt6  dr      fi  da 

2diW~~fdi  +  a*~di 
ds 

and  multiplying  the  former  by  2  j-  , 

dsd2s         2/4<fr     07    fdss 

*VV*-7*~*r\ 

.,       .  p  da  .    fds\* 

therefore  %  Tt  =  ~  2k/>  (di)  > 

-da 

= 


121.     To  obtain  a  formula  for  calculating  the  excentricity. 
We  have,  as  in  Art.  27, 


Differentiating  as  if  »•  were  constant,  and  -writing  —k--r 


,         r 

for    •-• 


=-kpdtdt)+7 
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Now  from  equation  (3) 


2)  /g       \  ds 

(r~l)dt' 


This  result  may  also  be  obtained  by  differentiating  the 
formula  h?  =  fj,a(l  —  e2),  and  substituting  the  expressions  for 

dh       ,  da 

•j-  and  -y-  ,  as  in  Art.  28. 

122.     To  obtain  a  formula  for  calculating  the  longitude  of 
perihelion. 

We  have,  as  in  Art.  29, 


Differentiating  as  if  r  and  9  were  constant,  and  writing 

7    dr  ds  r      d*r 

—  kp  -T.  -T.  tor  -T-O  , 


dr          2  .  -         .  J-nr      j    dr  ds        ,fi         ^      fl  t  p  \  dh 
dh  ,  ds 


or  since  7t  =  ~        dt  > 
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dr  v  d'uf      ,    a 


but  from  equation  (4)  of  Art.  22, 


dr  fQ         .      fM 

_cosec(0—  )-=?; 


therefore  -7-  = "  sm  (0  —  «")  T  • 

dt  e  J  dt 

123.  Before  proceeding  to  obtain  a  formula  for  calcu- 
lating the  epoch,  we  shall  express  the  results  of  the  preceding 
articles  in  a  form  convenient  for  application. 

If  u  denote  the  excentric  anomaly,  we  have 

r  —  OL  (1  —  e  cos  u) (1), 

nt  +  e  —  tar  =  M  —  e  sin  u (2). 

Hence  ( -y 

\dt 


a  \\  —  e  cos  u 

_  fjb  I  +  e  cos  u  $ 
a  1  —  ecosw' 

6?s  //I  +  e  cos  tA       .  a 

therefore       -7-  =  wa  A  /   -    ,  since  no?  =  v  w. 

a^  V  V1  ~-  e  cos  w/ 

From  (1),  by  differentiation  as  if  the  elements  were  inva- 
riable, we  have 

dr  du 
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and  from  the  equation 


dr 


dr 
equating  the  two  values  of  -j-  , 

(jit 

.     //i         v       ha    .       du 


From  (2),  by  differentiation  as  if  the  elements  were  inva- 
riable, we  have 

dt      \  ..  . 


-j-  =  -  U  —  &  cos 

du     n  ^ 


Hence,   by  substitution,   the   formulas  of  the  preceding 
articles  become 


da  ,     2  , 

-j-  =  -  2&pa2  (1  +  e  cos  u 

du 


N      //I  4-  e  cos  u\ 
)  A  /  (  -  -    , 

'  y   \l  -  e  cos  u)  ' 

//1  +  ecostA 
A  /  -  -    , 
V  Vl-ecosw/' 


-j     —        Arfft/fsiA/  v  •*•  ™"  v  J  COS  cv 

du 

dix  _     2/£jo« \/(l  —  e2)    .  //I  +  c  cos  iA 

rtw  e  \/  \1  —  e  cos  M/  " 

124.     To  obtain  a  formula  for  calculating  the  epoch. 
By  differentiating  the  formulae 

r  =  a  (1  —  ecos  w), 
—  -sr  =  u  —  e  sin  w, 


considering  the  elements,  first  variable,  and  then  invariable, 
and  comparing  the  results,  we  obtain 

du  de      .  .  da 

ae  sin  u  -=-  —  a  cos  u  -, —  (1  —  e  cos  w)  -y- , 


ndt  +  e.     Thus,  on  substituting  the  values  of  -^ ,   ~  ,   and 
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dn     de     d^  du  de 

t  -j-  +  -T-  — -j-  =  (1  —  e  cos  u)  -j sin  u  -r  : 

dt      dt      dt      v  '  dt  dt' 

whence,  eliminating  -y- 

efo,     o?e     c?-cr  __  cos  u  —  ede      (\  —  e  cos  w)2  da 
dt      dt      dt        e  sin  u     dt          ae  sin  u       dt  ' 

fl 77 

As  in  Art.  37,  we  may  omit  the  term  t  -7-  if  we  bear  in 

mind   that  the   mean    longitude   will   then  be   denoted  by 

de 
dt 

—j- ,  and  reducing,  we  obtain 

de      %Jcpa  ,„  *\        a          •»    •  //I  + 

-y-  =  — —  jl  —  VU  -  e)  —  e  cos  Mj  sin  w  A  /  (  -  ,  . 

du        e  V  U  ~ e  cos  w/ 

125.  The  formulae  of  the  preceding  articles  are  sufficient 
to  determine  the  elements  of  the  orbit  at  any  time,  and  being 
perfectly  general,  are  applicable  as  well  to  the  motion  of 
comets,  as  to  that  of  planets,  but  before  we  can  integrate 
them,  we  shall  require  a  knowledge  of  the  form  of  p.  Now 
the  analogy  of  the  terrestrial  atmosphere  would  lead  us  to 
suppose  that  if  the  sun  be  surrounded  by  an  ethereal  medium, 
its  density  decreases  as  the  distance  from  the  sun  increases. 
Moreover,  the  researches  of  Professor  Encke  on  the  comet 
which  bears  his  name,  seem  to  indicate  the  law  of  the  in- 
verse square.  We  will,  however,  assume  p  to  be  such  a 

function  of  r.  that  when  multiplied  by  A/(r — - —  -)  ,  and 

J  V  V1  ~~ e  cos  UJ 

developed  in  a  series  of  cosines  of  u  and  its  multiples,  it  takes 
the  form 

A  +  Be  cos  u  +  Cez  cos  2w  +  ... 
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Thus  our  formulas  become 

da 

-:-  =  —  2Jca?  {A  +  (A  +  B)  e  cos  u  +...}, 

de                  (  .               Be  1 

-7-  =  —  2&a  \A  cos  u  H (1  +  cos  2w)  +  ...  k 

dm         n7     (  A    •          Be   .  1 

e  ~T~  —  —  2&a  -j  -4  sin  w  H sm  2u  +  ,...>•  , 

>•  / 

-T-  =  ka  (Ae  sin  w  +  ...). 
aw 

126.  Supposing  the  orbit  nearly  circular,  to  examine  the 
effect  of  the  medium  upon  the  elements  of  the  orbit. 

Since  the  orbit  is  nearly  circular,  we  shall  neglect  squares 
and  higher  powers  of  e ;  thus  the  preceding  formulas  give 
on  integration 

a  =  const.  —  2&a2  [Au  +  (A  +  B)  e  sin  u}, 

«7     f  A  Be  f        sin  2 

e  =  const.  -  2ka  ^Aewu  +  —(u+  — 

{Be 
A  cos  u  H 
4 

e  =  const.  —  lea  Ae  cos  u. 

Hence  in  an  entire  revolution  of  the  planet,  the  mean 
distance  is  diminished  by  ±TrlcazA,  and  the  excentricity  by 
%7rkaBe,  while  the  longitudes  of  perihelion,  and  of  the  epoch 

remain  unchanged.  Also  from  the  formula  n  =  ~- ,  it  ap- 
pears that  the  mean  motion  is,  in  an  entire  revolution,  in- 
creased by 
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127.  We  have  already  remarked  that  no  traces  of  a 
resisting  medium  have  yet  been  discovered  in  the  motion  of 
the  planets :  but,  since  k  varies  inversely  as  the  mass  of  the 
body  acted  upon,  the  formulas  of  Art.  123  shew  that  such 
a  medium,  though  too  rare  to  influence  the  planets,  might  yet 
sensibly  affect  the  motions  of  comets,  in  consequence  of  the 
extreme  smallness  of  their  masses. 


C   P.  T. 


PROBLEMS. 


1.  SUPPOSING  in  the  Problem  of  the  Three  Bodies  the 
relative  orbit  of  two  of  the  bodies  to  be  a  circle  described 
uniformly,  obtain  equations  for  determining  the  motion  of  the 
third  body ;  and  transform  the  system  of  co-ordinates,  so  that 
the  plane  of  the  circular  orbit  being  that  of  xy,  the  axis  of  x 
shall  always  pass  through  the  two  bodies  in  that  plane. 

2.  Shew  that  the  plane  of  the  orbit  of  a  planet  revolves 
about  the  planet's  radius  vector  as  an  instantaneous  axis*. 

3.  A  particle  is  describing  an  orbit  round  a  centre  of 
force  which  is  any  function  of  the  distance,  and  is  acted  upon 
by  a  disturbing  force  which  is  always  perpendicular  to  the 
plane  of  the  instantaneous  orbit,  and  inversely  proportional 
to  the  distance  of  the  body  from  the  centre  of  the  principal 
force.     Prove  that  the  plane  of  the  instantaneous  orbit  re- 
volves uniformly  round  its  instantaneous  axis. 

4.  Find  when  the  curvature  of  the  instantaneous  orbit  of 
a  body,  acted  on  by  disturbing  forces,  is  the  same  as  that 
of  the  actual  orbit ;    and  shew  that  this  is  always  the  case 
when  the  only  disturbing  force  arises  from  the  action  of  a 
resisting  medium. 

*  In  this  and  the  following  problem,  the  plane  of  the  orbit  must  be  supposed 
to  have  no  angular  velocity  about  a  normal  to  itself.     See  note  to  Art.  19. 
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5.  If  R  be  expressed  on  the  one  hand  as  a  function  of 
rx,  0l9  and  z  (Art.  11),  and  on  the  other  as  a  function  of  r, 
0,  t,  and  £1,  0  being  measured  on  the  plane  of  the  orbit  from 
the  node,  prove  that 


and  obtain  a  formula  for  calculating  the  inclination. 

6.     If  R  be  expressed  as  a  function  of  t  and  the  usual 
elements,  obtain  the  formulae 


dt      h  sin  i  di  ' 

di_cot(0-n)dR 

dt~         h      ~  di  ' 

where  0  is  measured  on  the  plane  of  reference  as  far  as  the 
node,  and  thence  on  that  of  the  orbit,  and 


h  =  r2  [ -7-  -  2  sin2  -  — =-  )  . 
\dt  2  dtj 


i 
7.     The  central  force  being  -a  +  ~  ,  obtain  the  following 

equation  for  the  apsidal  motion 


d™  _  Vq(l  -ea)  /A'  cos  (0  - 
dt  ~ 


a,  e  and  &  being  elements  of  the  instantaneous  ellipse. 

8.  A  body  revolving  about  a  centre  of  force,  which  varies 
inversely  as  the  square  of  the  distance,  is  constantly  retarded 
by  a  small  constant  force;  find  the  alteration  of  the  major 
axis,  excentricity,  and  apse,  in  one  revolution. 

9—2 
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9.  When  the  disturbing  function  E  is  independent  of  6, 
find  expressions  for  -=-  and  -7-. 

m' 

If  ^  =  —  ,  these  expressions  give  variable  values  for  e 

and  -cr,  whereas  the  motion  of  the  body  actually  takes  place 
in  a  fixed  ellipse:  shew  this,  and  explain  the  apparent 
paradox. 

10.  A  planet  describes  an  orbit  under  the  action  of  a 
force  ^  tending  to  the    Sun,  p  not  being  quite   constant: 

obtain  the  following  equations  for  the  variations  of  the  ex- 
centricity and  longitude  of  perihelion ; 

=  —  cos  (6  —  OT), 


'CT  .       fn  . 

lie  -j—  =  —  sin  (u  —  CT). 

If  djjb  be  always  positive,  what  in  a  whole  revolution  is 
the  nature  of  its  effects  upon  the  excentricity  and  position  of 
the  major  axis  ? 

11.  If  the  equation  of  the  Moon's  orbit  be  reduced  to 
the  form 


shew  that  the  excentricity  and  longitude  of  perihelion  may 
be  found  from  the  equations 

de  /.  •     //i         \  dnr       /.        //i         \ 

^=  -/sin  (0-*r),     e  —  =/cos  ((9  -  «r). 

Apply  these  equations  to  find  e  and  -or,  when  /  is  a  small 
disturbing  force,  depending  only  upon  the  Moon's  distance 
from  the  Earth. 
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12.    Assuming  the  differential  equation  for  s  in  the  Lunar 
Theory  to  be 

d*s  -    (3     3        rt  x/1 

—  +  s  =  -  w25   -  +  -  cos  2  (0  - 


shew  that  if  7  be  the  longitude  of  the  Moon's  node, 

^  =  -  |  m2  {1  -  COS  2  (010  -  7)  -  COS  2  ((9-7) 

+  cos2(0-ro0)J. 

From  the  above  expression  for  -Jr  ,  find  the  ratio  of  the 

mean  motion  of  the  node  to  that  of  the  Moon,  taking  into 
account  terms  of  the  order  m4. 

13.  If  two  planets  disturbing  one  another  were  revolving 
in  periods  of  350  and  201  days,  what  form  of  terms  in  the 
disturbing  function  would  demand  examination  ? 

14.  The  periods  of  Venus  and  the  Earth  are  224*7  and 
365*256  days  respectively;  find  approximately  the  period  of 
the  long  inequality  arising  from  their  mutual  perturbations, 
the  important  term  in  the  disturbing  function  R  being  of  the 

form 

PeV2  cos  (13  (nt  +  e)  -  8  (rit  +  e')  -  3*7  -  2w'}. 

15.  The  radius  vector  of  a  planet  is  affected  with  a  small 
periodical  inequality  ;  shew  that  its  effect  may  be  represented 
by  continued  and  periodical  alterations  of  the  excentricity  and 

PT 

longitude  of  perihelion,  the  period  of  either  being   ^^  ^  , 

where  P  is  the  period  of  the  planet  and  T  that  of  the  in- 
equality. 
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16.  If  in  addition  to  the  force  of  the  Sun  on  a  planet 
there  be  a  small  force  tending  towards  the  Sun,  and  varying 
inversely   as   the  wth  power  of  the  distance  of  the   planet 
from  the  Sun,  prove  that  the  perihelion  of  the  orbit  will  have 
a  progressive  or  regressive  motion  according  as  m  is  greater 
or  less  than  2. 

Can  you  explain  this  result  by  reasoning  similar  to  that 
used  in  Airy's  Gravitation  ? 

17.  It  has  been  found  by  comparing  theory  with  ob- 
servation that  the  perihelion  of  Mercury  progresses  at  a  rate 
greater  by  a  than  that  due  to  the  attraction  of  known  bodies  : 

shew  that  this  increment  would  be  accounted  for  if  the  law  of 

/  i 

force  tending  to  the  Sun  were  ~  +  ^  ,  and  if  tf  =  ac4  A  /-  ,  the 

T  T  y'      C 

orbit  being  supposed  to  be  nearly  a  circle,  and  the  mean  dis- 
tance to  be  c. 

18.  The  central  force  acting  on  a  body  being 


shew  to  terms  inclusive  of  p'  and  the  square  of  the  excen- 
tricity,  that  the  motion  is  in  an  ellipse  revolving  uniformly 
about  the  focus. 

19.     Shew  by  means  of  the  formula 

da  _  2na*  dR 
dt        IJL     de 

that  the  chief  perturbation  of  the  axis  major  of  the  Moon's 
orbit  may  be  expressed  by  the  equation 


=  a    1  +  ^~T  -  ~\  cos  2  (nt  +  €  ~  n't  ~ 
2n  (n  —  n  ) 
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where  n  and  ri  are  the  mean  motions  of  the  Moon  and  Sun 
respectively. 

20.  A  satellite  revolving  in  an  ellipse  of  small  excen- 
tricity  is  disturbed  by  another  satellite  revolving  about  the 
same  primary  ;  find  approximately  the  variation  of  the  mean 
distance  and  the  motion  of  the  apse,  corresponding  to  the 
terms 

^  r*  [1  +  3  cos  J2  (n  -  ri)  t  +  e  -  e}] 

in  the  function  R,  having  given 


na 


-  e2  dR 


_ 
dt        fi,    de  '     dt  fjLe          de  * 

21.  Prove  that,  neglecting  periodical  variations,  the  ex- 
centricity  of  any  orbit   can   always  be  represented  by  the 
diagonal  of  a  parallelogram,  whose  sides  are  constant,  and 
angle  varies  uniformly. 

22.  Given  the  equations 

tan2  1  =  A?  +  N*  +  2N,NS  cos  (hj  +  Bl  -  S2), 


fan  o  _ 
" 


explain  the  nature  of  the  motion  of  the  node,  when  the  mini- 
mum inclination  is  zero. 

23.  Prove  that  as  far  as  secular  variations  only  are  con- 
cerned the  function  F  is  constant. 

24.  Considering  only  secular  variations,  obtain  the  fol- 
lowing equations  : 


/«      ^N  2/™tan.^  = 

\na      dt)  \na  dt  J 
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25.  If  the  squares  of  the  masses  of  two  mutually  dis- 
turbing planets  were  to  each  other  inversely  as  their  mean 
distances,  shew  that  the  nodes  would  oscillate  through  equal 
angles. 

26.  If  Mj  m,  m  be  the  masses  of  three  bodies  mutually 
attracting  according  to  the  law  of  gravity,  M  being-  much 
larger  than  m  or  m',  and  if  v,  v  be  the  velocities  of  m,  m  at 
distances  r,  r'  from  the  centre  of  Jf,  supposed  fixed,  shew  that 
the  equation  of  vis  viva  for  this  case  may  be  assumed  to  be 

2        t  '2     c*n/rfm      m      m>      m'\     * 
mv*  +  mv*  +  2if  -  ---  +  —  ,  --  ;    =  0, 
\2a      r      2a       r  J 

%a  and  2a  being  the  major  axes  of  the  instantaneous  ellipses 
of  m  and  m. 

27.  Infer  from  the  foregoing  equation  by  the  method  of 
the  variation  of  parameters  the  ratio  of  simultaneous  changes 
in  the  mean  distances  and  mean  motions  of  two  planets  mutu- 
ally disturbing. 

28.  If  r  be  the  true  radius  vector,  6^  the  projected  longi- 
tude, and  X  the  latitude  of  a  planet,  obtain  the  following 
equation  of  motion  : 

d\\z      i     dR 


29.  Obtain  the  following  equation  between  the  pertur- 
bations of  a  planet  in  longitude  and  radius  vector,  whatever 
be  the  law  of  force,  provided  it  be  central  and  a  function  of 
the  distance  only,  and  provided  such  a  function  as  E  can  be 
found  : 


dt  j     dt  dr 

-4^Sr-2r 

dr 
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where  F  denotes  the  central  force,  and  Ji  twice  the  sectorial 
area  described  by  the  undisturbed  planet  round  the  Sun. 

30.  If  the  orbits  of  two  planets  which  disturb  each  other 
be  very  nearly  circular,   shew  that  the  inequalities  of  the 
radius  vector  may  be  immediately  deduced  from  those  of  the 
longitude  by  means  of  the  equation 

d.W 

Tt 

31.  Integrate  the  equation 

-4r^  +  n*  •  r$r  =  2  {Pcos  (put  +  Q)}, 
determining  the   arbitrary  constants    so    that    Sr  =  0,    and 

7       r\ 

-4 —  =  0,  when  t  =  0 :  and  shew  that  for  small  values  of  tt 

Ptz 


the  case  of  p  =  1  being  included. 

32.     A  planet  moves  in  a  resisting  medium  of  which  the 
resistance 


apply  the  equation 


to  obtain  the  following,  in  which  e2  is  neglected  : 


df 


+  n*.rSr  +  n*.  rSr .  3e  cos  (nt  +  e  —  -or) 

+  2/^2a  [nt-+  —  e  sin  (nt  +  €  -  «r)}  =  0. 
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33.  The  co-ordinates  of  the  position  at  any  time  t  of  a 
disturbed  planet  being  x  +  &e,  y  +  §y,  z  +  $z,  reckoned  from 
the  Sun's  centre  as  a  fixed  origin,  and  referred  to  the  plane 
of  motion  at  a  given  epoch;  and  r  being  the  heliocentric 
distance,  x,  y  the  co-ordinates  of  the  position  which  the 
planet  would  have  had  at  the  time  t,  if  the  disturbance  had 
ceased  at  the  given  epoch  ;  obtain  the  following  equations  for 
determining  Sx,  Sy,  8z  to  the  first  order  of  the  disturbing 
force : 

^(«»)  ,/ 


in  which  JJL  is  the  sum  of  the  masses  of  the  Sun  and  planet, 
and  R'  is  put  for 


<m  being  the  mass,  and  a/,  y  ,  z  the  heliocentric  co-ordinates 
of  the  disturbing  planet. 

34.  Shew  that  the  effect  of  a  resisting  medium  on  the 
instantaneous  orbit  of  a  planet,  would  be  to  make  the  apsidal 
line  regrede  or  progrede,  according  as  the  planet  moved  from 
perihelion  to  aphelion,  or  from  aphelion  to  perihelion. 

35.  Two  small  planets  P,  $,  very  near  each  other,  re- 
volve about  the  Sun  in  orbits  very  nearly  circular,  and  make 
two  revolutions  about  each  other  while  they  make  one  revo- 
lution about  the  Sun.     Compare  the  sum  of  their  masses  with 
the  mass  of  the  Sun. 
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•  •    If  the  line  PQ  move  parallel  to  itself,  what  inference  do 
you  draw  ? 

36.  If  the  motion  of  a  planet  round  the  Sun  be  disturbed 
by  the  action  of  another  planet,  the  latter  being  supposed  to 
describe  a  circular  orbit  of  radius  a  with  uniform  velocity  n, 
obtain  the  following  exact  equation : 


2m' 

cos  a)  H 


where  r  is  the  radius  vector  of  the  disturbed  planet,  rt,  01?  z 
its  co-ordinates  referred  to  a  fixed  plane,  and  co  the  inclina- 
tion of  the  radii  vectores  of  the  disturbed  and  disturbing 
planets  to  each  other. 

37.  Prove  that,  if  the  periodic  times  of  a  disturbed  and 
disturbing  planet  are  not  commensurable,  the  secular  changes 
of  the  orbit  of  the  disturbed  planet  are  the  same  as  they 
would  be  if  the  mass  of  the  disturbing  planet  were  distributed 
over  its  orbit,  in  such  a  manner,  that  the  part  of  the  mass 
distributed  over  each  portion  of  the  orbit  should  be  propor- 
tional to  the  time  which  the  planet  actually  takes  to  describe 
that  portion. 


APPENDIX. 

ON  THE  FORM  OF  THE  EQUATIONS  OF  ART.  39. 

1.  ON  referring   to  Art.  39,  it  will   be  seen   that  the 
formulae  which  have  been  obtained  for  calculating  the  elements 
of  the  orbit  involve  only  partial  differential  coefficients  of  R 
with   respect   to  these   elements,    multiplied   by  coefficients 
which  do  not  contain  the  time  explicitly.     This  circumstance 
greatly  facilitates  their  application,  by  rendering  them  fit  for 
use  as  soon  as  the  partial  differential  coefficients  have  been 
calculated.     We  proceed  to  shew  that  this  remarkable  cha- 
racteristic is  not  restricted  to  the  particular  system  of  elements 
which  have  been  adopted,  but  may  be  attained  with   any 
system  of  elements  whatever.     The  discovery  of  this  is  due 
to  Lagrange. 

2.  If  the  motion  of  the  planet  be  referred  to  three  rect- 
angular axes  originating  in  the  centre  of  gravity  of  the  Sun, 
we  have  the  equations  of  motion  (see  Art.  9) 

d*x     fix     dR  ... 

~rW+^=^r (    '' 

ut  T  ttX 

cfy    ^_dE 

df  +  r1  ~  dy ()> 


f*     u*     dR 

W  +  ^=d^ (3)> 
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Let  a,  b,  c,  d,  e,  f  be  the  six  elements  introduced  by  in- 
tegrating these  equations  when  R  —  0,  and  for  -v- ,   -M- ,  -T- 

write  a?',  y',  z:  then  a;',  ?/',  and  z  can  be  expressed  as  functions 
of  t  and  the  elements ;  hence 

dx'  da     dx   db 


where  in  f-T-J  the  elements  are  supposed  constant. 

If  in  equation  (1)  we  put  R  equal  to  0,  we  have 
<dx^ 


,       f  x      (x\ 

therefore  --= —    —r- }  =  -j- ; 

dt      \  dt  J      dx' 

dR     dx   da     dx   db 
therefore  _=-_  +  _-  + 


and  similar  equations  hold  for  -, —  and  -, — . 

Now  since  R  is  a  function  of  x,  y,  and  £, 

dR  _  dR  dx     dR   dy     dR  dz^ 
da  ~  dx  da      dy    da      dz    da 

_  /dx  dx'      dy  dy_      dz_  dz  \  da 
"~  \da  ^a      da  ~da      da   da  )  dt 

/dx  dx'      dy  dy       dz_  dz'\  db 
\da   db       da   db       da   db )  dt 

(dx  dx       dy  dy       dz_  dz  \  dc 
(da'dc'^'da'cic'da  lie)  dt 
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3.     We  may  eliminate  -—  from  this  expression :  for,  sup- 
at 

posing  x,  y,  and  z  expressed  as  functions  of  t  and  the  elements, 
we  have 

dx  _  (dx\      dx  da      dx  db 

~di  ~  (dt)  +  da  dt  +  db  di 

but  by  the   principles  of  the  method  of  the  Variation   of 
Parameters 

dx     fdx\ 

dt  ~(dtj  ' 

,       »  dx  da     dx  db     dx  dc 

therefore          -7-  -y  +-77  -j-  +-. — r  +  ...=0. 
da  at      ao  at      dc  dt 

a.    .,    ,       dy  da     du  db     dy  dc 
Similarly,    •£  -rr-K-*  ^r  +  -^  -77+  ...  =  0, 


-  - 

da  dt      db  dt        c  dt 

dz  da      dz   db      dz  dc 


Multiplying  these  equations  by  -?-  ,  -jj-  ,  &c.5  and  adding, 

we  obtain 

_  (dx  dx      dy  dy      dz  dz'\  da 
\da  da      da  da      da  da)  dt 

(dx  dx      dy  dy      dz  dz\  db 
*(dbda*db  da  +db~da)  Jt 

(dx  dx      dy  dy       dz  dz'\  dc 
\dc  da      dc  da      dc  da)  dt 
+  ............... 

If  this  expression  be  subtracted  from  that  for  -7—  in  Art.  2. 
the  latter  may  be  written 

dR      r     7n  db      r      n  dc 
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r     7n      dx  dx      dx  dx      dy  dy      dy  dy 

where  [a,  0]  =  -j-  -jr  —  -JT  -j — \-  -j    -jrr  —  -JT  ~r~ 

L       J      da  db      db  da      da  db      db  da 

dz  dz      dz  dz 

Similarly, 

dR      ri     -i  da     n     -.  dc  . 


4.  By  successive  elimination  between  these   equations, 

we   can   obtain  expressions  for    -7-,   -7-,  &c.,    in  terms  of 

at      fit 

/77?      /77? 

-7-  ,   -37- ,  &c.,  [a,  5],  [a,  c],  &c. :  if,  then,  we  can  shew  that 

[a,  J],  [a,  c],  &c.,  are  independent  of  the  time  explicitly,  it  will 

follow  that  this  is  also  the  case  with  the  coefficients  of  -y— , 

da 

dR     p        .     ,,  .        r     da     db     p 

-77- ,  &c.,  in  the  expressions  for  -r  j   T  »  &c- 
ao  at      at 

5.  To   shew  that  [a,  b]   t*   independent  of  the  time  ex- 
plicitly. 

Let  F  =  -  ;  then  the  equations  of  motion  give 


dtj       dx 

Now  differentiating  with  respect  to  t  only  so  far  as  it  occurs 
explicitly, 

dr-     ,-,      dx  d  fdx\  .  dx'  d 
dt  da  o 


_dx_d^  (<M\__dx_d_  (fa\ 
~dbdt(da)     ~daJt  (db) 
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_  dx  d  (dx\      dx^  dx^ 

~~  ~~i      Ti    I       T7  I     '        / 1    ~7 


dbda 


dx  d  fdx\  _  dx  dx' 
~"db  da  \dt)  ~~~da  ~db 


V 
+ 


=  __ 

~~dadb\dx)      db  da\dx) 


__dy_  <    (_  \ 
da  db  \  dy        db  da  \  dy  ) 


dz   d  (dV\      dz  d  (dV\ 
da  db  \dz  )      db  da\dz) 

dx  d  (dV\      dy  d  (dV\     ^z_  d_  (^\ 
=  dadx  \db)  ^tady  \dbl  +  da  dz  \db  ) 

dx  d  fdV\      dy  d  fdV\      dz  d  (dV\ 
~~dbdx  \da)  ~~aldy  \da)  ~Tb~dz  \da) 


da  db     db  da 

Hence  [«,  Z>]  does  not  contain  the  time  explicitly.  The 
same  is  of  course  true  of  [a,  c],  [b,  c],  &c.  It  follows,  then, 
that  whatever  system  of  elements  be  adopted,  we  can  always 
express  their  differential  coefficients  in  terms  of  the  partial  dif- 
ferential coefficients  of  E  with  respect  to  them,  multiplied  by 
coefficients  which  do  not  involve  the  time  explicitly. 

6.  From  the  formula  of  Art.  3  of  this  Appendix,  which 
is  due  to  Lagrange,  those  of  Chapter  II.  may  be  deduced :  for 
this  we  refer  to  Pontecoulant's  Systeme  du  Monde,  Tom.  I. 
p.  542. 
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ON    THE    GEOMETRICAL    INTERPRETATION    OF  THE   FORMULA 

FOR  THE  SECULAR  VARIATIONS   OF  THE  NODE 

AND   INCLINATION. 

7.  In  Art.  84  we  have  shewn  that  in  consequence  of  the 
secular  variations  of  the  node  and  inclination,  the  normal  to 
the  plane  of  the  orbit  moves  uniformly,  so  as  to  generate  in 
space  a  fixed  right  circular  cone.  This  result  may  also  be 
obtained  somewhat  differently:  we  will  here  indicate  the 
process  *. 

The  equations  to  be  interpreted  are 

p  =  tan  i  sin  ft  =  N^  sin  (hjt  +  8X)  +  JV2  sin  S2 , 
q  =  tan  i  cos  ft  =  ^  cos  (hjk  +  SJ  +  N2  cos  S2. 

Employing  the  figure  and  construction  of  Art.  84,  it  will 
be  found  that 

cot  QA  0  =  tan  *  sin  ft,   cot  QL  0  =  tan  i  cos  ft. 

These  are  the  geometric  representations  of  p  and  q.  Now 
it  may  be  shewn  by  Spherical  Trigonometry  that 

n  A  n     tan  /sin  o>  +  tan  p  (sin  co  cos  6  +  sec  /cos  o>  sin  0} 

COt  ty  A.  U  = 7 T~- -7; — — , 

1  —  tan  1  tan  p  cos  u 

QJ  0  _  tan  /cos  6)  +  tan  p  (cos  co  cos  6  —  sec  /sin  <a  sin  0) 
1  —  tan /tan  p  cos  6 

These  expressions  are  rigorous ;  if  we  neglect  the  product 
tan2/  tan  p,  and  higher  products  of  tan/  and  tan  p,  they 
become 

cot  QA  0  =  tan  /sin  G>  +  tan  p  sin  (6  +  o>), 

cot  QLO  —  tan  /cos  &>  +  tan  p  cos  (0  +  «). 

*  This  method  as  well  as  that  of  Art.  84,  is  due  to  Mr  Freeman,  of  St  John's 
College. 

C.  P.  T.  10 
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From  these  equations  the  interpretation  follows  as  in  Art. 
84.  The  advantage  of  this  method  is  that  it  affords  a  geo- 
metric representation  of  p  and  q:  on  the  other  hand,  in 
the  method  of  Art.  84,  the  trigonometrical  reduction  is 
simpler. 

ON  THE  METHODS  OF  CALCULATING  THE   MASSES  OF 
THE  PLANETS. 

8.  There  are  in  general  two  methods  of  determining  the 
masses  of  the  planets ;  either  by  observations  on  a  satellite, 
when  the  planet  is  accompanied  by  a  satellite  ;  or  by  compar- 
ing the  inequalities  produced  in  their  motion  by  their  mutual 
action,  as   deduced   from   observation,   with   the   same   ine- 
qualities calculated  from  theory.     The  secular  variations  are 
best  adapted  to  give  the  most  exact  results ;  but  these  are  not 
yet  known  with  sufficient  accuracy  to  allow  of  this  use.     We 
are  therefore  obliged  to  recur  to  the  periodic  variations,  and, 
by  combining  a  vast  number  of  observations,  gather  from 
them  the  most  probable  results*. 

9.  When  the  planet  is  accompanied  by  a  satellite  the 
formula  for  calculating  its  mass  may  be  obtained  as  follows  : 

Let  M,  m,  m  be  the  masses  of  the  Sun,  the  planet,  and  the 
satellite :  P,  P'  the  periodic  times  of  the  planet  about  the  Sun, 
and  the  satellite  about  the  planet ;  a,  a  the  mean  distances  of 
the  planet  from  the  Sun,  and  the  satellite  from  the  planet. 
Then  we  have 

o__4 

r.    p'=- 


m  +  mf      Pzan 
therefore 


*  Ponte"coulant's  Systime  du  Monde,  Tome  m.  p.  340. 
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or  approximately, 


This  equation  gives  the  mass  of  the  planet  when  that  of 
its  satellite  is  known.  If  the  latter  be  neglected,  the  formula 
becomes 

m      P*a* 


10.  In  the  case  of  the  Earth,  this  method  is  not  suf- 
ficiently exact,  but  the  following  may  be  employed.  The 
attraction  of  the  Earth  on  a  body  at  its  surface,  in  the  parallel 

of  which  the  square  of  the  sine  of  the  latitude  is  -  ,  is  very 

o 

nearly  the  same  as  if  the  Earth  were  condensed  into  its 
centre.  (See  Pratt's  Figure  of  the  Earth,  Art.  89.)  Let  then 

sin2  1  —  -  ,  g  =  the  Earth's  attraction  on  a  body  at  its  surface  in 

o 

latitude  I,  b  the  mean  radius  of  the  Earth,  E  the  mass  of  the 
Earth,  M  the  mass  of  the  Sun,  P  the  length  of  the  year,  and 
a  the  mean  radius  of  the  Earth's  orbit.  Then 

_E  2™* 

-  = 


E      gVP*      gP*  /b\5 
-     =  =   _    - 


/b\5 

therefore  -^  =  ,   ,  ,  =  ?_  (  -  )  , 

M      4wV     4?r26  \a)  ' 

where          -  =  sine  of  Sun's  parallax  =  sin  8".57. 


11.  For  the  methods  of  calculating  the  mass  of  the 
Moon  we  refer  to  Ponte'coulant's  Systeme  du  Monde,  Tome 
IV.  p.  651. 
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Tables  of  the  numerical  values  of  the  masses  of  the 
planets,  and  of  the  elements  of  their  orbits,  will  be  found  in 
Herschel's  Outlines  of  Astronomy,  pp.  693  et  seq. 


THE   END. 
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1.   Algebra. 

For  the  Use  of  Colleges  and  Schools. 

Third  Edition.      542  pp.  (1862). 
Strongly  bound  in  cloth.       7s.  6d. 

This  work  contains  all  the  propositions  which 
are  usually  included  in  elementary  treatises  on 
Algebra,  and  a  large  number  of  Examples  for 
Exercise.  The  author  has  sought  to  render  the 
work  easily  intelligible  to  students  without  im- 
pairing the  accuracy  of  the  demonstrations,  or 
contracting  the  limits  of  the  subject.  The  Ex- 
amples have  been  selected  with  a  view  to  illustrate 
every  part  of  the  subject,  and  as  the  number  of 
them  is  about  sixteen  hundred  and  fifty,  it  is 
hoped  they  will  supply  ample  exercise  for  the 
student.  Each  set  of  Examples  has  been  carefully 
arranged,  commencing  with  very  simple  exercises, 
and  proceeding  gradually  to  those  which  are  less 
obvious. 

2.   Plane  Trigonometry. 

For  Schools  and  Colleges. 

Second  Edition.    279  pp.  (1860),    Cm.  8vo. 
Strongly  bound  in  cloth.      5s. 

The  design  of  this  work  has  been  to  render  the 
subject  intelligible  to  beginners,  and  at  the  same 
time  to  afford  the  student  the  opportunity  of  ob- 
taining all  the  .information  which  he  will  require 
on  this  branch  of  Mathematics.  Each  chapter  is 
followed  by  a  set  of  Examples ;  those  which  are 
entitled  Miscellaneous  JSxamples,  together  with 
a  few  in  some  of  the  other  sets,  may  be  advan- 
tageously reserved  by  the  student  for  exercise 
after  he  has  made  some  progress  in  the  subject. 
As  the  Text  and  Examples  of  the  present  work  have 
been  tested  by  considerable  experience  in  teaching, 
the  hope  is  entertained  that  they  will  be  suitable  for 
imparting  a  sound  and  comprehensive  knowledge 
of  Plane  Trigonometry,  together  with  readiness  in 
the  application  of  this  knowledge  to  the  solution 
of  problems.  In  the  Second  Edition  the  hints  for 
the  solution  of  the  Examples  have  been  consider- 
ably increased. 


4.   The    Integral    Calculus 
and  its  Applications. 

With  numerous  Examples. 

Second  Edition.     342  pp.  (1862). 

Crown  8vo.  cloth.      10*.  Gd. 

In  writing  the  present  Treatise  on  the  Integral 
Calculus,  the  object  has  been  to  produce  a  work  at 
once  elementary  and  complete— adapted  for  the 
use  of  beginners,  and  sufficient  for  the  wants  of 
advanced  students.  In  the  selection  of  the  propo- 
sitions, and  in  the  mode  of  establishing  them,  the 
author  has  endeavoured  to  exhibit  fully  and  clearly 
the  principles  of  the  subject,  and  to  illustrate 
all  their  most  important  results.  In  order  that 
the  student  may  find  in  the  volume  all  that  he 
requires,  a  large  collection  of  Examples  for  ex- 
ercise has  been  appended  to  the  different  chapters. 


5.  The  Differential  Calculus 

With  numerous  Examples. 

Third  Edition,  398  pp.  (1860). 
Crown  8vo.  cloth,  10*.  6d. 

This  work  is  intended  to  exhibit  a  comprehen- 
sive view  of  the  Differential  Calculus  on  the 
method  of  Limits.  In  the  more  elementary  por- 
tions, explanations  have  been  given  in  considerable 
detail,  with  the  hope  that  a  reader  who  is  without 
the  assistance  of  a  tutor  may  be  enabled  to  acquire 
a  competent  acquaintance  with  the  subject.  More 
than  one  investigation  of  a  theorem  has  been 
frequently  given,  because  it  is  believed  that  the 
student  derives  advantage  from  viewing  the  same 
proposition  under  different  aspects,  and  that  in 
order  to  succeed  in  the  examinations  which  he 
may  have  to  undergo,  he  should  be  prepared  for 
a  considerable  variety  in  the  order  of  arranging 
the  several  branches  of  the  subject,  and  for 
a  corresponding  variety  in  the  mode  of  demon- 
stration. 


3,  Spherical  Trigonometry.     6. 

For  the  Use  of  Colleges  and  Schools. 

112  pp.  Crown  8vo.  (1859). 
Strongly  bound  in  cloth.      4s.  Qd. 

This  work  is  constructed  on  the  same  plan  as 
the  Treatise  on  Plane  Trigonometry,  to  which  it 
is  intended  as  a  sequel.  Considerable  labour  has 
been  expended  on  the  text  in  order  to  render  it 
comprehensive  and  accurate,  and  the  Examples, 
which  have  been  chiefly  selected  from  University 
and  College  Papers,  have  all  been  carefully  veri- 
fied. 


Analytical  Statics. 

With  numerous  Examples. 

Second  Edition.      330  pp.  (1858). 
Crown  8vo.  cloth.     10s.  6d. 


In  this  work  will  be  found  all  the  propositions 
which  usually  appear  in  treatises  on  Theoretical 
Statics.  To  the  different  chapters  Examples  are 
appended,  which  have  been  selected  principally 
from  the  University  and  College  Examination 
Papers ;  these  will  furnish  ample  exercise  in  the 
application  of  the  principles  of  the  subject. 
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7.  An  Elementary  Treatise 
on  the  Theory  of  Equa- 
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With  a  Collection  of  Examples. 
Crown  8vo.  cloth.      7s.  Gd. 

This  treatise  contains  all  the  propositions  which 
are  usually  included  in  elementary  treatises  on  the 
Theory  of  Equations,  together  with  a  collection  of 
Examples  for  exercise.  It  may  be  read  by  those 
who  are  familiar  with  Algebra,  since  no  higher 
knowledge  is  assumed,  except  in  Arts.  175,  267, 
308—314,  which  may  be  postponed  by  those  who 
are  not  acquainted  with  De  Moivre's  Theorem  in 
Trigonometry.  This  work  may  in  fact  be  regarded 
as  a  sequel  to  that  on  Algebra  by  the  same  writer, 
and  accordingly  the  student  has  occasionally  been 
referred  to  the  treatise  an  Algebra  for  preliminary 
information  on  some  topics  here  discussed. 

8.  Plane  Co-Ordinate  Geo- 

metry 

AS  APPLIED  TO  THE  STRAIGHT  LINE  AND 

THE  CONIC  SECTIONS ; 

With  numerous  Examples. 

Third  and  Cheaper  Edition. 
Crown  8vo.  cloth.      Is.  6d. 

This  Treatise  exhibits  the  subject  in  a  simple 
manner  for  the  benefit  of  beginners,  and  at  the 
same  time  includes  in  one  volume  all  that  students 
usually  require.  The  Examples  at  the  end  of 
each  chapter  will,  it  is  hoped,  furnish  sufficient 
exercise,  as  they  have  been  carefully  selected  with 
the  view  of  illustrating  the  most  important  points, 
and  have  been  tested  by  repeated  experience  with 
pupils.  In  consequence  of  the  demand  for  the 
work  proving  much  greater  than  had  been  origi- 
nally anticipated,  a  large  number  of  copies  of  the 
Third  Edition  has  been  printed,  and  a  consider- 
able reduction  effected  in  the  price. 


9.  Examples  of  Analytical 
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mensions. 

76  pp.  (1858).     Crown  8vo.  cloth.     4*. 

A  collection  of  examples  in  illustration  of  Ana- 
lytical Geometry  of  Three  Dimensions  has  long 
been  required  both  by  students  and  teachers,  and 
the  present  work  is  published  with  the  view  of 
supplying  the  want. 


10.  History  of  the  Progress 
of  the  Calculus  of  Vari- 
ations. 

During  the  Nineteenth  Century. 

It  is  of  importance  that  those  who  wish  to  culti- 
vate any  subject  jnay  be  able  to  ascertain  what 
results  have  already  been  obtained,  and  thus 
reserve  their  strength  for  difficulties  which  have 
not  yet  been  conquered.  And  those  who  merely 
desire  to  ascertain  the  present  state  of  a  subject 
without  any  purpose  of  original  investigation  will 
often  find  that  the  study  of  the  past  history  of  that 
subject  assists  them  materially  in  obtaining  a 
sound  and  extensive  knowledge  of  the  condition 
which  it  has  attained.  The  Author  has  en- 
deavoured in  this  work  to  ascertain  distinctly 
what  has  been  effected  in  the  Progress  of  the 
Calculus,  and  to  form  some  estimate  of  the  manner 
in  which  it  has  been  effected  :  accordingly,  unless 
the  contrary  is  distinctly  stated,  it  may  be  assumed 
that  any  treatise  or  memoir  relating  to  the  Cal- 
culus of  Variations  which  is  described  in  this  work 
has  undergone  thorough  examination  and  study. 
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Archdeacon  of  Calcutta,  late  Fellow  of  Gonville 
and  Caius  College,  Cambridge. 

A  Treatise  on  Attractions, 
!La  Place's  Functions,  and 
the  Figure  of  the  Earth. 

Second    Edition.       Crown    8vo.    126    pp. 
(1861).     cloth.     6s.  6d. 

In  the  present  Treatise  the  author  has  en- 
deavoured to  supply  the  want  of  a  work  on  a 
subject  of  great  importance  and  high  interest — 
La  Place's  Coefficients  and  Functions  and  the  cal- 
culation of  the  Figure  of  the  earth  by  means  of 
his  remarkable  analysis.  No  student  of  the  higher 
branches  of  Physical  Astronomy  should  be  ignorant 
of  La  Place's  analysis  and  its  result — "  a  calculus," 
says  Airy,  "the  most  singular  in  its  nature  and 
the  most  powerful  in  its  application  that  has  ever 
appeared." 


By    P.    G.    TAIT,    M.A.,    and 
W.   J.   STEELE,   B.A. 

Late  Fellows  of  St.  Peter's  College,  Cambridge. 

Dynamics  of  a  Particle. 

With  numerous  Examples. 
304  pp.  (1856).    Crown  8vo.  cloth.    10s.  6^. 

In  this  Treatise  will  be  found  all  the  ordinary 
propositions  connected  with  the  Dynamics  of 
Particles  which  can  be  conveniently  deduced  with- 
out the  use  of  D'Alembert's  Principles.  Through- 
out, the  book  will  be  found  a  number  of  illustrative 
Examples  introduced  in  the  text,  and  for  the  most 
part  completely  worked  out;  others,  with  occa- 
sional solutions  or  hints  to  assist  the  student  are 
appended  to  each  Chapter, 


By  REV.  8.  PARKINSON,  B.D. 

Fellow  and  Prselector  of  St.  John's  Coll.  Camb. 

1.  Elementary  Treatise  on 

Mechanics. 

With  a  Collection  of  Examples. 

Second  Edition.      345  pp.  (1861). 
Crown  8vo.  cloth.     9*.  6d. 

The  Author  has  endeavoured  to  render  the 
present  volume  suitable  as  a  Manual  for  the  junior 
classes  in  Universities  and  the  higher  classes  in 
Schools.  With  this  object  there  have  been  in- 
cluded in  it  those  portions  of  theoretical  Mechanics 
which  can  be  conveniently  investigated  without 
the  Differential  Calculus,  and  with  one  or  two 
short  exceptions  the  student  is  not  presumed  to 
require  a  knowledge  of  any  branches  of  Mathe- 
matics beyond  the  elements  of  Algebra,  Geometry, 
and  Trigonometry.  A  collection  of  Problems  and 
Examples  has  been  added,  chiefly  taken  from  the 
Senate-House  and  College  Examination  Papers — 
which  will,  it  is  trusted,  be  found  useful  as  an 
exercise  for  the  student.  In  the  Second  Edition 
several  additional  propositions  have  been  incorpo- 
rated in  the  work  for  the  purpose  of  rendering  it 
more  complete,  and  the  Collection  of  Examples 
and  Problems  has  been  largely  increased. 

2.  A   Treatise    on    Optics. 

304  pp.  (1859).     Crown  8,vo.     10s.  6d. 

The  present  work  may  be  regarded  as  a  new 
edition  of  the  Treatise  on  Optics,  by  the  Rev. 
W.  N.  Griffin,  which  being  some  time  ago  out  of 
print,  was  very  kindly  and  liberally  placed  at  the 
disposal  of  the  author.  The  author  has  freely 
used  the  liberty  accorded  to  him,  and  has  re- 
arranged the  matter  with  considerable  alterations 
and  additions — especially  in  those  parts  which  re- 
quired more  copious  explanation  and  illustration  to 
render  the  work  suitable  for  the  present  course  of 
reading  in  the  University.  A  collection  of  Ex- 
amples and  Problems  has  been  appended,  which 
are  sufficiently  numerous  and  varied  in  character 
to  afford  an  useful  exercise  for  the  student :  for 
the  greater  part  of  them  recourse  has  been  had  to 
the  Examination  Papers  set  in  the  University  and 
the  several  Colleges  during  the  last  twenty  years. 

Subjoined  to  the  copious  Table  of  Contents  the 
author  has  venturned  to  indicate  an  elementary 
course  of  reading  not  unsuitable  for  the  require- 
ments of  the  First  Three  Days  in  the  Cambridge 
Senate  House  Examinations. 


By  J.  B.  PHEAR,  M.A. 
Fellow  &  late  Mathematical  Lecturer  of  Clare  Coll. 

Elementary  Hydrostatics. 

With  numerous  Examples  and  Solutions. 

Second   Edition.     156  pp.  (1857). 
Crown  8vo.  cloth.      5s.  Qd. 

"An  excellent  Introductory  Book.  The  defi- 
nitions are  very  clear ;  the  descriptions  and  ex- 
planations are  sufficiently  full  and  intelligible; 
the  investigations  are  simple  and  scientific.  The 
examples  greatly  enhance  its  value." — ENGLISH 
JOURNAL  OF  EDUCATION. 

This  Edition  contains  147  Examples,  and  solu- 
tions to  all  these  examples  are  given  at  the  end  of 
the  book. 


By  G.  B.  AIRY,  M.A. 

Astronomer  Royal. 

1.    Mathematical  Tracts 

On  the  Lunar  and  Planetary  Theories,  Figure 
of  the  Earth,  the  Undulatory  Theory  of 
Optics,  $c. 

Fourth  Edition.   400  pp.  (1858).    8vo.    15*. 


2.  Theory  of  Errors  of  Ob- 
servations and  the  Com- 
bination of  Observations. 

103  pp.  (1861).     Crown  8vo.  6*.  6d. 

In  order  to  spare  astronomers  and  observers  in 
natural  philosophy  the  confusion  and  loss  of  time 
which  are  produced  by  referring  to  the  ordinary 
treatises  embracing  both  branches  of  Probabilities, 
the  author  has  thought  it  desirable  to  draw  up 
this  work,  relating  only  to  Errors  of  Observation, 
and  to  the  rules  derivable  from  the  consideration 
of  these  Errors,  for  the  Combination  of  the  Results 
of  Observations.  The  Author  has  thus  also  the 
advantage  of  entering  somewhat  more  fully  into 
several  points  of  interest  to  the  observer,  than  can 
possibly  be  done  in  a  General  Theory  of  Proba- 
bilities. 

By  GEORGE  BOOLE,  D.C.L. 

Professor  of  Mathematics  in  the  Queen's  Univer- 
sity, Ireland. 

1.  Differential  Equations. 

468  pp.  (1859).     Crown  8vo.  cloth.     14s. 

The  Author  has  endeavoured  in  this  treatise  to 
convey  as  complete  an  account  of  the  present  state 
of  knowledge  on  the  subject  of  Differential  Equa- 
tions as  was  consistent  with  the  idea  of  a  work 
intended,  primarily,  for  elementary  instruction. 
The  object  has  been  first  of  all  to  meet  the  wants 
of  those  who  had  no  previous  acquaintance  with 
the  subject,  and  also  not  quite  to  disappoint  others 
who  might  seek  for  more  advanced  information. 
The  earlier  sections  of  each  chapter  contain  that 
kind  of  matter  which  has  usually  been  thought 
suitable  for  the  beginner,  while  the  latter  ones  are 
devoted  either  to  an  account  of  recent  discovery, 
or  to  the  discussion  of  such  deeper  questions  of 
principle  as  are  likely  to  present  themselves  to  the 
reflective  student  in  connection  with  the  methods 
and  processes  of  his  previous  course. 


2.  The   Calculus  of  Finite 
Differences. 

248  pp.  (1860).    Crown  8vo.  cloth.    10s.  6d. 

In  this  work  particular  attention  has  been  paid 
to  the  connexion  of  the  methods  with  those  of  the 
Differential  Calculus — a  connexion  which  in  some 
instances  involves  far  more  than  a  merely  formal 
analogy.  The  work  is  in  some  measure  designed 
as  a  sequel  to  the  Author's  Treatise  on  Differential 
Hquations,  and  it  has  been  composed  on  the  same 
plan, 


By  EDWARD  JOHN  ROUTS,  M.A. 

Fellow  and  Assistant  Tutor  of  St.  Peter's  Co.llege, 
Cambridge. 

Dynamics  of  a  System  of 
Rigid  Bodies. 

With  numerous  Examples. 
336  pp.  (1860).    Crown  8vo.  cloth.    10s.  §d. 

CONTENTS  :  Chap.  I.  Of  Moments  of  Inertia.— 
II.  D'Alembert's  Principle.— III.  Motion  about 
a  Fixed  Axis.— IV.  Motion  in  Two  Dimensions.— 

V.  Motion  of  a  Rigid  Body  in  Three  Dimensions. — 

VI.  Motion  of  a  Flexible  String.— VII.  Motion  of 
a  System  of  Rigid  Bodies.— VIII.  Of  Impulsive 
Forces. — IX,  Miscellaneous  Examples. 

The  numerous  Examples  which  will  be  found  at 
the  end  of  each  chapter  have  been  chiefly  selected 
from  the  Examination  Papers  set  in  the  University 
and  Colleges  of  Cambridge  during  the  last  few 
years. 


By  N.  M.  FERRERS,  M.A. 

Fellow  and  Mathematical  Lecturer  of  Gonville  and 
Caius  College,  Cambridge. 

An  Elementary  Treatise  on 
Trilinear  Co-Ordinates. 

The  Method  of  Reciprocal  Polars,  and  the 
Theory  of  Projectiles. 

154  pp.  (1861).     Crown  8vo.  cloth.     6*.  6d. 

The  object  of  the  Author  in  writing  on  this 
subject  has  mainly  been  to  place  it  on  a  basis 
altogether  independent  of  the  ordinary  Cartesian 
System,  instead  of  regarding  it  as  only  a  special 
form  of  abridged  Notation.  A  short  chapter  on 
Determinants  has  been  introduced. 


By  W.  H.  DREW,  M.A. 

Second  Master  of  Blackheath  School. 

Geometrical     Treatise     on 
Conic  Sections. 

With  a  copious  Collection  of  Examples. 
Second  Edition.    Crown  8vo.  cloth.    4*.  6d. 

In  this  work  the  subject  of  Conic  Sections  has 
been  placed  before  the  student  in  such  a  form  that, 
it  is  hoped,  after  mastering  the  elements  of 
Euclid,  he  may  find  it  an  easy  and  interesting 
continuation  of  his  geometrical  studies.  With 
a  view  also  of  rendering  the  work  a  complete 
Manual  of  what  is  required  at  the  Universities, 
there  have  been  either  embodied  into  the  text,  or 
inserted  among  the  examples,  every  book-work 
question,  problem,  and  rider,  which  has  been 
proposed  in  the  Cambridge  examinations  up  to  the 
present  time. 


By  G.  H.  PUCKLE,   M.A. 
Principal    of  Windermere    College. 

Conic    Sections  and   Alge- 
braic Geometry. 

With  numerous  Easy  Examples 


Second  Edition.     264  pp.  (1856). 
Crown  8vo.     7s.  G<?. 

This  book  has  been  written  with  special  reference 
to  those  difficulties  and  misapprehensions  which 
commonly  beset  the  student  when  he  commences. 
With  this  object  in  view,  the  earlier  part  of  the 
subject  has  been  dwelt  on  at  length,  and  geo- 
metrical and  numerical  illustrations  of  the  analysis 
have  been  introduced.  The  Examples  appended 
to  each  section  are  mostly  of  an  elementary  de- 
scription. The  work  will,  it  is  hoped,  be  found  to 
contain  all  that  is  required  by  the  upper  classes  of 
schools  and  by  the  generality  of  students  at  the 
Universities. 


By    R.    D.    BEASLET,    M.A. 
Head  Master  of  Grantham  School. 

Plane  Trigonometry. 

AN  ELEMENTARY  TREATISE. 
With  a,  numerous   Collection   of  Examples. 

106  pp.   (1858),   strongly  bound  in  cloth. 
3s.  Qd. 

This  Treatise  is  specially  intended  for  use  in 
Schools.  The  choice  of  matter  has  been  chiefly 
guided  by  the  requirements  of  the  three  days' 
Examination  at  Cambridge,  with  the  exception' of 
proportional  parts  in  logarithms,  which  have  been 
omitted.  About  four  hundred  examples  have 
been  added,  mainly  collected  from  the  Examina- 
tion Papers  of  the  last  ten  years,  and  great  pains 
have  been  taken  to  exclude  from  the  body  of  the 
work  any  which  might  dishearten  a  beginner  by 
their  difliculty. 


By  J.  C.  SNOWBALL,  M.A. 
Late  Fellow  of  St.  John's  College,  Cambridge. 

Plane  and  Spherical  Trigo- 
nometry. 

With  the  Construction  and  Use  of  Tables  of 
Logarithms. 

Ninth  Edition.  240  pp.  (1857).  Crown  8vo. 
7*.  U. 

In  preparing  a  new  edition,  the  proofs  of  some 
of  the  more  important  propositions  have  been 
rendered  more  strict  and  general ;  and  a  con- 
siderable addition,  of  more  than  two  hundred 
examples,  taken  principally  from  the  questions  in 
the  Examinations  of  Colleges  and  the  University, 
has  been  made  to  the  collection  of  Examples  and 
Problems  for  practice. 


By  B.  DRAKE,  M.A. 
Late  Fellow  of  King's  College,  Cambridge. 

1.  Demosthenes  on  the 

Crown. 

With  English  Notes. 
Second    Edition.     To    which    is    prefixed 

-ZESCHINES    AGAINST    CTESIPHON.        "With 

English  Notes. 

287  pp.  (1860).     Fcap.  8vo.  cloth.    5s. 

The  first  edition  of  the  late  Mr.  Drake's  edition 
of  Demosthenes  de  Corona  having  met  with  con- 
siderable acceptance  in  various  Schools,  and  a  new 
edition  being  called  for,  the  Oration  of  JSschines 
against  Ctesiphon,  in  accordance  with  the  wishes 
of  many  teachers,  has  been  appended  with  useful 
notes  by  a  competent  scholar. 

2.  -ffischyli  Eumenides. 

With   English    Verse    Translation,    Copious 

Introduction,  and  Notes. 
"Mr.  Drake's  ability  as  a  critical  Scholar  is 
known  and  admitted.  In  the  edition  of  the 
Eumenides  before  us  we  meet  with  him  also  in  the 
capacity  of  a  Poet  and  Historical  Essayist.  The 
translation  is  flowing  and  melodious,  elegant  and 
scholarlike.  The  Greek  Text  is  well  printed  :  the 
notes  are  clear  and  useful." — GUARDIAN. 

By  J.  BROOK  SMITH,  M.A. 

St.  John's  College,  Cambridge. 

Arithmetic  in  Theory  and 
Practice. 

For  Advanced  Pupils. 
PART  I.     Crown  8vo.  cloth.      3s.  6d. 

This  work  forms  the  first  part  of  a  Treatise 
on  Arithmetic,  in  which  the  Author  has  en- 
deavoured, from  very  simple  principles,  to  explain 
in  a  full  and  satisfactory  manner  all  the  important 
processes  in  that  subject. 

The  proofs  have  in  all  cases  been  given  in  a  form 
entirely  arithmetical:  for  the  author  does  not 
think  that  recourse  ought  to  be  had  to  Algebra 


long  and  perplexing. 

At  the  end  of  every  chapter  several  examples 
have  been  worked  out  at  length,  in  which  the  best 
practical  methods  of  operation  have  been  carefully 
pointed  out. 

By  the  Rev.  G.  F.  CHILD E,  M.A. 
Mathematical  Professor  in  the  South  African  Coll. 

Singular  Properties  of  the 
Ellipsoid 

And  Associated  Surfaces  of  the  Nth  Degree. 

152  pp.  (1861).     8vo.  boards.     !Qs.*6d. 

As  the  title  of  this  volume  indicates,  its  object  is 
to  develope  peculiarities  in  the  Ellipsoid;  and 
further,  to  establish  analogous  properties  in  un- 
limited congeneric  series  of  which  this  remarkable 
surface  is  a  constituent. 


By  C.   W.  UNDERWOOD,  M.A. 

Vice-Principal    of   the    Collegiate   Institution 
Liverpool. 

A  Short  Manual  of 
Arithmetic. 

Fcp.  8vo.  96pp.  (1860).  limp  cloth.  2*.  Qd. 

The  object  aimed  at  by  the  Compiler  of  this 
Manual  is  to  bring  before  junior  students  so  much 
of  the  Theory  of  Arithmetic  as  may  be  fairly  ex- 
pected of  them,  and  to  present  it  in  such  a  form 
that  the  study  of  the  Science  may  become  to  some 
extent  a  mental  training.  It  is  rather  a  Grammar 
of  Arithmetic  than  a  treatise  on  that  subject,  and 
should  for  the  most  part  be  committed  to  memory. 
It  will  be  found  well  adapted  for  vivd  voce  ex- 
amination, and  enable  candidates  to  prepare  them- 
selves for  the  Local  University  Examinations. 
The  Definitions  are  briefly  and  carefully  worded. 
Each  rule  is  stated  so  as  to  include  the  proof  of  it 
where  this  was  possible. 


Senate -House  Mathe- 
matical Problems. 

With  Solutions. 

1848-51.    By  FERRERS  and  JACKSON.  8vo.   15s.  6d. 
1848-51.     (RIDERS).    By  JAMESON.    8vo.    7*.  6d. 
1854.    By  WALTON  and  MACKENZIE.   8vo.    10s.  6d. 
1857.    By  CAMPION  and  WALTON.     8vo.    8s.  6d. 
1860.    By  ROUTH  and  WATSON.  Crown  8vo.  7s.  6d. 

The  above  books  contain  Problems  and  Ex- 
amples which  have  been  set  in  the  Cambridge 
Senate-house  Examinations  at  various  periods 
during  the  last  twelve  years,  together  with  Solu- 
tions of  the  same.  The  Solutions  are  in  all  cases 
given  by  the  Examiners  themselves  or  under  their 
sanction. 


By  H.  A.  MORGAN,  M.A. 

Fellow  of  Jesus  College,  Cambridge. 

A  Collection  of  Mathemat- 
ical Problems  and  Ex- 
amples, with  Answers. 

190  pp.  (1858).     Crown  8vo.     6s.  6d. 

This  book  contains  a  number  of  problems, 
chiefly  elementary,  in  the  Mathematical  subjects 
usually  read  at  Cambridge.  They  have  been 
selected  from  the  papers  set  during  late  years  at 
Jesus  College.  Very  few  of  them  are  to  be  met 
with  in  other  collections,  and  by  far  the  larger 
number  are  due  to  some  of  the  most  distinguished 
Mathematicians  in  the  University. 


By  JOHN  E.  £.  MA  TOE,  M.A. 

Fellow  and  Classical  Lecturer  of  St.  John's  College, 
Cambridge. 

1.    Juvenal. 

With  English  Notes. 
464  pp.  (1854).  Crown  8vo.  cloth.    10*.  6d. 

"  A  School  edition  of  Juvenal,  which,  for  really 
ripe  scholarship,  extensive  acquaintance  with 
Latin  literature,  and  familiar  knowledge  with 
Continental  criticism,  ancient  and  modern,  is  un- 
surpassed, we  do  not  say  among  English  School- 
books,  but  among  English  editions  generally." — 
EDINBURGH  REVIEW. 


2.  Cicero's  Second  Philippic 

With  English  Notes. 
168  pp.  (1861).     Fcp.  8vo.  cloth.      5*. 

The  Text  is  that  of  Halm's  2nd  edition,  (Leipsig, 
Weidmann,  1858),  with  some  corrections  from 
Madvig's  4th  Edition  (Copenhagen,  1858) .  Halm's 
Introduction  has  been  closely  translated,  with 
some  additions.  His  notes  have  been  curtailed, 
omitted,  or  enlarged,  at  discretion;  passages  to 
which  he  gives  a  bare  reference,  are  for  the  most 
part  printed  at  length  ;  for  the  Greek  extracts  an 
English  version  has  been  substituted.  A  large 
body  of  notes,  chiefly  grammatical  and  historical, 
has  been  added  from  various  sources.  A  list  of 
books  useful  to  the  student  of  Cicero,  a  copious 
Argument,  and  an  Index  to  the  introduction  and 
notes,  complete  the  book. 


The  Cambridge  Year  Book 

AND  UNIVERSITY  ALMANACK 
For  1862. 

Crown  8vo.  228  pp.  price  Is.  6d. 

The  specific  features  of  this  annual  publication 
will  be  obvious  at  a  glance,  and  its  value  to 
teachers  engaged  in  preparing  students  for,  and  to 
parents  who  are  sending  their  sons  to,  the  Uni- 
versity, and  to  the  public  generally  will  be  clear. 

1.  The  whole  mode  of  proceeding  in  entering 
a  student  at  the  University  and  at  any  particular 


2.  The  course  of  the  studies  as  regulated  by  the 
University  examinations,  the  manner  of  these  ex- 
aminations, and  the  specific  subjects  and  times  for 
the  year  1862  are  given. 

3.  A  complete  account  of  all  Scholarships  and 
Exhibitions  at  the  several  Colleges,  their  value,  and 
the  means  by  which  they  are  gained. 

4.  A  brief  summary  of  all  Graces  of  the  Senate, 
Degrees  conferred  during  the  year  1861,  and  Uni- 
versity news  generally  ate  given. 

5.  The  Regulations  for  the  LOCAL  EXAMINATION 
of  those  who  are  not  members  of  the  University, 
to  be  held  this  year,  with  the  names  of  the  books 
oh  which  the  Examination  will  be  based,  and  the 
date  on  which  the  Examination  will  be  held. 


%  C.  MERIVALE,  B.D. 
Author  of  "A  History  of  Rome,"  &c. 

Sallust. 


With  English  Notes. 

Second  Edition.      172  pp.   (1858). 
8vo.     4s.  6d. 


Fcap. 


"  This  School  edition  of  Sallust  is  precisely  what 
the  School  edition  of  a  Latin  author  ought  to  be. 
No  useless  words  are  spent  in  it,  and  no  words 
that  could  be  of  use  are  spared.  The  text  has 
been  carefully  collated  with  the  best  editions.  It 
is  printed  in  a  large  bold  type,  which  manifests 
a  just  regard  for  the  young  eyes  that  are  to  work 
upon  it  :  under  the  text  there  flows  through  every 
page  a  full  current  of  extremely  well-selected 
annotations."—  THE  EXAMINER. 

The  "  CATILINA"  and  "  JUGUB.THA"  may  be  had 
separately,  price  2s.  6d.  each,  bound  in  cloth. 


Cambridge  University 

Examination  Papers, 

1860-1. 

A  Collection  of  all  the  Papers  Set  at  the 
Examinations  for  the  Degrees,  the  various 
Triposes,  and  the  Theological  Certificates 
in  the  University,  with  List  of  Candidates 
Examined  and  of  those  Approved,  and  an 
Index  to  the  Subjects. 

Previous  Examination,  1860— Previous  Exam- 
ination, 1861— B.A.  Degree  Examination,  1860— 
B.A.  Degree  Examination,  Jan.  and  May,  1861— 
Bachelor  of  Laws  Examination,  1860 — Bachelor  of 
Laws  Examination,  1861— Bachelor  of  Medicine 
Examinations,  1861— Classical  Tripos,  1861— Moral 
Sciences  Tripos,  1861— Natural  Sciences  Tripos, 
1861— Smith's  Prizes,  1861— Chancellor's  Medals 
for  Legal  and  Classical  Studies,  1861. 


By  J.  HERBERT  LATHAM,  M.A. 
Civil  Engineer. 

The  Construction  of 
"Wrought-Iron   Bridges. 

Embracing  the  Practical  Application  of  the 
Principles  of  Mechanics  to  Wrought-Iron 
Girder-  Work. 

"The  great  merit  of  this  book  is  that  it  deals 
with  practice  more  than  theory.  All  the  calcula- 
tions in  the  book  connected  with  the  strength  of 
girders  are  based  upon  their  actual  application 
which  abounds  in  practical  investigations  into 
girder-work  in  all  its  bearings,  and  will  be 
welcomed  as  one  of  the  most  valuable  contributions 
yet  made  to  this  important  branch  of  engineering." 
— ATHEN^UM. 


By  J.  WRIGHT,  M.A. 
Head  Master  of  Sutton  Coldfield  School. 

1.    Help  to  Latin  Grammar. 

With    Easy    Exercises,    both    English   and 

Latin,  Questions  and  Vocabulary. 
175  pp.  (1855).  Crown  8vo.  cloth.  4s.  6d. 
"  This  book  aims  at  helping  the  learner  to  over- 
step the  threshold  difficulties  of  the  Latin  Gram- 
mar ;  and  never  was  there  a  better  a:;d  offered 
alike  to  teacher  and  scholar  in  that  arduous  pass. 
The  style  is  at  once  familiar  and  strikingly  simple 
and  lucid ;  and  the  explanations  precisely  hit  the 
difficulties,  and  thoroughly  explain  them.  It  will 
also  much  facilitate  the  acquirement  of  English 
Grammar." — ENGLISH  JOURNAL  OF  EDUCATION. 


3.    Hellenica. 

A  FIRST  GREEK  READING  BOOK. 

From    Diodorus    and    Thucydides.        With 
Vocabulary. 

Second  Edition.      150  pp.   (1851).     Fcap. 
8vo.  cloth.     3s.  6d. 

In  the  last  twenty  chapters  of  this  volume, 
Thucydides  sketches  the  rise  and  progress  of  the 
Athenian  Empire  in  so  clear  a  style  and  in  such 
simple  language,  that  the  author  doubts  whether 
any  easier  or  more  instructive  passages  can  be  se- 
lected for  the  use  of  the  pupil  who  is  commencing 
Greek. 


4.  Vocabulary  and  Ex- 

ercise,s  on  "  The  Seven 

Kings  of  Rome." 

94  pp.  (1857).     Crown  8vo.  cloth.      2s.  6d. 

The  Vocabulary  is  published  apart  from  the 
Text  in  order  to  suit  the  views  of  those  who 
prefer  their  pupils  to  consult  a  general  dictionary, 
but  it  may  also  be  had  bound  together  with  the 
"  SEVEN  KINGS  OF  ROME,"  if  preferred.  As  the 
aim  of  the  Text  is  to  teach  the  elements  of  grammar, 
so  the  Exercises  are  intended  to  test  the  Pupil's 
knowledge  of  grammar.  Indeed  there  is  hardly 
an  ordinary  Latin  construction  which  is  not  illus- 
trated in  the  text,  explained  in  the  notes,  and 
proved  in  the  Exercises. 


2.  The  Seven  Kings  of  Rome 

A  First  Latin  Reading  Book,  abridged 
from  Livy,  by  the  omission  of  difficult 
passages,  with  Notes  and  Index. 

Second  Edition,      138  pp.   (1857).     Fcap. 
8vo.  cloth.     3s. 

This  work  is  intended  to  supply  the  pupil  with 
an  easy  Construing-book,  which  may,  at  the  same 
time,  be  made  the  vehicle  for  instructing  him  in 
the  rules  of  grammar  and  principles  of  compo- 
sition. These  branches  of  the  study  of  Latin  seem 
to  the  author  to  have  hitherto  been  kept  too  much 
apart.  Boys  have  construed  their  Delectus,  or  Eu- 
tropius,  or .  Nepos,  and  have  gone  elsewhere  for 
their  grammatical  exercises.  Nor  can  this  be  won- 
dered at.  An  educated  man  must  feel  positively 
ashamed  of  taking  his  pupils  away  from  our  good 
English  authors,  and  setting  before  him  instead 
a  Delectus  or  Eutropius.  He  therefore  skims  over 
them  as  lightly,  and  escapes  from  them  as  quickly 
as  possible,  and  has  recourse  for  his  composition 
lesson  to  one  of  the  many  Exercise-books  which 
swarm  from  our  educational  press.  To  remedy 
these  evils  this  book  has  been"  published.  Here 
Livy  tells  his  own  pleasant  stories  in  his  own 
pleasant  words.  What  is  omitted,  is  that  which 
no  one  can  wish  a  beginner  to  learn,  and  which 
may  be  better  learnt  elsewhere.  Let  Livy  be  the 
master  to  teach  a  boy  Latin,  not  some  English 
collector  of  sentences,  and  he  will  not  be  found 
a  dull  one. 


Elementary  Statics. 

By  the  Rev.  GEORGE  RA  WLINSON, 

Professor  of  Applied  Sciences,  Elphinstone  Coll., 
Bombay. 

Edited  by  the  Rev.  E.  STURGES,  M.A. 

Rector  of  Kencott,  Oxfordshire. 
(150  pp.)  1860.     Crown  8vo.  cloth.    4s.  Qd. 

This  work  is  published  under  the  authority  of 
H.  M.  Secretary  of  State  for  India  for  use  in  the 
Government  Schools  and  Colleges  in  India. 


By  P.  FROST,  Jun.,  M.A. 
Late  Fellow  of  St.  John's  College,  Cambridge. 

Thucydides.     Book  VI. 

With  English  Notes,  Map  and  Index. 
8vo.  cloth.     Is.  Qd. 

It  has  been  attempted  in  this  work  to  facilitate 
the  attainment  of  accuracy  in  translation.  With 
this  end  in  view  the  Text  has  been  treated  gram- 
matically. 


By  EDWARD  TRUING,  M.A. 
Head  Master  of  Uppingham  Grammar  School. 

1.    Elements    of    Grammar 
Taught  in  English. 

With  Questions. 

Third  Edition.      136  pp.    (1860).      Demy 
18mo.     2s. 

2.    The  Child's  English 
Grammar. 

New    Editton.      86    pp.    (1859).      Demy 
18mo.     1*. 

The  Author's  effort  in  these  two  books  has  been 
to  point  out  the  broad,  beaten,  every-day  path, 
carefully  avoiding  digressions  into  the  byeways 
and  eccentricities  of  language.  This  Work  took 
its  rise  from  questionings  in  National  Schools,  and 
the  whole  of  the  first  part  is  merely  the  writing 
out  in  order  the  answers  to  questions  which  have 
been  used  already  with  success.  The  study  of 
Grammar  in  English  has  been  much  neglected, 
nay  by  some  put  on  one  side  as  an  impossibility. 
There  was  perhaps  much  ground  for  this  opinion, 
in  the  medley  of  arbitrary  rules  thrown  before  the 
student,  which  applied  indeed  to  a  certain  number 
of  instances,  but  would  not  work  at  all  in  many 
others,  as  must  always  be  the  case  when  principles 
are  not  put  forward  in  a  language  full  of  ambigui- 
ties. The  present  work  does  not,  therefore,  pre- 
tend to  be  a  compendium  of  idioms,  or  a  philo- 
logical treatise,  but  a  Grammar.  Or  in  other 
words,  its  intention  is  to  teach  the  learner  how  to 
speak  and  write  correctly,  and  to  understand  and 
explain  the  speech  and  writings  of  others.  Its 
success,  not  only  in  National  Schools,  from  practi- 
cal work  in  which  it  took  its  rise,  but  also  in 
classical  schools,  is  full  of  encouragement. 


3.  A  First  Latin  Construing 
Book. 

104  pp.  (1855).     Fcap.  8vo.     2s.  6d. 

This  Construing  Book  is  drawn  up  on  the  same 
sort  of  graduated  scale  as  the  Author's  English 
Grammar.  Passages  out  of  the  best  Latin  Poets 
are  gradually  built  up  into  their  perfect  shape. 
The  few  words  altered,  or  inserted  as  the  passages 
go  on,  are  printed  in  Italics.  It  is  hoped  by  this 
plan  that  the  learner,  whilst  acquiring  the  rudi- 
ments of  language,  may  store  his  mind  with  good 
poetry  and  a  good  vocabulary. 

4.   School  Songs. 

A   COLLECTION  OP  SONGS  FOR  SCHOOLS. 

WITH   THE  MUSIC   ARRANGED   FOR  FOUR 
VOICES. 

Edited  by  Rev.  E.  THRING  %  H.  SICCIUS. 

Music  Size.     7*.  §d. 


By  C.  J.   VAUGHAN,  D.D. 
Head  Master  of  Harrow  School. 

St.   Paul's   Epistle   to   the 
Ho  mans. 

The  Greek  Text  with  English  Notes. 
Second  Edition.   Crown  Svo.cl.    (1861).  5s. 

By  dedicating  this  work  to  his  elder  Pupils  at 
Harrow,  the  Author  hopes  that  he  sufficiently 
indicates  what  is  and  what  is  not  to  be  looked  for 
in  it.  He  desires  to  record  his  impression,  de- 
rived from  the  experience  of  many  years,  that  the 
Epistles  of  the  New  Testament,  no  less  than  the 
Gospels,  are  capable  of  furnishing  useful  and 
solid  instruction  to  the  highest  classes  of  our 
Public  Schools.  If  they  are  taught  accurately,  not 
controversially ;  positively,  not  negatively ;  au- 
thoratively,  yet  not  dogmatically ;  taught  with 
close  and  constant  reference  to  their  literal  mean- 
ing, to  the  connexion  of  their  parts,  to  the  sequence 
of  their  argument,  as  well  as  to  their  moral  and 
spiritual  instruction  ;  they  will  interest,  they  will 
inform,  they  will  elevate ;  they  will  inspire  a  reve- 
rence for  Scripture  never  to  be  discarded,  they 
will  awaken  a  desire  to  drink  more  deeply  of  the 
Word  of  God,  certain  hereafter  to  be  gratified  and 
fulfilled. 


Notes  for  Lectures  on  Confirma- 
tion: "With  Suitable  Prayers.  By 
C.  J.  VAUGHAN,  D.D.  Fourth  Edition. 
70  pp.  (1862).  Fcp.  8vo.  1*.  6d. 

The  Church  Catechism  Illustrated 
and  Explained.  By  ARTHUR 
RAMSAY,  M.A.  18mo.  cloth. 

Hand-Book  to  Butler's  Analogy. 
By  C.  A.  SWAINSON,  M.A.  55  pp. 
(1856).  Crown  8vo.  Is.  6d. 

History  of  the  Christian  Church 
during  the  First  Three  Cen- 
turies, and  the  Reformation  in 
England.  By  W.  SIMPSON,  M.A. 
307  pp.  (1857).  Fcp.  8vo.  cloth.  5s. 

Analysis  of  Paley's  Evidences  of 
Christianity.  By  CHARLES  H. 
CROSSE,  M.A.  115  pp.  (1855).  18mo. 
3s.  6rf. 
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UNIFORMLY  PRINTED    AND    BOUND. 

This  Series  of  Theological  Manuals  has  been  published  with  the  aim  of 
supplying  Books  concise,  comprehensive,  and  accurate,  convenient  for 
the  Student,  and  yet  interesting  to  the  general  reader. 


1.  History  of  the  Christian  Church 
during  the  Middle  Ages,  By 
ARCHDEACON  HARDWICK.  Second 
Edition.  482  pp.  (1861).  With  Maps. 
Crown  8vo.  cloth.  10s.  6d. 


This  Volume  claims  to  be  regarded  as  an  integral 
and  independent  treatise  on  the  Mediaeval  Church. 
The  History  commences  with  the  time  of  Gregory 
the  Great,  because  it  is  admitted  on  all  hands  that 
his  pontificate  became  a  turning-point,  not  only  in 
the  fortunes  of  the  Western  tribes  and  nations,  but 
of  Christendom  at  large.  A  kindred  reason  has 
suggested  the  propriety  of  pausing  at  the  year 
1520,  —  the  year  when  Luther,  having  been  ex- 
truded from  those  Churches  that  adhered  to  the 
Communion  of  the  Pope,  established  a  provisional 
form  of  government  and  opened  a  fresh  era  in  the 
history  of  Europe. 


2.  History  of  the  Christian  Church 

during  the  Reformation.  By 
ARCHDEACON  HARDWICK.  459  pp. 
(1856).  Crown  8vo.  cloth.  10s.  6d. 

This  Work  forms  a  Sequel  to  the  Author's  Book 
on  The  Middle  Ages.  The  Author's  wish  has  been 
to  give  the  reader  a  trustworthy  version  of  those 
stirring  incidents  which  mark  the  Reformation 
period. 

3.  History  of  the  Book  of  Common 
Prayer.  "With  a  Rationale  of  its  Offi- 
ces.    By  FRANCIS  PROCTER,  M.A. 

Fifth  Edition.  464  pp.  (1860).  Crown 
8vo.  cloth.  10s.  6d. 


The  Subject  of  this  Book  has  been  already 
treated  by  numerous  writers  of  distinction.  When 
the  present  series  of  Manuals  was  projected,  it  did 
not  appear  that  any  one  of  the  existing  volumes 
taken  singly  was  available  for  the  desired  object. 


In  the  course  of  the  last  twenty  years  the  whole 
question  of  liturgical  knowledge  has  been  reopened 
with  great  learning  and  accurate  research,  and  it  is 
mainly  with  the  view  of  epitomizing  their  extensive 
publications,  and  correcting  by  their  help  the 
errors  and  misconceptions  which  had  obtained 
currency,  that  the  present  volume  has  been  put 
together. 


4.  History  of  the  Canon  of  the  New 
Testament  during  the  First 
Four  Centuries.  By  BROOKE 
FOSS  WESTCOTT,  M.A.  594  pp. 
(1855).  Crown  8vo.  cloth.  12*.  6d. 

The  Author  has  endeavoured  to  connect  the 
histoiy  of  the  New  Testament  Canon  with  the 
growth  and  consolidation  of  the  Catholic  Church, 
and  to  point  out  the  relation  existing  between  the 
amount  of  evidence  for  the  authenticity  of  its  com- 
ponent parts  and  the  whole  mass  of  Christian 
literature.  Such  a  method  of  inquiry  will  convey 
both  the  truest  notion  of  the  connexion  of  the 
written  Word  with  the  living  Body  of  Christ,  and 
the  surest  conviction  of  its  divine  authority. 


5.  Introduction  to  the  Study  of  the 
GOSPELS.  By  BROOKE  FOSS 
WESTCOTT,  M.A.  458  pp.  (1860). 
Crown  8vo.  cloth.  10s.  6d. 

The  title  of  this  book  will  explain  the  chief  aim 
which  the  Author  had  in  view.  It  is  intended  to 
be  an  Introduction  to  the  Study  of  the  Gospels. 
The  Author  has  therefore  confired  himself  in  many 
cases  to  the  mere  indication  of  lines  of  thought  and 
inquiry  from  the  conviction  that  truth  is  felt  to  be 
more  precious  in  proportion  as  it  is  opened  to  us  by 
our  own  work.  In  a  subject  which  involves  so 
vast  a  literature  much  must  have  been  overlooked ; 
but  the  Author  has  made  it  a  point  at  least  to  study 
the  researches  of  the  great  writers,  and  consciously 
to  neglect  none. 


BOOKS  SUITABLE  FOR  SCHOOL  PRIZES. 

KEPT  IN  VAKIOUS  BINDINGS  BY  THE  PUBLISHERS. 
THE  GOLDEN  TREASURY  OF  THE  BEST  SONGS  AND  LYRICAL 

POEMS  IN  THE  ENGLISH  LANGUAGE.    Cloth,  4*.  Gd.,  morocco  plain,  7s.  Gd.,  extra, 
10s.  6d. 

THE  CHILDREN'S  GARLAND  FROM  THE  BEST  POETS.     Cloth, 

4s.  Gd.  morocco  plain,  7s.  Gd.,  extra,  10s.  Gd. 

BUNYAN'S  PILGRIM'S  PROGRESS.     Cloth,  4s.  6d.,  morocco  plain, 

7s.  Gd.,  extra,  10s.  6d. 

TOM  BROWN'S  SCHOOL  DAYS.     BY  AN  OLD  BOY.     Seventh  Edition. 

Fcap.  8vo.  cloth,  5s.,  calf,  9s.,  morocco,  11s. 

THE  HEROES;  or  GREEK  FAIRY  TALES.     By  CHARLES  KINGSLEY, 

Rector  of  Eversley.    Second  Edition,  with  EIGHT  ILLUSTRATIONS.    Imperial  16mo.  cloth, 
gilt  leaves,  5s.,  calf,  9s.  Gd.,  morocco,  12s. 

DAY1D,    KING    OF    ISRAEL;    Readings    for   the    Young.       With 

Six  ILLUSTRATIONS.    By  J.  WEIGHT,  M.A.    Imp.  16mo.  cloth,  5s.,  calf,  9s.  Gd.,  mor.  12s. 

LITTLE   ESTELLA   AND  OTHER  FAIRY   TALES.      Imp.    16mo. 

cloth,  5s.,  calf,  9s.  Gd.,  mor.  12s. 

PAULI'S    PICTURES    OF    OLD    ENGLAND.      Crown    8vo.    cloth, 

8s.  Gd.,  calf,  13s.  6d.,  mor.  16s.  Gd. 

DAYS   OF   OLD:       STORIES  FROM   OLD   ENGLISH   HISTORY 

of  the  Druids,  Anglo-Saxons,   and  the  Crusaders.     By  the  Author  of  RUTH  AND  HER 
FRIENDS.    Cloth,  5s.,  calf,  9s.  6d.,  mor.  12s. 

RUTH  AND  HER  FRIENDS.       A  Story  for  Girls.       Third  Edition. 

Cloth,  5s.,  calf,  9s.  Gd.,  mor.  12s. 

OTJR  YEAR:    A  Child's  Book  in  Prose   and  Verse.      By  the   Author 

of  "John  Halifax."    Numerous  Illustrations.     Royal  16mo.  cloth,  gilt  leaves,  5s.,  calf, 
9s.  Gd.,  mor.  12s.. 

EARLY  EGYPTIAN  HISTORY  FOR  THE  YOUNG.     With  Descrip- 

tion  of  the  Tombs  and  Monuments.    Fcp.  8vo.  cloth,  5s.,  calf,  9s.  Gd.,  mor.  12s. 

WESTWARD  HO !     THE  ADVENTURES  OF  SIR  AMYAS  LEIGH 

in  the  Reign  of  Elizabeth.     Third  Edition.     By  CHARLES  KINGSLEY.    Crown  8vo.  cloth,  6s. 
calf,  11s.,  mor.  14s. 

TWO  YEARS  AGO.     By  CHARLES  KINGSLEY.     Third  Edition.     Crown 

8vo.  cloth,  6s.,  calf,  11s.,  mor.  14s. 

THE    RECOLLECTIONS    OF    GEOFFRY   HAMLYN.      By    HENRY 

KINGSLEY.    SECOND  EDITION.    Crown  8vo.  cloth,  6s.,  calf,  11s.,  mor.  14s. 

GLAUCUS ;  or,  WONDERS  OF  THE  SHORE.     By  CHARLES  KINGSLEY. 

ILLUSTRATED  EDITION,  containing  Coloured  Illustrations  of  the  objects  mentioned  in  the 
Work.    Imp.  16mo.  cloth,  gilt  leaves,  7s.  Gd.,  calf,  12s..  mor.  14s.  Gd. 

ESSAYS,  CHIEFLY  ON  ENGLISH  POETS.     By  DAVID  MASSON,  M.A. 

8vo.  cloth,  12s.  Gd.,  calf,  18s.,  mor.  22s.  Gd. 

ARCHER   BUTLER'S    HISTORY    OF    ANCIENT    PHILOSOPHY. 

2  vols.  8vo.  cloth,  II.  5s.,  calf,  11.  16s.,  mor.  21.  5s. 

HISTORY  OF  THE  CANON  OF  THE  NEW  TESTAMENT.  By 

B.  F.  WESTCOTT,  M.A.  Third  Edition.    Crown  8vo.  cloth,  12s.  Gd.t  calf,  17s.  Gd.,  mor.  II.  Gs. 

INTRODUCTION    TO    THE    STUDY    OF    THE    GOSPELS.      By 

B.  F.  WESTCOTT,  M.A.    Crown  8vo.  cloth,  10s.  Gd.,  calf,  15s.  Gd.,  mor.  18s.  Gd. 

GEORGE  BRIMLEY'S  ESSAYS.      Second  Edition.     Fcp.  8vo.  cloth, 

5s.,  calf,  9s.,  mor.  11s. 

THE    PLATONIC   DIALOGUES   FOR   ENGLISH   READERS.     By 

W.  WHEWELL,  D.D.    3  vols.  Fcp.  Svo.  cloth.    21s.  Gd.,  calf,  35s.,  mor.  42s.  Gd. 
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MACMILLAN   &  CO.'S 

ELEMENTARY   SCHOOL    CLASS   BOOKS. 
Euclid  for  Colleges  and  Schools. 

By    I.    TODHUNTER,    M.A.,    F.R.S.,    Fellow    and    Principal   Mathematical 

Lecturer  of  St.   John's   College,    Cambridge.  [Nearly  Ready. 

***  This  work,   which  will  form   a  volume  of  Hacmillan  and   Co.'s  ELEMENTARY 

SCHOOL  CLASS  BOOKS,  will  be  handsomely  printed  in  I8mo.  and  all  the  volumes 

of  the  series  will  be  published  at  a  low  price,  to  ensure  an  extensive  sale  in  the 

Schools  of  the  United  Kingdom  and  the  Colonies. 

The  following  will  form  early  volumes  in  this  Series  : 

An  Elementary  History  of  the  Prayer-Book. 

By  FRANCIS   PROCTER,   M.A.,   Vicar  of  Witton,   Norfolk,   late   Fellow  of 
St.  Catharine's  College,  Camhridge.  [Nearly  Ready. 

The  School  Class  Book  of  Arithmetic. 

By  BARNARD  SMITH,   M.A.,    Fellow  of  St.   Peter's  College,   Camhridge. 


Elementary  Latin  Grammar. 


By  H.  J.  ROBY,  M.A.,  Assistant  Master  in  Dulwich  College,  formerly  Fellow 
of  St.  John's  College. 

The  Bible   Word-Book. 

A  Glossary  of  old  English  Bible  Words,  with  Illustrations.  By  J.  EASTWOOD, 
M.A.,  St.  John's  College,  Cambridge,  and  W.  ALDIS  WRIGHT,  M.A.,  Trinity 
College,  Cambridge. 

*#*  Other  Volumes  will  be  announced  in  due  course. 

FORTHCOMING    BOOKS. 

An  Elementary  Treatise  on  the  Planetary  Theory. 

By  C.  H.  H.  CHEYNE,  B.A.,  Scholar  of  St.  John's  College.         [In  the  Press. 

An  Elementary  Treatise  on  Natural  Philosophy. 

By  WILLIAM  THOMSON,  LL.D.,  F.R.S.,  late  Fellow  of  St.  Peter's  College, 
Cambridge,  Professor  of  Natural  Philosophy  in  the  University  of  Glasgow ;  and 
PETER  GUTHRIE  TAIT,  M.A.,  late  Fellow  of  St.  Peter's  College,  Cambridge, 
Professor  of  Natural  Philosophy  in  the  University  of  Edinburgh.  With  numerous 
Illustrations.  [Preparing. 

An  Elementary  Treatise  on  Quaternions. 

With  numerous  Examples.  By  P.  G.  TAIT,  M.A.,  Professor  of  Natural  Phi- 
losophy in  the  University  of  Edinburgh.  [Preparing. 

A   Treatise  on  Geometry  of  Three  Dimensions. 

By  PERCIVAL  FROST,  M.A.,  St.  John's  CoUege,  and  JOSEPH  WOLSTEN- 
HOLME,  M.A.,  Christ's  College,  Cambridge.  [In  the  Press. 

\*  The  First  Portion  has  been  issued  for  the  convenience  of  Cambridge  Students. 

First  Book  of  Algebra.     For  Schools. 

By  J.  C.  W.  ELLIS,  M.A.,  and  P.  M.  CLARK,  M.A.,  Sidney  Sussex  College, 
Cambridge.  [Preparing. 

The  New  Testament  in  the   Original  Greek. 

Text  revised  by  B.  F.  WESTCOTT,  M.A.,  and  F.  J.  HORT,  M.A.,  formerly 
Fellows  of  Trinity  College. 

Aristotelis  de  Rhetorica. 

\Vitl.  "N"r>fpa  n-nrl  Introduction.     Bv  E.  M.  COPE.  M.A..  Fellow  and  Assistant 
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